Specific capture and whole-genome sequencing of viruses from clinical samples

Detection of enterovirus RNA in pancreas and lymphoid tissues of organ donors with type 1 diabetes

Aims/hypothesis

The nPOD-Virus group collaboratively applied innovative technologies to detect and sequence viral RNA in pancreas and other tissues from organ donors with type 1 diabetes. These analyses involved the largest number of pancreas samples collected to date.

Methods

We analysed pancreas, spleen, pancreatic lymph nodes, and duodenum samples from the following donor groups: a) donors with type 1 diabetes (n=71), with (n=35) or without (n=36) insulin-containing islets, (b) donors with single or double islet autoantibody positivity without diabetes (n=22) and c) autoantibody-negative donors without diabetes (control donors) (n=74). Five research laboratories participated in this collaborative effort using approaches for unbiased discovery of RNA viruses (two RNA-Seq platforms), targeted detection of Enterovirus A-D species using RT-PCR, and tests for virus growth in cell-culture.

Results

Direct RNA-Seq did not detect virus signal in pancreas samples, whereas RT-PCR detected enterovirus RNA confirmed by sequencing in low amounts in pancreas samples in three of the five donor groups, namely donors with type 1 diabetes with insulin-containing islets, 16% (5/32) donors being positive, donors with single islet autoantibody positivity with 53% (8/15) donors being positive, and non-diabetic donors with 8% (4/49) being enterovirus RNA positive. Detection of enterovirus RNA was significantly more frequent in single islet autoantibody-positive donors compared to donors with type 1 diabetes with insulin-deficient islets (p-value <0.001) and control donors (p-value 0.004). In some donors, pancreatic lymph nodes were also positive. RT-PCR detected enterovirus RNA also in spleen of a small number of donors and virus enrichment in susceptible cell lines before RT-PCR resulted in much higher rate in spleen positivity, particularly in donors with type 1 diabetes. Interestingly, the enterovirus strains detected did not cause a typical lytic infection, possibly reflecting their persistence-prone nature.

Conclusions/interpretation

This was the largest coordinated effort to examine the presence of enterovirus RNA in pancreas of organ donors with type 1 diabetes, using a multitude of assays. These findings are consistent with the notion that both the subjects with type 1 diabetes and those with islet autoantibodies may carry a low-grade enterovirus infection in the pancreas and lymphoid tissues.

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.

Persistent Low-Level Variants in a Subset of Viral Genes Are Highly Predictive of Poor Outcome in Immunocompromised Patients With Cytomegalovirus Infection

BACKGROUND

Human cytomegalovirus (HCMV) is the most common and serious opportunistic infection after solid organ and hematopoietic stem cell transplantation. In this study, we used whole-genome HCMV data to investigate viral factors associated with the clinical outcome.

METHODS

We sequenced HCMV samples from 16 immunocompromised pediatric patients with persistent viremia. Eight of the 16 patients died of complications due to HCMV infection. We also sequenced samples from 35 infected solid organ adult recipients, of whom 1 died with HCMV infection.

RESULTS

We showed that samples from both groups have fixed variants at resistance sites and mixed infections. Next-generation sequencing also revealed nonfixed variants at resistance sites in most of the patients who died (6/9). A machine learning approach identified 10 genes with nonfixed variants in these patients. These genes formed a viral signature that discriminated patients with HCMV infection who died from those who survived with high accuracy (area under the curve = 0.96). Lymphocyte numbers for a subset of patients showed no recovery posttransplant in the patients who died.

CONCLUSIONS

We hypothesize that the viral signature identified in this study may be a useful biomarker for poor response to antiviral drug treatment and indirectly for poor T-cell function, potentially identifying early those patients requiring nonpharmacological interventions.

Exploring viral diversity and metagenomics in livestock: insights into disease emergence and spillover risks in cattle

Cattle have a significant impact on human societies in terms of both economics and health. Viral infections pose a relevant problem as they directly or indirectly disrupt the balance within cattle populations. This has negative consequences at the economic level for producers and territories, and also jeopardizes human health through the transmission of zoonotic diseases that can escalate into outbreaks or pandemics. To establish prevention strategies and control measures at various levels (animal, farm, region, or global), it is crucial to identify the viral agents present in animals. Various techniques, including virus isolation, serological tests, and molecular techniques like PCR, are typically employed for this purpose. However, these techniques have two major drawbacks: they are ineffective for non-culturable viruses, and they only detect a small fraction of the viruses present. In contrast, metagenomics offers a promising approach by providing a comprehensive and unbiased analysis for detecting all viruses in a given sample. It has the potential to identify rare or novel infectious agents promptly and establish a baseline of healthy animals. Nevertheless, the routine application of viral metagenomics for epidemiological surveillance and diagnostics faces challenges related to socioeconomic variables, such as resource availability and space dedicated to metagenomics, as well as the lack of standardized protocols and resulting heterogeneity in presenting results. This review aims to provide an overview of the current knowledge and prospects for using viral metagenomics to detect and identify viruses in cattle raised for livestock, while discussing the epidemiological and clinical implications.

A Review of Probe-Based Enrichment Methods to Inform Plant Virus Diagnostics

Modern diagnostic techniques based on DNA sequence similarity are currently the gold standard for the detection of existing and emerging pathogens. Whilst individual assays are inexpensive to use, assay development is costly and carries risks of not being sensitive or specific enough to capture an increasingly diverse range of targets. Sequencing can provide the entire nucleic acid content of a sample and may be used to identify all pathogens present in the sample when the depth of coverage is sufficient. Targeted enrichment techniques have been used to increase sequence coverage and improve the sensitivity of detection within virus samples, specifically, to capture sequences for a range of different viruses or increase the number of reads from low-titre virus infections. Vertebrate viruses have been well characterised using in-solution hybridisation capture to target diverse virus families. The use of probes for genotyping and strain identification has been limited in plants, and uncertainty around sensitivity is an impediment to the development of a large-scale virus panel to use within regulatory settings and diagnostic pipelines. This review aims to compare significant studies that have used targeted enrichment of viruses to identify approaches to probe design and potential for use in plant virus detection and characterisation.

A novel fast hybrid capture sequencing method for high-efficiency common human coronavirus whole-genome acquisition

Rapid and accurate sequencing of the entire viral genome, coupled with continuous monitoring of genetic changes, is crucial for understanding the epidemiology of coronaviruses. We designed a novel method called micro target hybrid capture system (MT-Capture) to enable whole-genome sequencing in a timely manner. The novel design of probes used in target binding exhibits a unique and synergistic "hand-in-hand" conjugation effect. The entire hybrid capture process is within 2.5 hours, overcoming the time-consuming and complex operation characteristics of the traditional liquid-phase hybrid capture (T-Capture) system. By designing specific probes for these coronaviruses, MT-Capture effectively enriched isolated strains and 112 clinical samples of coronaviruses with cycle threshold values below 37. Compared to multiplex PCR sequencing, it does not require frequent primer updates and has higher compatibility. MT-Capture is highly sensitive and capable of tracking variants.IMPORTANCEMT-Capture is meticulously designed to enable the efficient acquisition of the target genome of the common human coronavirus. Coronavirus is a kind of virus that people are generally susceptible to and is epidemic and infectious, and it is the virus with the longest genome among known RNA viruses. Therefore, common human coronavirus samples are selected to evaluate the accuracy and sensitivity of MT-Capture. This method utilizes innovative probe designs optimized through probe conjugation techniques, greatly shortening the time and simplifying the handwork compared with traditional hybridization capture processes. Our results demonstrate that MT-Capture surpasses multiplex PCR in terms of sensitivity, exhibiting a thousandfold increase. Moreover, MT-Capture excels in the identification of mutation sites. This method not only is used to target the coronaviruses but also may be used to diagnose other diseases, including various infectious diseases, genetic diseases, or tumors.

Evaluation of a custom designed hybridisation assay for whole genome sequencing of human adenoviruses direct from clinical samples

BACKGROUND

Human Adenoviruses are a common cause of disease and can cause significant morbidity and mortality in immunocompromised patients. Nosocomial transmission events can occur with whole genome sequencing playing a crucial role. This study evaluates the performance of a custom designed SureSelectXT target enrichment assay based on 14 adenovirus genomes for sequencing direct from clinical samples.

METHODS

Modifications were made to the SureSelectXT low input protocol to enhance performance for viral targets. Consensus sequences were generated using an in-house designed three stage bioinformatics pipeline. We assessed, percentage of on target reads, average depth of coverage and percentage genome coverage to determine assay performance across a range of sample matrices.

RESULTS

Whole genome sequences were successfully generated for 91.6 % of samples assessed. Adenovirus DNA concentration was a good indicator of enrichment success. Highly specific enrichment was observed with only 6 % of samples showing < 50 % on target reads. Respiratory and faecal samples performed well where bloods showed higher levels of non-specific enrichment likely confounded by low adenovirus DNA concentrations. Protocol performance did not appear impacted by Adenovirus type or species.

CONCLUSION

Overall performance of this modified SureSelectXT protocol appears in line with previously published works although there are some confounding factors requiring further investigation. The use of a small RNA bait set has the potential to reduce associated costs which can be prohibitive.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Applying targeted gene hybridization capture to viruses with a focus to SARS-CoV-2

Although next-generation sequencing technologies are advancing rapidly, many research topics often require selective sequencing of genomic regions of interest. In addition, sequencing low-titre viruses is challenging, especially for coronaviruses, which are the largest RNA viruses. Prior to sequencing, enrichment of viral particles can help to significantly increase target sequence information as well as avoid large sequencing efforts and, consequently, can increase sensitivity and reduce sequencing costs. Targeting nucleic acids using capture by hybridization is another efficient method that can be performed by applying complementary probes (DNA or RNA baits) to directly enrich genetic information of interest while removing background non-target material. In studies where sequence capture by hybridization has been applied to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, most authors agree that this technique is useful to easily access sequence targets in complex samples. Furthermore, this approach allows for complete or near-complete sequencing of the viral genome, even in samples with low viral load or poor nucleic acid integrity. In addition, this strategy is highly efficient at discovering new variants by facilitating downstream investigations, such as phylogenetics, epidemiology, and evolution. Commercial kits, as well as in-house protocols, have been developed for enrichment of viral sequences. However, these kits have multiple variations in procedure, with differences in performance. This review compiles and describes studies in which hybridization capture has been applied to SARS-CoV-2 variant genomes.

Whole Genome Sequencing of Epizootic Hemorrhagic Disease Virus

Next-generation sequencing (NGS) technologies are continuously being developed and are becoming a more cost-effective tool for the characterization of viral genomes. Whole genome sequencing of segmented viruses, such as epizootic hemorrhagic disease virus (EHDV), provides insights into the molecular epidemiology as well as such viral evolutionary mechanisms as genetic reassortment. Here, we present a detailed method for obtaining full genome sequence data for EHDV using Illumina technology. The protocol includes details from RNA extraction and purification, the synthesis of cDNA, sequencing library preparation, to genome assembly.

Towards a post-pandemic future for global pathogen genome sequencing

Pathogen genome sequencing has become a routine part of our response to active outbreaks of infectious disease and should be an important part of our preparations for future epidemics. In this Essay, we discuss the innovations that have enabled routine pathogen genome sequencing, as well as how genome sequences can be used to understand and control the spread of infectious disease. We also explore the impact of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic on the field of pathogen genomics and outline the challenges we must address to further improve the utility of pathogen genome sequencing in the future.

Is Xenopus laevis introduction linked with Ranavirus incursion, persistence and spread in Chile?

Ranaviruses have been associated with amphibian, fish and reptile mortality events worldwide and with amphibian population declines in parts of Europe. Xenopus laevis is a widespread invasive amphibian species in Chile. Recently, Frog virus 3 (FV3), the type species of the Ranavirus genus, was detected in two wild populations of this frog near Santiago in Chile, however, the extent of ranavirus infection in this country remains unknown. To obtain more information about the origin of ranavirus in Chile, its distribution, species affected, and the role of invasive amphibians and freshwater fish in the epidemiology of ranavirus, a surveillance study comprising wild and farmed amphibians and wild fish over a large latitudinal gradient (2,500 km) was carried out in 2015-2017. In total, 1,752 amphibians and 496 fish were tested using a ranavirus-specific qPCR assay, and positive samples were analyzed for virus characterization through whole genome sequencing of viral DNA obtained from infected tissue. Ranavirus was detected at low viral loads in nine of 1,011 X. laevis from four populations in central Chile. No other amphibian or fish species tested were positive for ranavirus, suggesting ranavirus is not threatening native Chilean species yet. Phylogenetic analysis of partial ranavirus sequences showed 100% similarity with FV3. Our results show a restricted range of ranavirus infection in central Chile, coinciding with X. laevis presence, and suggest that FV3 may have entered the country through infected X. laevis, which appears to act as a competent reservoir host, and may contribute to the spread the virus locally as it invades new areas, and globally through the pet trade.

A large population sample of African HIV genomes from the 1980s reveals a reduction in subtype D over time associated with propensity for CXCR4 tropism

We present 109 near full-length HIV genomes amplified from blood serum samples obtained during early 1986 from across Uganda, which to our knowledge is the earliest and largest population sample from the initial phase of the HIV epidemic in Africa. Consensus sequences were made from paired-end Illumina reads with a target-capture approach to amplify HIV material following poor success with standard approaches. In comparisons with a smaller 'intermediate' genome dataset from 1998 to 1999 and a 'modern' genome dataset from 2007 to 2016, the proportion of subtype D was significantly higher initially, dropping from 67% (73/109), to 57% (26/46) to 17% (82/465) respectively (p < 0.0001). Subtype D has previously been shown to have a faster rate of disease progression than other subtypes in East African population studies, and to have a higher propensity to use the CXCR4 co-receptor ("X4 tropism"); associated with a decrease in time to AIDS. Here we find significant differences in predicted tropism between A1 and D subtypes in all three sample periods considered, which is particularly striking the 1986 sample: 66% (53/80) of subtype D env sequences were predicted to be X4 tropic compared with none of the 24 subtype A1. We also analysed the frequency of subtype in the envelope region of inter-subtype recombinants, and found that subtype A1 is over-represented in env, suggesting recombination and selection have acted to remove subtype D env from circulation. The reduction of subtype D frequency over three decades therefore appears to be a result of selective pressure against X4 tropism and its higher virulence. Lastly, we find a subtype D specific codon deletion at position 24 of the V3 loop, which may explain the higher propensity for subtype D to utilise X4 tropism.

Profiling genome-wide recombination in Epstein Barr virus reveals type-specific patterns and associations with endemic-Burkitt lymphoma

BACKGROUND

Endemic Burkitt lymphoma (eBL) is potentiated through the interplay of Epstein Barr virus (EBV) and holoendemic Plasmodium falciparum malaria. To better understand EBV's biology and role in eBL, we characterized genome-wide recombination sites and patterns as a source of genetic diversity in EBV genomes in our well-defined population of eBL cases and controls from Western Kenya.

METHODS

EBV genomes representing 54 eBL cases and 32 healthy children from the same geographic region in Western Kenya that we previously sequenced were analyzed. Whole-genome multiple sequence alignment, recombination analyses, and phylogenetic inference were made using multiple alignment with fast Fourier transform, recombination detection program 4, and molecular evolutionary genetics analysis.

RESULTS

We identified 28 different recombination events and 71 (82.6%) of the 86 EBV genomes analyzed contained evidence of one or more recombinant segments. Associated recombination breakpoints were found to occur in a total of 42 different genes, with only 7 (16.67%) being latent genes. Recombination events were major drivers of clustering within genome-wide phylogenetic trees. The occurrence of recombination segments was comparable between genomes from male and female participants and across age groups. More recombinant segments were found in EBV type 1 genomes (p = 6.4e - 06) and the genomes from the eBLs (p = 0.037). Two recombination events were enriched in the eBLs; event 47 (OR = 4.07, p = 0.038) and event 50 (OR = 14.24, p = 0.012).

CONCLUSIONS

EBV genomes have extensive evidence of recombination likely acquired progressively and cumulatively over time. Recombination patterns display a heterogeneous occurrence rate across the genome with enrichment in lytic genes. Overall, recombination appears to be a major evolutionary force impacting EBV diversity and genome structure with evidence of the association of specific recombinants with eBL.

ProbeTools: designing hybridization probes for targeted genomic sequencing of diverse and hypervariable viral taxa

BACKGROUND

Sequencing viruses in many specimens is hindered by excessive background material from hosts, microbiota, and environmental organisms. Consequently, enrichment of target genomic material is necessary for practical high-throughput viral genome sequencing. Hybridization probes are widely used for enrichment in many fields, but their application to viral sequencing faces a major obstacle: it is difficult to design panels of probe oligo sequences that broadly target many viral taxa due to their rapid evolution, extensive diversity, and genetic hypervariability. To address this challenge, we created ProbeTools, a package of bioinformatic tools for generating effective viral capture panels, and for assessing coverage of target sequences by probe panel designs in silico. In this study, we validated ProbeTools by designing a panel of 3600 probes for subtyping the hypervariable haemagglutinin (HA) and neuraminidase (NA) genome segments of avian-origin influenza A viruses (AIVs). Using in silico assessment of AIV reference sequences and in vitro capture on egg-cultured viral isolates, we demonstrated effective performance by our custom AIV panel and ProbeTools' suitability for challenging viral probe design applications.

RESULTS

Based on ProbeTool's in silico analysis, our panel provided broadly inclusive coverage of 14,772 HA and 11,967 NA reference sequences. For each reference sequence, we calculated the percentage of nucleotide positions covered by our panel in silico; 90% of HA and NA references sequences had at least 90.8 and 95.1% of their nucleotide positions covered respectively. We also observed effective in vitro capture on a representative collection of 23 egg-cultured AIVs that included isolates from wild birds, poultry, and humans and representatives from all HA and NA subtypes. Forty-two of forty-six HA and NA segments had over 98.3% of their nucleotide positions significantly enriched by our custom panel. These in vitro results were further used to validate ProbeTools' in silico coverage assessment algorithm; 89.2% of in silico predictions were concordant with in vitro results.

CONCLUSIONS

ProbeTools generated an effective panel for subtyping AIVs that can be deployed for genomic surveillance, outbreak prevention, and pandemic preparedness. Effective probe design against hypervariable AIV targets also validated ProbeTools' design and coverage assessment algorithms, demonstrating their suitability for other challenging viral capture applications.

High genetic stability of co-circulating human adenovirus type 31 lineages over 59 years

Type 31 of human adenovirus species A (HAdV-A31) is a significant pathogen primarily associated with diarrhoea in children but also with life-threatening disseminated disease in allogeneic haematopoietic stem cell transplant (HSCT) recipients. Nosocomial outbreaks of HAdV-A31 have been frequently described. However, the evolution of HAdV-A31 has not been studied in detail. The evolution of other HAdV types is driven either by intertypic recombination, where different types exchange genome regions, or by immune escape selection of neutralisation determinants. Complete genomic HAdV-A31 sequences from sixty diagnostic specimens of the past 18 years (2003-21) were generated, including fourteen specimens of a presumed outbreak on two HSCT wards. Additionally, twenty-three complete genomes from GenBank were added to our phylogenetic analysis as well as in silico generated and previously published restriction fragment polymorphism (RFLP) data. Phylogenetic analysis of eighty-three genomes indicated that HAdV-A31 evolved slowly with six lineages co-circulating. The two major lineages were lineage 1, which included the prototype from 1962 and nine recent isolates, and lineage 2, which split into four sublineages and included most isolates from 2003 to 2021. The average nucleotide identity within lineages was high (99.8 per cent) and identity between lineages was 98.7 and 99.2 per cent. RFLP data allowed the construction of a lower-resolution phylogeny with two additional putative lineages. Surprisingly, regions of higher diversity separating lineages were found in gene regions coding for non-structural and minor capsid proteins. Intertypic recombinations were not observed, but the phylogeny of lineage 3 was compatible with an interlineage recombination event in the fibre gene. Applying the phylogenetic analysis to the presumed nosocomial outbreak excluded two suspected transmission events and separated it into two different, simultaneous outbreaks caused by different sublineages of lineage 2. However, due to the high nucleotide identity within HAdV-A31 lineages, the proof of infection chains remains debatable. This in-depth study on the molecular phylogeny of HAdV-A31 highlights the high genetic stability of co-circulating HAdV-A31 lineages over almost six decades. It also supports the epidemiological hypothesis that HAdV-A31 circulates as an etiological agent of a childhood disease infecting immunologically naive patients without strong positive selection of immune escape variants and recombinants.

EBV genome variations enhance clinicopathological features of nasopharyngeal carcinoma in a non-endemic region

Nasopharyngeal carcinoma (NPC) is caused by infection with Epstein–Barr virus (EBV) and endemic in certain geographic regions. EBV lytic gene, BALF2, closely associates with viral reactivation and BALF2 gene variation, the H‐H‐H strain, causes NPC in endemic region, southern China. Here, we investigate whether such EBV variations also affect NPC in a non‐endemic region, Japan. Viral genome sequencing with 47 EBV isolates of Japanese NPC were performed and compared with those of other EBV‐associated diseases from Japan or NPC in Southern China. EBV genomes of Japanese NPC are different from those of other diseases in Japan or endemic NPC; Japanese NPC was not affected by the endemic strain (the BALF2 H‐H‐H) but frequently carried the type 2 EBV or the strain with intermediate risk of endemic NPC (the BALF2 H‐H‐L). Seven single nucleotide variations were specifically associated with Japanese NPC, of which six were present in both type 1 and 2 EBV genomes, suggesting the contribution of the type 2 EBV‐derived haplotype. This observation was supported by a higher viral titer and stronger viral reactivation in NPC with either type 2 or H‐H‐L strains. Our results highlight the importance of viral strains and viral reactivation in the pathogenesis of non‐endemic NPC.

Do Epstein–Barr Virus Mutations and Natural Genome Sequence Variations Contribute to Disease?

Most of the world’s population is infected by the Epstein–Barr virus (EBV), but the incidence of the diseases associated with EBV infection differs greatly in different parts of the world. Many factors may determine those differences, but variation in the virus genome is likely to be a contributing factor for some of the diseases. Here, we describe the main forms of EBV genome sequence variation, and the mechanisms by which variations in the virus genome are likely to contribute to disease. EBV genome deletions or polymorphisms can also provide useful markers for monitoring disease. If some EBV strains prove to be more pathogenic than others, this suggests the possible value of immunising people against infection by those pathogenic strains.

Deep sequencing of HIV-1 reveals extensive subtype variation and drug resistance after failure of first-line antiretroviral regimens in Nigeria

BACKGROUND

Deep sequencing could improve understanding of HIV treatment failure and viral population dynamics. However, this tool is often inaccessible in low- and middle-income countries.

OBJECTIVES

To determine the genetic patterns of resistance emerging in West African HIV-1 subtypes during first-line virological failure, and the implications for future antiretroviral options.

PATIENTS AND METHODS

Participants were selected from a Nigerian cohort of people living with HIV who had failed first-line ART and subsequently switched to second-line therapy. Whole HIV-1 genome sequences were generated from first-line virological failure samples with Illumina MiSeq. Mutations detected at ≥2% frequency were analysed and compared by subtype.

RESULTS

HIV-1 sequences were obtained from 101 participants (65% female, median age 30 years, median 32.9 months of nevirapine- or efavirenz-based ART). Thymidine analogue mutations (TAMs) were detected in 61%, other core NRTI mutations in 92% and NNRTI mutations in 99%. Minority variants (<20% frequency) comprised 18% of all mutations. K65R was more prevalent in CRF02_AG than G subtypes (33% versus 7%; P = 0.002), and ≥3 TAMs were more common in G than CRF02_AG (52% versus 24%; P = 0.004). Subtype G viruses also contained more RT cleavage site mutations. Cross-resistance to at least one of the newer NNRTIs, doravirine, etravirine or rilpivirine, was predicted in 81% of participants.

CONCLUSIONS

Extensive drug resistance had accumulated in people with West African HIV-1 subtypes, prior to second-line ART. Deep sequencing significantly increased the detection of resistance-associated mutations. Caution should be used if considering newer-generation NNRTI agents in this setting.

In action-an early warning system for the detection of unexpected or novel pathogens

Proactive approaches in preventing future epidemics include pathogen discovery prior to their emergence in human and/or animal populations. Playing an important role in pathogen discovery, high-throughput sequencing (HTS) enables the characterization of microbial and viral genetic diversity within a given sample. In particular, metagenomic HTS allows the unbiased taxonomic profiling of sequences; hence, it can identify novel and highly divergent pathogens such as viruses. Newly discovered viral sequences must be further investigated using genomic characterization, molecular and serological screening, and/or invitro and invivo characterization. Several outbreak and surveillance studies apply unbiased generic HTS to characterize the whole genome sequences of suspected pathogens. In contrast, this study aimed to screen for novel and unexpected pathogens in previously generated HTS datasets and use this information as a starting point for the establishment of an early warning system (EWS). As a proof of concept, the EWS was applied to HTS datasets and archived samples from the 2018–9 West Nile virus (WNV) epidemic in Germany. A metagenomics read classifier detected sequences related to genome sequences of various members of Riboviria. We focused the further EWS investigation on viruses belonging to the families Peribunyaviridae and Reoviridae, under suspicion of causing co-infections in WNV-infected birds. Phylogenetic analyses revealed that the reovirus genome sequences clustered with sequences assigned to the species Umatilla virus (UMAV), whereas a new peribunyavirid, tentatively named ‘Hedwig virus’ (HEDV), belonged to a putative novel genus of the family Peribunyaviridae. In follow-up studies, newly developed molecular diagnostic assays detected fourteen UMAV-positive wild birds from different German cities and eight HEDV-positive captive birds from two zoological gardens. UMAV was successfully cultivated in mosquito C6/36 cells inoculated with a blackbird liver. In conclusion, this study demonstrates the power of the applied EWS for the discovery and characterization of unexpected viruses in repurposed sequence datasets, followed by virus screening and cultivation using archived sample material. The EWS enhances the strategies for pathogen recognition before causing sporadic cases and massive outbreaks and proves to be a reliable tool for modern outbreak preparedness.

NK-/T-cell lymphomas

Natural killer/T-cell lymphoma (NKTL) is a sub-type of Epstein-Barr virus (EBV)-related non-Hodgkin lymphomas common in Asia and Latin America but rare elsewhere. Its pathogenesis is complex and incompletely understood. Lymphoma cells are transformed from NK- or T-cells, sometimes both. EBV-infection and subsequent genetic alterations in infected cells are central to NKTL development. Hemophagocytic syndrome is a common complication. Accurate staging is important to predict outcomes but there is controversy which system is best. More than two-thirds of NKTL lympohmas are localized at diagnosis, are frequently treated with radiation therapy only and have 5-year survival of about 70 percent. Persons with advanced NKTLs receive radiation therapy synchronously or metachronously with diverse multi-drug chemotherapy typically including L-asparginase with 5-year survival of about 40 percent. Some persons with widespread NKTL receive chemotherapy only. There are few data on safety and efficacy of high-dose therapy and a haematopoietic cell autotransplant. Immune therapies, histone deacetylase (HDAC)-inhibitors and other drugs are in early clinical trials. There are few randomized controlled clinical trials in NKTLs and no therapy strategy is clearly best; more effective therapy(ies) are needed. Some consensus recommendations are not convincingly evidence-based. Mechanisms of multi-drug resistance are considered. We discuss these issues including recent advances in our understanding of and therapy of NKTLs.

Varicella-Zoster Virus-Genetics, Molecular Evolution and Recombination

This chapter first details the structure, organization and coding content of the VZV genome to provide a foundation on which the molecular evolution of the virus can be projected. We subsequently describe the evolution of molecular profiling approaches from restriction fragment length polymorphisms to single nucleotide polymorphism profiling to modern day high-throughput sequencing approaches. We describe how the application of these methodologies led to our current model of VZV phylogeograpy including the number and structure of geographic clades and the role of recombination in reshaping these.

Human Papillomavirus Detection by Whole-Genome Next-Generation Sequencing: Importance of Validation and Quality Assurance Procedures

Next-generation sequencing (NGS) yields powerful opportunities for studying human papillomavirus (HPV) genomics for applications in epidemiology, public health, and clinical diagnostics. HPV genotypes, variants, and point mutations can be investigated in clinical materials and described in previously unprecedented detail. However, both the NGS laboratory analysis and bioinformatical approach require numerous steps and checks to ensure robust interpretation of results. Here, we provide a step-by-step review of recommendations for validation and quality assurance procedures of each step in the typical NGS workflow, with a focus on whole-genome sequencing approaches. The use of directed pilots and protocols to ensure optimization of sequencing data yield, followed by curated bioinformatical procedures, is particularly emphasized. Finally, the storage and sharing of data sets are discussed. The development of international standards for quality assurance should be a goal for the HPV NGS community, similar to what has been developed for other areas of sequencing efforts including microbiology and molecular pathology. We thus propose that it is time for NGS to be included in the global efforts on quality assurance and improvement of HPV-based testing and diagnostics.

Using Whole Genome Sequences to Investigate Adenovirus Outbreaks in a Hematopoietic Stem Cell Transplant Unit

A recent surge in human mastadenovirus (HAdV) cases, including five deaths, amongst a haematopoietic stem cell transplant population led us to use whole genome sequencing (WGS) to investigate. We compared sequences from 37 patients collected over a 20-month period with sequences from GenBank and our own database of HAdVs. Maximum likelihood trees and pairwise differences were used to evaluate genotypic relationships, paired with the epidemiological data from routine infection prevention and control (IPC) records and hospital activity data. During this time period, two formal outbreaks had been declared by IPC, while WGS detected nine monophyletic clusters, seven were corroborated by epidemiological evidence and by comparison of single-nucleotide polymorphisms. One of the formal outbreaks was confirmed, and the other was not. Of the five HAdV-associated deaths, three were unlinked and the remaining two considered the source of transmission. Mixed infection was frequent (10%), providing a sentinel source of recombination and superinfection. Immunosuppressed patients harboring a high rate of HAdV positivity require comprehensive surveillance. As a consequence of these findings, HAdV WGS is being incorporated routinely into clinical practice to influence IPC policy contemporaneously.

Viral metagenomics reveals the presence of novel Zika virus variants in Aedes mosquitoes from Barbados

Background

The Zika virus (ZIKV) epidemic of 2015/2016 spread throughout numerous countries. It emerged in mainland Latin America and spread to neighboring islands, including the Caribbean island of Barbados. Recent studies have indicated that the virus must have already been circulating in local mosquito populations in Brazil for almost 2 years before it was identified by the World Health Organization in 2015. Metagenomic detection assays have the potential to detect emerging pathogens without prior knowledge of their genomic nucleic acid sequence. Yet their applicability as vector surveillance tools has been widely limited by the complexity of DNA populations from field-collected mosquito preparations. The aim of this study was to investigate local vector biology and characterize metagenomic arbovirus diversity in Aedes mosquitoes during the ongoing 2015/2016 ZIKV epidemic.

Methods

We performed a short-term vector screening study on the island of Barbados during the ongoing 2015/2016 ZIKV epidemic, where we sampled local Aedes mosquitoes. We reanalyzed mosquito viral microbiome data derived from standard Illumina MiSeq sequencing to detect arbovirus sequences. Additionally, we employed deep sequencing techniques (Illumina HiSeq) and designed a novel bait capture enrichment assay to increase sequencing efficiency for arbovirus sequences from complex DNA samples.

Results

We found that Aedes aegypti seemed to be the most likely vector of ZIKV, although it prevailed at a low density during the observed time period. The number of detected viruses increased with sequencing depth. Arbovirus sequence enrichment of metagenomic DNA preparations allowed the detection of arbovirus sequences of two different ZIKV genotypes, including a novel one. To our knowledge, this is the first report of the S3116W mutation in the NS5 gene region of ZIKV polyprotein.

Conclusions

The metagenomic arbovirus detection approach presented here may serve as a useful tool for the identification of epidemic-causing arboviruses with the additional benefit of enabling the collection of phylogenetic information on the source. Apart from detecting more than 88 viruses using this approach, we also found evidence of novel ZIKV variants circulating in the local mosquito population during the observed time period.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04840-0.

Comprehensive Evolutionary Analysis of Complete Epstein-Barr Virus Genomes from Argentina and Other Geographies

The sequence variability of the Epstein–Barr virus has been extensively studied throughout previous years in isolates from various geographic regions and consequent variations at both genetic and genomic levels have been described. However, isolates from South America were underrepresented in these studies. Here, we sequenced 15 complete EBV genomes that we analyzed together with publicly available raw NGS data for 199 EBV isolates from other parts of the globe by means of a custom-built bioinformatic pipeline. The phylogenetic relations of the genomes, the geographic structure and variability of the data set, and the evolution rates for the whole genome and each gene were assessed. The present work contributes to overcoming the scarcity of complete EBV genomes from South America and is the most comprehensive geography-related variability study, which involved determining the actual contribution of each EBV gene to the geographic segregation of the entire genome. Moreover, to the best of our knowledge, we established for the first time the evolution rate for the entire EBV genome based on a host–virus codivergence-independent assumption and assessed their evolution rates on a gene-by-gene basis, which were related to the encoded protein function. Considering the evolution of dsDNA viruses with a codivergence-independent approach may lay the basis for future research on EBV evolution. The exhaustive bioinformatic analysis performed on this new dataset allowed us to draw a novel set of conclusions regarding the genome evolution of EBV.

Evolutionary History of Endogenous Human Herpesvirus 6 Reflects Human Migration out of Africa

Human herpesvirus 6A and 6B (HHV-6) can integrate into the germline, and as a result, ∼70 million people harbor the genome of one of these viruses in every cell of their body. Until now, it has been largely unknown if 1) these integrations are ancient, 2) if they still occur, and 3) whether circulating virus strains differ from integrated ones. Here, we used next-generation sequencing and mining of public human genome data sets to generate the largest and most diverse collection of circulating and integrated HHV-6 genomes studied to date. In genomes of geographically dispersed, only distantly related people, we identified clades of integrated viruses that originated from a single ancestral event, confirming this with fluorescent in situ hybridization to directly observe the integration locus. In contrast to HHV-6B, circulating and integrated HHV-6A sequences form distinct clades, arguing against ongoing integration of circulating HHV-6A or “reactivation” of integrated HHV-6A. Taken together, our study provides the first comprehensive picture of the evolution of HHV-6, and reveals that integration of heritable HHV-6 has occurred since the time of, if not before, human migrations out of Africa.

Human Cytomegalovirus Genome Diversity in Longitudinally Collected Breast Milk Samples

Reactivation and shedding of human cytomegalovirus (HCMV) in breast milk during lactation is highly frequent in HCMV-seropositive mothers. This represents a key transmission route for postnatal HCMV infection and can lead to severe disease in preterm neonates. Little is known about HCMV strain composition or longitudinal intrahost viral population dynamics in breast milk from immunocompetent women. We performed HCMV-specific target enrichment and high-throughput sequencing of 38 breast milk samples obtained in Germany between days 10 and 60 postpartum from 15 mothers with HCMV DNA lactia, and assembled HCMV consensus sequences de novo. The genotype distribution and number of HCMV strains present in each sample were determined by quantifying genotype-specific sequence motifs in 12 hypervariable viral genes, revealing a wide range of genotypes (82/109) for these genes in the cohort and a unique, longitudinally stable strain composition in each mother. Reactivation of up to three distinct HCMV strains was detected in 8/15 of mothers, indicating that a representative subset of the woman’s HCMV reservoir might be locally reactivated early during lactation. As described previously, nucleotide diversity of samples with multiple strains was much higher than that of samples with single strains. Breast milk as a main source of postnatal mother-to-infant transmission may serve as a repository for viral diversity and thus play an essential role in the natural epidemiology of HCMV.

Identification of UL69 Gene and Protein in Cytomegalovirus-Transformed Human Mammary Epithelial Cells

A growing body of evidence addressing the involvement of human cytomegalovirus (HCMV) in malignancies had directed attention to the oncomodulation paradigm. HCMV-DB infected human mammary epithelial cells (HMECs) in culture showed the emergence of clusters of rapidly proliferating, spheroid-shaped transformed cells named CTH (CMV-Transformed HMECs) cells. CTH cells assessment suggests a direct contribution of HCMV to oncogenesis, from key latent and lytic genes activating oncogenic pathways to fueling tumor evolution. We hypothesized that the presence of HCMV genome in CTH cells is of pivotal importance for determining its oncogenic potential. We previously reported the detection of a long non-coding (lnc) RNA4.9 gene in CTH cells. Therefore, we assessed here the presence of UL69 gene, located nearby and downstream of the lncRNA4.9 gene, in CTH cells. The HCMV UL69 gene in CTH cells was detected using polymerase chain reaction (PCR) and sequencing of UL69 gene was performed using Sanger method. The corresponding amino acid sequence was then blasted against the UL69 sequence derived from HCMV-DB genome using NCBI Protein BLAST tool. A 99% identity was present between the nucleotide sequence present in CTH cells and HCMV-DB genome. UL69 transcript was detected in RNA extracts of CTH cells, using a reverse transcription polymerase chain reaction (RT-PCR) assay, and pUL69 protein was identified in CTH lysates using western blotting. Ganciclovir-treated CTH cells showed a decrease in UL69 gene detection and cellular proliferation. In CTH cells, the knockdown of UL69 with siRNA was assessed by RT-qPCR and western blot to reveal the impact of pUL69 on HCMV replication and CTH cell proliferation. Finally, UL69 gene was detected in breast cancer biopsies. Our results indicate a close link between the UL69 gene detected in the HCMV-DB isolate used to infect HMECs, and the UL69 gene present in transformed CTH cells and tumor biopsies, further highlighting a direct role for HCMV in breast tumor development.

The influence of human genetic variation on Epstein-Barr virus sequence diversity

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

Next-generation diagnostics: virus capture facilitates a sensitive viral diagnosis for epizootic and zoonotic pathogens including SARS-CoV-2

BACKGROUND

The detection of pathogens in clinical and environmental samples using high-throughput sequencing (HTS) is often hampered by large amounts of background information, which is especially true for viruses with small genomes. Enormous sequencing depth can be necessary to compile sufficient information for identification of a certain pathogen. Generic HTS combining with in-solution capture enrichment can markedly increase the sensitivity for virus detection in complex diagnostic samples.

METHODS

A virus panel based on the principle of biotinylated RNA baits was developed for specific capture enrichment of epizootic and zoonotic viruses (VirBaits). The VirBaits set was supplemented by a SARS-CoV-2 predesigned bait set for testing recent SARS-CoV-2-positive samples. Libraries generated from complex samples were sequenced via generic HTS (without enrichment) and afterwards enriched with the VirBaits set. For validation, an internal proficiency test for emerging epizootic and zoonotic viruses (African swine fever virus, Ebolavirus, Marburgvirus, Nipah henipavirus, Rift Valley fever virus) was conducted.

RESULTS

The VirBaits set consists of 177,471 RNA baits (80-mer) based on about 18,800 complete viral genomes targeting 35 epizootic and zoonotic viruses. In all tested samples, viruses with both DNA and RNA genomes were clearly enriched ranging from about 10-fold to 10,000-fold for viruses including distantly related viruses with at least 72% overall identity to viruses represented in the bait set. Viruses showing a lower overall identity (38% and 46%) to them were not enriched but could nonetheless be detected based on capturing conserved genome regions. The internal proficiency test supports the improved virus detection using the combination of HTS plus targeted enrichment but also points to the risk of cross-contamination between samples.

CONCLUSIONS

The VirBaits approach showed a high diagnostic performance, also for distantly related viruses. The bait set is modular and expandable according to the favored diagnostics, health sector, or research question. The risk of cross-contamination needs to be taken into consideration. The application of the RNA-baits principle turned out to be user friendly, and even non-experts can easily use the VirBaits workflow. The rapid extension of the established VirBaits set adapted to actual outbreak events is possible as shown for SARS-CoV-2. Video abstract.

Next-generation diagnostics: virus capture facilitates a sensitive viral diagnosis for epizootic and zoonotic pathogens including SARS-CoV-2

BACKGROUND

The detection of pathogens in clinical and environmental samples using high-throughput sequencing (HTS) is often hampered by large amounts of background information, which is especially true for viruses with small genomes. Enormous sequencing depth can be necessary to compile sufficient information for identification of a certain pathogen. Generic HTS combining with in-solution capture enrichment can markedly increase the sensitivity for virus detection in complex diagnostic samples.

METHODS

A virus panel based on the principle of biotinylated RNA baits was developed for specific capture enrichment of epizootic and zoonotic viruses (VirBaits). The VirBaits set was supplemented by a SARS-CoV-2 predesigned bait set for testing recent SARS-CoV-2-positive samples. Libraries generated from complex samples were sequenced via generic HTS (without enrichment) and afterwards enriched with the VirBaits set. For validation, an internal proficiency test for emerging epizootic and zoonotic viruses (African swine fever virus, Ebolavirus, Marburgvirus, Nipah henipavirus, Rift Valley fever virus) was conducted.

RESULTS

The VirBaits set consists of 177,471 RNA baits (80-mer) based on about 18,800 complete viral genomes targeting 35 epizootic and zoonotic viruses. In all tested samples, viruses with both DNA and RNA genomes were clearly enriched ranging from about 10-fold to 10,000-fold for viruses including distantly related viruses with at least 72% overall identity to viruses represented in the bait set. Viruses showing a lower overall identity (38% and 46%) to them were not enriched but could nonetheless be detected based on capturing conserved genome regions. The internal proficiency test supports the improved virus detection using the combination of HTS plus targeted enrichment but also points to the risk of cross-contamination between samples.

CONCLUSIONS

The VirBaits approach showed a high diagnostic performance, also for distantly related viruses. The bait set is modular and expandable according to the favored diagnostics, health sector, or research question. The risk of cross-contamination needs to be taken into consideration. The application of the RNA-baits principle turned out to be user friendly, and even non-experts can easily use the VirBaits workflow. The rapid extension of the established VirBaits set adapted to actual outbreak events is possible as shown for SARS-CoV-2. Video abstract.

Intra-host changes in Kaposi sarcoma-associated herpesvirus genomes in Ugandan adults with Kaposi sarcoma

Intra-host tumor virus variants may influence the pathogenesis and treatment responses of some virally-associated cancers. However, the intra-host variability of Kaposi sarcoma-associated herpesvirus (KSHV), the etiologic agent of Kaposi sarcoma (KS), has to date been explored with sequencing technologies that possibly introduce more errors than that which occurs in the viral population, and these studies have only studied variable regions. Here, full-length KSHV genomes in tumors and/or oral swabs from 9 Ugandan adults with HIV-associated KS were characterized. Furthermore, we used deep, short-read sequencing using duplex unique molecular identifiers (dUMI)–random double-stranded oligonucleotides that barcode individual DNA molecules before library amplification. This allowed suppression of PCR and sequencing errors to ~10⁻⁹/base as well as afforded accurate determination of KSHV genome numbers sequenced in each sample. KSHV genomes were assembled de novo, and rearrangements observed were confirmed by PCR and Sanger sequencing. 131-kb KSHV genome sequences, excluding major repeat regions, were successfully obtained from 23 clinical specimens, averaging 2.3x10⁴ reads/base. Strikingly, KSHV genomes were virtually identical within individuals at the point mutational level. The intra-host heterogeneity that was observed was confined to tumor-associated KSHV mutations and genome rearrangements, all impacting protein-coding sequences. Although it is unclear whether these changes were important to tumorigenesis or occurred as a result of genomic instability in tumors, similar changes were observed across individuals. These included inactivation of the K8.1 gene in tumors of 3 individuals and retention of a region around the first major internal repeat (IR1) in all instances of genomic deletions and rearrangements. Notably, the same breakpoint junctions were found in distinct tumors within single individuals, suggesting metastatic spread of rearranged KSHV genomes. These findings define KSHV intra-host heterogeneity in vivo with greater precision than has been possible in the past and suggest the possibility that aberrant KSHV genomes may contribute to aspects of KS tumorigenesis. Furthermore, study of KSHV with use of dUMI provides a proof of concept for utilizing this technique for detailed study of other virus populations in vivo.

Kaposi sarcoma (KS) is a leading cancer in sub-Saharan Africa and in persons with HIV co-infection. Kaposi sarcoma-associated herpesvirus (KSHV, also referred to as human herpesvirus-8, or HHV-8) is the etiologic agent of KS, but the factors that contribute to the development of KS, which occurs in only a small subset of infected individuals, remain largely unknown. While strain differences or mutations in other tumor viruses are known to affect the risk and progression of their associated cancers, whether genetic variation in KSHV is important to the natural history of KS is unclear. Most studies of KSHV diversity have only characterized ~4% of its 165-kb genome, and the observed variation in some studies is likely to have been impacted by PCR or cloning artifacts. To precisely define genomic diversity of KSHV in vivo, we evaluated full-length viral genomes (except the internal repeat regions) using a technique that greatly lowers sequencing error rates and thus measures genomic diversity much more accurately than previous studies. In addition, we extended our analyses to look for potential tumor-specific changes in the KSHV genomes by examining viruses in both tumor and non-tumor tissues. To these ends, we performed highly sensitive, single-molecule sequencing of whole KSHV genomes in paired KS tumors and oral swabs from 9 individuals with KS. We found that KSHV genomes were virtually identical within the 9 individuals, with no evidence of quasispecies formation or multi-strain infection. However, KSHV genome aberrations and gene-inactivating mutations were found to be common in KS tumors, often impacting the same genes and genomic regions across individuals. Whether theses mutations influence KS tumorigenesis or result from genomic instability commonly found in tumors warrants further study. Lastly, aberrant KSHV genomes were found to be shared by distinct tumors within individuals, suggesting the capacity of KS tumor cells to metastasize and seed new lesions.

A global phylogenetic analysis of Japanese tonsil-derived Epstein-Barr virus strains using viral whole-genome cloning and long-read sequencing

Epstein-Barr virus (EBV) establishes lifelong latent infection in the majority of healthy individuals, while it is a causative agent for various diseases, including some malignancies. Recent high-throughput sequencing results indicate that there are substantial levels of viral genome heterogeneity among different EBV strains. However, the extent of EBV strain variation among asymptomatically infected individuals remains elusive. Here, we present a streamlined experimental strategy to clone and sequence EBV genomes derived from human tonsillar tissues, which are the reservoirs of asymptomatic EBV infection. Complete EBV genome sequences, including those of repetitive regions, were determined for seven tonsil-derived EBV strains. Phylogenetic analyses based on the whole viral genome sequences of worldwide non-tumour-derived EBV strains revealed that Asian EBV strains could be divided into several distinct subgroups. EBV strains derived from nasopharyngeal carcinoma-endemic areas constitute different subgroups from a subgroup of EBV strains from non-endemic areas, including Japan. The results could be consistent with biased regional distribution of EBV-associated diseases depending on the different EBV strains colonizing different regions in Asian countries.

Whole-Genome Sequencing of Mycobacterium tuberculosis Directly from Sputum Samples

Whole-genome sequencing is a powerful, high-resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography, and mutations associated with antimicrobial resistance. The ability to sequence pathogen genomes directly from clinical specimens, without the requirement for in vitro culturing, is attractive in terms of time- and labor-saving, especially in the case of slow growing pathogens, such as Mycobacterium tuberculosis. However, clinical samples typically contain too low levels of pathogen nucleic acid, plus relatively high levels of human and natural microbiota DNA/RNA, to make this a viable option. Using a combination of whole-genome enrichment and deep sequencing, which has been proven to be a nonmutagenic approach, we can capture all known variations found within M. tuberculosis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of M. tuberculosis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

Human cytomegalovirus multiple-strain infections and viral population diversity in haematopoietic stem cell transplant recipients analysed by high-throughput sequencing

Human cytomegalovirus (HCMV) is an important opportunistic pathogen in allogeneic haematopoietic stem cell transplant (HSCT) recipients. High-throughput sequencing of target-enriched libraries was performed to characterise the diversity of HCMV strains present in this high-risk group. Forty-four HCMV-DNA-positive plasma specimens (median viral input load 321 IU per library) collected at defined time points from 23 HSCT recipients within 80 days of transplantation were sequenced. The genotype distribution for 12 hypervariable HCMV genes and the number of HCMV strains present (i.e. single- vs. multiple-strain infection) were determined for 29 samples from 16 recipients. Multiple-strain infection was observed in seven of these 16 recipients, and five of these seven recipients had the donor (D)/recipient (R) HCMV-serostatus combination D + R + . A very broad range of genotypes was detected, with an intrahost composition that was generally stable over time. Multiple-strain infection was not associated with particular virological or clinical features, such as altered levels or duration of antigenaemia, development of acute graft-versus-host disease or increased mortality. In conclusion, despite relatively low viral plasma loads, a high frequency of multiple-strain HCMV infection and a high strain complexity were demonstrated in systematically collected clinical samples from this cohort early after HSCT. However, robust evaluation of the pathogenic role of intrahost viral diversity and multiple-strain infection will require studies enrolling larger numbers of recipients.

Whole genome sequencing of Theileria parva using target capture

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0-33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.

Evaluation of targeted next-generation sequencing for detection of equine pathogens in clinical samples

Equine infectious disease outbreaks may have profound economic impact, resulting in losses of millions of dollars of revenue as a result of horse loss, quarantine, and cancelled events. Early and accurate diagnosis is essential to limit the spread of infectious diseases. However, laboratory detection of infectious agents, especially the simultaneous detection of multiple agents, can be challenging to the clinician and diagnostic laboratory. Next-generation sequencing (NGS), which allows millions of DNA templates to be sequenced simultaneously in a single reaction, is an ideal technology for comprehensive testing. We conducted a proof-of-concept study of targeted NGS to detect 62 common equine bacterial, viral, and parasitic pathogens in clinical samples. We designed 264 primers and constructed a bioinformatics tool for the detection of targeted pathogens. The designed primers were able to specifically detect the intended pathogens. Results of testing 27 clinical samples with our targeted NGS assay compared with results of routine tests (assessed as a group) yielded positive percent agreement of 81% and negative percent agreement of 83%, overall agreement of 81%, and kappa of 0.56 (moderate agreement). This moderate agreement was likely the result of low sensitivity of some primers. However, our NGS assay successfully detected multiple pathogens in the clinical samples, including some pathogens missed by routine techniques.

Alphaherpesvirus Genomics: Past, Present and Future

Alphaherpesviruses, as large double-stranded DNA viruses, were long considered to be genetically stable and to exist in a homogeneous state. Recently, the proliferation of high-throughput sequencing (HTS) and bioinformatics analysis has expanded our understanding of herpesvirus genomes and the variations found therein. Recent data indicate that herpesviruses exist as diverse populations, both in culture and in vivo, in a manner reminiscent of RNA viruses. In this review, we discuss the past, present, and potential future of alphaherpesvirus genomics, including the technical challenges that face the field. We also review how recent data has enabled genome-wide comparisons of sequence diversity, recombination, allele frequency, and selective pressures, including those introduced by cell culture. While we focus on the human alphaherpesviruses, we draw key insights from related veterinary species and from the beta- and gamma-subfamilies of herpesviruses. Promising technologies and potential future directions for herpesvirus genomics are highlighted as well, including the potential to link viral genetic differences to phenotypic and disease outcomes.

Whole-genome Sequencing Provides Data for Stratifying Infection Prevention and Control Management of Nosocomial Influenza A

Background

Influenza A virus causes annual epidemics in humans and is associated with significant morbidity and mortality. Haemagglutinin (HA) and neuraminidase (NA) gene sequencing have traditionally been used to identify the virus genotype, although their utility in detecting outbreak clusters is still unclear. The objective of this study was to determine the utility, if any, of whole-genome sequencing over HA/NA sequencing for infection prevention and control (IPC) in hospitals.

Methods

We obtained all clinical samples from influenza (H1N1)-positive patients at the Great Ormond Street Hospital between January and March 2016. Samples were sequenced using targeted enrichment on an Illumina MiSeq sequencer. Maximum likelihood trees were computed for both whole genomes and concatenated HA/NA sequences. Epidemiological data was taken from routine IPC team activity during the period.

Results

Complete genomes were obtained for 65/80 samples from 38 patients. Conventional IPC analysis recognized 1 outbreak, involving 3 children, and identified another potential cluster in the haemato-oncology ward. Whole-genome and HA/NA phylogeny both accurately identified the previously known outbreak cluster. However, HA/NA sequencing additionally identified unrelated strains as part of this outbreak cluster. A whole-genome analysis identified a further cluster of 2 infections that had been previously missed and refuted suspicions of transmission in the haemato-oncology wards.

Conclusions

Whole-genome sequencing is better at identifying outbreak clusters in a hospital setting than HA/NA sequencing. Whole-genome sequencing could provide a faster and more reliable method for outbreak monitoring and supplement routine IPC team work to allow the prevention of transmission.

Epstein-Barr Virus Genomes Reveal Population Structure and Type 1 Association with Endemic Burkitt Lymphoma

Endemic Burkitt lymphoma (eBL), the most prevalent pediatric cancer in sub-Saharan Africa, is distinguished by its inclusion of Epstein-Barr virus (EBV). In order to better understand the impact of EBV variation in eBL tumorigenesis, we improved viral DNA enrichment methods and generated a total of 98 new EBV genomes from both eBL cases (n = 58) and healthy controls (n = 40) residing in the same geographic region in Kenya. Using our unbiased methods, we found that EBV type 1 was significantly more prevalent in eBL patients (74.5%) than in healthy children (47.5%) (odds ratio = 3.24, 95% confidence interval = 1.36 to 7.71, P = 0.007), as opposed to similar proportions in both groups. Controlling for EBV type, we also performed a genome-wide association study identifying six nonsynonymous variants in the genes EBNA1, EBNA2, BcLF1, and BARF1 that were enriched in eBL patients. In addition, viruses isolated from plasma of eBL patients were identical to their tumor counterparts consistent with circulating viral DNA originating from the tumor. We also detected three intertypic recombinants carrying type 1 EBNA2 and type 2 EBNA3 regions, as well as one novel genome with a 20-kb deletion, resulting in the loss of multiple lytic and virion genes. Comparing EBV types, viral genes displayed differential variation rates as type 1 appeared to be more divergent, while type 2 demonstrated novel substructures. Overall, our findings highlight the complexities of the EBV population structure and provide new insight into viral variation, potentially deepening our understanding of eBL oncogenesis.IMPORTANCE Improved viral enrichment methods conclusively demonstrate EBV type 1 to be more prevalent in eBL patients than in geographically matched healthy controls, which previously underrepresented the prevalence of EBV type 2. Genome-wide association analysis between cases and controls identifies six eBL-associated nonsynonymous variants in EBNA1, EBNA2, BcLF1, and BARF1 genes. Analysis of population structure reveals that EBV type 2 exists as two genomic subgroups and was more commonly found in female than in male eBL patients.

A Comparison of Whole Genome Sequencing of SARS-CoV-2 Using Amplicon-Based Sequencing, Random Hexamers, and Bait Capture

Genome sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is increasingly important to monitor the transmission and adaptive evolution of the virus. The accessibility of high-throughput methods and polymerase chain reaction (PCR) has facilitated a growing ecosystem of protocols. Two differing protocols are tiling multiplex PCR and bait capture enrichment. Each method has advantages and disadvantages but a direct comparison with different viral RNA concentrations has not been performed to assess the performance of these approaches. Here we compare Liverpool amplification, ARTIC amplification, and bait capture using clinical diagnostics samples. All libraries were sequenced using an Illumina MiniSeq with data analyzed using a standardized bioinformatics workflow (SARS-CoV-2 Illumina GeNome Assembly Line; SIGNAL). One sample showed poor SARS-CoV-2 genome coverage and consensus, reflective of low viral RNA concentration. In contrast, the second sample had a higher viral RNA concentration, which yielded good genome coverage and consensus. ARTIC amplification showed the highest depth of coverage results for both samples, suggesting this protocol is effective for low concentrations. Liverpool amplification provided a more even read coverage of the SARS-CoV-2 genome, but at a lower depth of coverage. Bait capture enrichment of SARS-CoV-2 cDNA provided results on par with amplification. While only two clinical samples were examined in this comparative analysis, both the Liverpool and ARTIC amplification methods showed differing efficacy for high and low concentration samples. In addition, amplification-free bait capture enriched sequencing of cDNA is a viable method for generating a SARS-CoV-2 genome sequence and for identification of amplification artifacts.

Advanced molecular surveillance approaches for characterization of blood borne hepatitis viruses

Defining genetic diversity of viral infections directly from patient specimens is the ultimate goal of surveillance. Simple tools that can provide full-length sequence information on blood borne viral hepatitis viruses: hepatitis C, hepatitis B and hepatitis D viruses (HCV, HBV and HDV) remain elusive. Here, an unbiased metagenomic next generation sequencing approach (mNGS) was used for molecular characterization of HCV infections (n = 99) from Israel which yielded full-length HCV sequences in 89% of samples, with 7 partial sequences sufficient for classification. HCV genotypes were primarily 1b (68%) and 1a (19%), with minor representation of genotypes 2c (1%) and 3a (8%). HBV/HDV coinfections were characterized by suppressed HBV viral loads, resulting in sparse mNGS coverage. A probe-based enrichment approach (xGen) aiming to increase HBV and HDV coverage was validated on a panel of diverse genotypes, geography and titers. The method extended HBV genome coverage a median 61% (range 8–84%) and provided orders of magnitude boosts in reads and sequence depth for both viruses. When HBV-xGen was applied to Israeli samples, coverage was improved by 28–73% in 4 samples and identified HBV genotype A1, A2, D1 specimens and a dual B/D infection. Abundant HDV reads in mNGS libraries yielded 18/26 (69%) full genomes and 8 partial sequences, with HDV-xGen only providing minimal extension (3–11%) of what were all genotype 1 genomes. Advanced molecular approaches coupled to virus-specific capture probes promise to enhance surveillance of viral infections and aid in monitoring the spread of local subtypes.

Assessing Genome-Wide Diversity in European Hantaviruses through Sequence Capture from Natural Host Samples

Research on the ecology and evolution of viruses is often hampered by the limitation of sequence information to short parts of the genomes or single genomes derived from cultures. In this study, we use hybrid sequence capture enrichment in combination with high-throughput sequencing to provide efficient access to full genomes of European hantaviruses from rodent samples obtained in the field. We applied this methodology to Tula (TULV) and Puumala (PUUV) orthohantaviruses for which analyses from natural host samples are typically restricted to partial sequences of their tri-segmented RNA genome. We assembled a total of ten novel hantavirus genomes de novo with very high coverage (on average >99%) and sequencing depth (average >247×). A comparison with partial Sanger sequences indicated an accuracy of >99.9% for the assemblies. An analysis of two common vole (Microtus arvalis) samples infected with two TULV strains each allowed for the de novo assembly of all four TULV genomes. Combining the novel sequences with all available TULV and PUUV genomes revealed very similar patterns of sequence diversity along the genomes, except for remarkably higher diversity in the non-coding region of the S-segment in PUUV. The genomic distribution of polymorphisms in the coding sequence was similar between the species, but differed between the segments with the highest sequence divergence of 0.274 for the M-segment, 0.265 for the S-segment, and 0.248 for the L-segment (overall 0.258). Phylogenetic analyses showed the clustering of genome sequences consistent with their geographic distribution within each species. Genome-wide data yielded extremely high node support values, despite the impact of strong mutational saturation that is expected for hantavirus sequences obtained over large spatial distances. We conclude that genome sequencing based on capture enrichment protocols provides an efficient means for ecological and evolutionary investigations of hantaviruses at an unprecedented completeness and depth.

Fetal inheritance of chromosomally integrated human herpesvirus 6 predisposes the mother to pre-eclampsia

Pre-eclampsia (typically characterized by new-onset hypertension and proteinuria in the second half of pregnancy) represents a major determinant of the global burden of disease^1,2. Its pathophysiology involves placental dysfunction, but the mechanism is unclear. Viral infection can cause organ dysfunction, but its role in placentally related disorders of human pregnancy is unknown³. We addressed this using RNA sequencing metagenomics^4-6 of placental samples from normal and complicated pregnancies. Here, we show that human herpesvirus 6 (HHV-6, A or B) RNA was detected in 6.1% of cases of pre-eclampsia and 2.2% of other pregnancies. Fetal genotyping demonstrated that 70% of samples with HHV-6 RNA in the placenta exhibited inherited, chromosomally integrated HHV-6 (iciHHV-6). We genotyped 467 pre-eclampsia cases and 3,854 controls and found an excess of iciHHV-6 in the cases (odds ratio of 2.8, 95% confidence intervals of 1.4-5.6, P = 0.008). We validated this finding by comparing iciHHV-6 in a further 740 cases with controls from large-scale population studies (odds ratio of 2.5, 95% confidence intervals of 1.4-4.4, P = 0.0013). We conclude that iciHHV-6 results in the transcription of viral RNA in the human placenta and predisposes the mother to pre-eclampsia.

Human Cytomegalovirus Genomes Sequenced Directly From Clinical Material: Variation, Multiple-Strain Infection, Recombination, and Gene Loss

The genomic characteristics of human cytomegalovirus (HCMV) strains sequenced directly from clinical pathology samples were investigated, focusing on variation, multiple-strain infection, recombination, and gene loss. A total of 207 datasets generated in this and previous studies using target enrichment and high-throughput sequencing were analyzed, in the process enabling the determination of genome sequences for 91 strains. Key findings were that (i) it is important to monitor the quality of sequencing libraries in investigating variation; (ii) many recombinant strains have been transmitted during HCMV evolution, and some have apparently survived for thousands of years without further recombination; (iii) mutants with nonfunctional genes (pseudogenes) have been circulating and recombining for long periods and can cause congenital infection and resulting clinical sequelae; and (iv) intrahost variation in single-strain infections is much less than that in multiple-strain infections. Future population-based studies are likely to continue illuminating the evolution, epidemiology, and pathogenesis of HCMV.

Whole-Genome Sequencing of Chlamydia trachomatis Directly from Human Samples

Whole-genome sequencing is a powerful, high-resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography and mutations associated with antimicrobial resistance. The ability to sequence pathogen genomes directly from clinical specimens, without the requirement for in vitro culturing, is attractive in terms of time- and labor-saving, especially in the case of slow growing, or obligate intracellular pathogens, such as Chlamydia trachomatis. However clinical samples typically contain too low levels of pathogen nucleic acid, plus relatively high levels of human and natural microbiota DNA/RNA, to make this a viable option. Using a combination of whole-genome enrichment and deep sequencing, which has been proven to be a nonmutagenic approach, we can capture all known variation found within C. trachomatis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of C. trachomatis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

Optimization of sample preparation for culture-independent sequencing of Bordetella pertussis

Bordetella pertussis, the aetiological agent of whooping cough, is re-emerging globally despite widespread vaccination. B. pertussis is highly infectious and, prior to vaccination programmes, was the leading cause of infant mortality. The WHO estimated that over 600 000 deaths are prevented annually by pertussis vaccination, but B. pertussis infection was still responsible for over 63 000 deaths globally in 2013. The re-emergence of B. pertussis has been linked to strains with inactive or absent major virulence factors included in vaccines such as pertactin, pertussis toxin and filamentous haemagglutinin. Thus, the molecular surveillance of currently circulating strains is critical in understanding and controlling B. pertussis. Such information provides data on strains to inform control measures and the identification of future vaccine antigens. Current surveillance and typing methods for B. pertussis rely on the availability of clinical isolates. However, since the 1990s, the majority of pertussis cases have been diagnosed by PCR, where an isolate is not needed. The rapid decline in the availability of B. pertussis isolates impacts our ability to monitor this infection. The growing uptake of next-generation sequencing (NGS) has offered the opportunity for culture-independent genome sequencing and typing of this fastidious pathogen. Therefore, the objective of the study was to optimize respiratory sample preparation, independent of culture, in order to type B. pertussis using NGS. The study compared commercial depletion kits and specimen-processing methods using selective lysis detergents. The goal was to deplete human DNA, the major obstacle for sequencing a pathogen directly from a clinical sample. Samples spiked with a clinically relevant amount of B. pertussis were used to provide comparison between the different methods. Commercial depletion kits including the MolYsis, Qiagen Microbiome and NEBNext Kits were tested. Previously published methods, for Saponin and TritonX-100, were also trialled as a depletion. The ratio of B. pertussis to human DNA was determined by real-time PCR for ERV3 and IS481 (as markers of human and B. pertussis DNA, respectively), then samples were sequenced using the Illumina NextSeq 500 platform. The number of human and B. pertussis sequenced reads were then compared between treatments. The results showed that commercial kits reduced the human DNA present, but also reduced the concentration of target B. pertussis. However, selective lysis with Saponin treatment resulted in almost undetectable levels of human DNA, with minimal loss of target B. pertussis DNA. Sequencing read depth improved five-fold in reads to B. pertussis. Our investigation delivered a potential protocol that will enable the public health laboratory surveillance of B. pertussis in the era of culture-independent testing.