Next Article in Journal
Rice Dwarf Virus Small RNA Profiles in Rice and Leafhopper Reveal Distinct Patterns in Cross-Kingdom Hosts
Previous Article in Journal
Follow-Up of Viral Parameters in FeLV- or FIV-Naturally Infected Cats Treated Orally with Low Doses of Human Interferon Alpha
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 11
Issue 9
10.3390/v11090846
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Deep-Sequencing Workflow for the Fast and Efficient Generation of High-Quality African Swine Fever Virus Whole-Genome Sequences

by
Jan H. Forth
1,
Leonie F. Forth
1,
Jacqueline King
1,
Oxana Groza
2,
Alexandra Hübner
1,
Ann Sofie Olesen
3,
Dirk Höper
1,
Linda K. Dixon
4,
Christopher L. Netherton
4,
Thomas Bruun Rasmussen
5,6,
Sandra Blome
1,*,
Anne Pohlmann
1 and
Martin Beer
1
1
Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Südufer 10, 17493 Greifswald - Insel Riems, Germany
2
Animal Health Diagnostic Laboratory, Republican Center for Veterinary Diagnostic, str. Murelor, 3, MD-2051, mun. Chișinău, Moldova
3
Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, DK-1870 Frederiksberg C, Denmark
4
The Pirbright Institute, Ash Road, Pirbright, Woking GU24 0NF, UK
5
DTU National Veterinary Institute, Technical University of Denmark, Lindholm, DK-4771 Kalvehave, Denmark
6
Department of Virus & Microbiological Special Diagnostics, Statens Serum Institut, DK-2300 Copenhagen S, Denmark
*
Author to whom correspondence should be addressed.
Viruses 2019, 11(9), 846; https://doi.org/10.3390/v11090846
Submission received: 5 August 2019 / Revised: 5 September 2019 / Accepted: 9 September 2019 / Published: 11 September 2019
(This article belongs to the Section Animal Viruses)
Browse Figures
Review Reports Versions Notes

Abstract

:
African swine fever (ASF) is a severe disease of suids caused by African swine fever virus (ASFV). Its dsDNA genome (170–194 kbp) is scattered with homopolymers and repeats as well as inverted-terminal-repeats (ITR), which hamper whole-genome sequencing. To date, only a few genome sequences have been published and only for some are data on sequence quality available enabling in-depth investigations. Especially in Europe and Asia, where ASFV has continuously spread since its introduction into Georgia in 2007, a very low genetic variability of the circulating ASFV-strains was reported. Therefore, only whole-genome sequences can serve as a basis for detailed virus comparisons. Here, we report an effective workflow, combining target enrichment, Illumina and Nanopore sequencing for ASFV whole-genome sequencing. Following this approach, we generated an improved high-quality ASFV Georgia 2007/1 whole-genome sequence leading to the correction of 71 sequencing errors and the addition of 956 and 231 bp at the respective ITRs. This genome, derived from the primary outbreak in 2007, can now serve as a reference for future whole-genome analyses of related ASFV strains and molecular approaches. Using both workflow and the reference genome, we generated the first ASFV-whole-genome sequence from Moldova, expanding the sequence knowledge from Eastern Europe.

1. Introduction

Over the last decade, tremendous progress has been made in the field of DNA sequencing [1]. With the introduction of modern next-generation sequencing (NGS) platforms—e.g., Roche 454, Illumina, Ion Torrent and PacBio—the amount of generated sequences has dramatically increased while the costs have rapidly decreased [1,2]. Especially in the field of emerging infectious diseases, whole-genome sequencing using NGS has become an essential and widespread tool in understanding molecular epidemiology and pathogen evolution [3,4].
One emerging animal pathogen for which the need for whole-genome sequences has massively increased in recent years is African swine fever virus (ASFV) [5].
ASFV was first described and is endemic in sub-Saharan Africa [6] where it can be detected in wild African suids—e.g., forest hog, red river hog and bushpigs—and is transmitted in a sylvatic cycle between warthogs and Ornithodoros soft ticks [7]. Since only distantly related giant viruses with double-stranded DNA genomes are known so far, ASFV is the only member of its genus Asfivirus and family Asfarviridae [5]. ASFV causes a disease in domestic pigs and European wild boar called African swine fever (ASF), which is characterised by high fever, haemorrhages, respiratory, gastro-intestinal and neurological signs and cyanosis [8,9]. The animals show a lethality of up to 100% depending on the strain-specific virulence and usually succumb to the infection within 7 to 10 days after infection [9]. Neither vaccines nor treatments are available [10]; thus, stopping ASFV spread relies solely on strict quarantine measures including the culling of infected herds, establishment of restriction zones and movement control leading to enormous socio-economic consequences in affected countries [11,12].
From Africa, ASFV was introduced into Europe (Portugal) two times and subsequently spread to other European countries; e.g., Spain, France, Italy, Belgium and the Netherlands [13,14]. After the virus was eradicated from most affected countries (with the exception of Sardinia, where ASFV is still endemic) [11], ASFV was reintroduced in Georgia in 2007 [15]. Since then, the virus has subsequently spread through wild boar and domestic pig populations of Eastern Europe (Georgia, Armenia, Russia, Belarus, Moldova, Ukraine) [13,16,17,18] and the European Union (Bulgaria, Hungary, Czech Republic, Slovakia, Romania, Estonia, Latvia, Lithuania and Poland) [18]. More recently, the virus spread to Western Europe (Belgium) [19] and reached Asia, causing devastating outbreaks (China, Cambodia, Mongolia, Vietnam, North Korea, Laos and Myanmar) [13,20].
The sequencing of partial genomes by Sanger sequencing has been used for decades to identify viral genotypes and trace virus introductions through molecular epidemiology, especially in Africa where ASFV is widely distributed and different genotypes occur [21,22,23,24,25,26,27]. The viral genome consists of one linear molecule of double-stranded DNA with a length of 170–194 kbp [28]. It shows a low mutation rate typical for DNA viruses that leads to a very low genetic variability, especially in close geographic regions [23,29,30]. However, the identification of attenuated viruses with large deletions, inversions or duplications as well as novel genotypes in Africa might be explained by mechanisms leading to genetic reorganisations such as homologous recombination, as was suggested and described for other large dsDNA viruses [31]. Therefore, the analysis of partial sequences quickly becomes ineffective for detailed analyses of phylogeny, molecular epidemiology, or virus evolution, as observed especially in Eastern Europe [19].
The first ASFV whole-genome sequence was published in 1995 and generated from the Spanish cell-culture-adapted ASFV strain BA71V using Sanger sequencing [32,33]. Since then, only a few additional genomes have been published using the same technique [34,35] (Table 1). In 2009, the first ASFV genome sequence was published using the Roche 454 GS FLX NGS technique in combination with Sanger sequencing [36]. Subsequently, different platforms have been used to generate ASFV whole-genome sequences including Illumina HiScanSQ, MiSeq, HiSeq, NextSeq500 and PacBio [37,38,39,40,41,42,43,44,45,46] (Table 1).
Since the ASFV genome includes extensive homopolymer and repeat regions including inverted terminal repeats (ITR) of variable length [28], platform-specific limitations especially in homopolymer sequencing need to be considered [47,48]. In particular, the short read data generated by sequencing platforms of the second generation—e.g., Roche 454, Illumina and IonTorrent—need to be handled with care. Used for mapping against a reference sequence, genomic reorganisations such as inversions or duplications might be missed, and the quality of the consensus-sequence strongly relies on the sequence used as reference. Furthermore, when used for assembly, the small reads (200–300 bp) might have only small overlapping regions and therefore can be misassembled and—especially in homopolymer and repeat regions—the quality of the consensus-sequence might be low.
Therefore, data on sequence quality such as coverage, mappings and alignments are crucial to assess the results of whole genome sequencing and the analysis of genetic variations including single nucleotide variants (SNV) analyses and molecular epidemiology experiments [3].
In particular, reference sequences frequently used as a basis for experiments and analyses [45,49,50] need careful validation. Especially for ASFV, where only little data on gene expression and protein translation are available and most of the between 150 and 167 ORFs have only been predicted [28,37], reliable whole genome sequences are the basis needed for transcriptome and proteome analyses and up-to-date annotations.
Here, we describe a sequencing protocol combining ASFV-specific target enrichment, Illumina sequencing and long-read Nanopore sequencing to generate high-quality complete ASFV genome sequences. We used the described protocol for resequencing of the ASFV Georgia 2007/1 genome, leading to the correction of 71 homopolymer errors and other sequencing artefacts, a longer ITR region and ultimately updated annotations [37] of the corrected genome sequence. Using this high-quality sequence as a reference, we used target enrichment in combination with Illumina sequencing and sequenced ASFV directly from organ tissue, thereby providing an ASFV whole-genome sequence from an outbreak in Moldova.

2. Materials and Methods

2.1. Virus Cultivation

The original ASFV Georgia 2007/1 isolate from Georgia reported in 2011 [37] was passaged once on porcine bone marrow cells (PBMs) prior to sequencing.
To generate high molecular weight (HMW)-DNA to facilitate long reads through Nanopore sequencing, ASFV Georgia 2007/1 was passaged a second time on porcine primary peripheral blood monocytes (PBMCs) as described elsewhere [54].

2.2. DNA Extraction

For Illumina sequencing, DNA was extracted from PBMC cell culture supernatant using the NucleoMag® VET Kit (Macherey Nagel, Düren, Germany) according to the manufacturer’s instructions.
The Moldovan organ sample (spleen), collected from infected domestic pigs, was homogenised in a 2 mL reaction tube containing one 5 mm steel bead and 200 µL phosphate buffered saline (PBS) in a TissueLyser II (Qiagen, Hilden, Germany) for 3 min at 30 Hz. After centrifugation at 10,000× g for 3 min, DNA was extracted from the supernatant as described above.
For Nanopore sequencing, HMW-DNA was extracted using a customised protocol. PBMCs [54] were infected with ASFV Georgia 2007/1 at an MOI of 0.1 and harvested after 24 h. After one freeze-thaw-cycle of infected PBMCs, cell debris was removed by centrifugation at 1000× g for 10 min at 4 °C and the virus particles were concentrated using ultracentrifugation through a sucrose cushion. Briefly, 6 mL 40% sucrose (in PBS) were carefully overlaid with the virus-containing supernatant (30 mL) and centrifuged for one hour at 50,000× g at 4 °C. The resulting pellet was resuspended in 100 µL PBS and heat-inactivated for 30 min at 75 °C. Subsequently, the virus suspension was mixed with 3 mL TEN buffer (40 mM Tris-HCl, 1 mM EDTA pH 8.0, 150 mM NaCl) and 1.5 mL sarcosyl buffer (3% sodium lauroyl sarcosinate, 75 mM Tris-HCl pH 8.0, 25 mM EDTA) and incubated at 65 °C for 15 min. The DNA was extracted using equal volumes of phenol (Roth) in a first step, phenol/chloroform/isoamyl alcohol (Roth) in a second step and finally chloroform (Roth)/isoamyl alcohol (Merck) (24:1) in a final step. Precipitation of the DNA was carried out using 0.1 volumes of sodium acetate buffer (3 M, pH 4.8) and 2.5 volumes of ethanol at −20 °C for 16 h. The precipitated DNA was pelleted and washed with 70% ethanol and air-dried before final resuspension in 70 µL EB buffer (Qiagen).

2.3. Illumina Sequencing

Samples were prepared for and analysed on an Illumina MiSeq platform in the 300 bp paired-end mode, as previously described [55].

2.4. Nanopore Sequencing

Library preparation was carried out using the Rapid Barcoding Sequencing Kit (SQK-RBK004, Oxford Nanopore Technologies, Oxford, UK) according to the manufacturer’s instructions. The finished library was loaded onto a R9.5.1 MinION flow cell (FLO-MIN107, Oxford Nanopore Technologies). Sequencing was conducted under the standard settings using a MinION (Mk1B, Oxford Nanopore Technologies) in combination with a MinIT (MNT-001, Oxford Nanopore Technologies). Signal information (fast5) was basecalled with MinIT and Guppy v2.1.3 (Oxford Nanopore Technologies).

2.5. Target Enrichment

For the enrichment of ASFV-specific sequences prior to Illumina MiSeq sequencing, the myBaits® kit with 1–20K unique baits (Arbor Bioscience, Ann Arbor, MI, USA) was used according to the manufacturer’s instructions (Figure 1). Briefly, RNA baits (11,958 unique baits) were designed in silico to cover the ASFV Georgia 2007/1 sequence (FR682468.1) and the ASFV Ken06.bus sequence (KM111295) with a three-fold coverage by Arbor Bioscience. After DNA extraction and standard library preparation [55], 100–500 ng of indexed libraries were heat-denatured at 95 °C and, after cooling to hybridisation temperature (65 °C), adapter sequences were blocked using Illumina-specific blocking oligos supplied by the kit (Arbor Bioscience). Denatured and blocked libraries were combined with the pre-heated hybridisation mix containing the ASFV-specific biotinylated RNA-baits and incubated at the pre-set hybridisation temperature (65 °C) for 16 h for bait to target hybridisation. Subsequently, the hybridised bait–target mix was mixed with MyOne™ Streptavidin C1 beads (Thermo Fisher Scientific, Waltham, MA, USA). After immobilisation on a magnetic rack and washing using pre-heated washing buffer (65 °C), target sequences were eluted from the beads by adding 10 mM Tris-HCl, 0.05% TWEEN®-20 solution (pH 8.0–8.5) and incubating for 5 min at 95 °C. Enriched libraries were amplified using the AccuPrime™ Taq DNA Polymerase System (Thermo Fisher Scientific) according to the manufacturer’s instructions for 14 cycles prior to Illumina MiSeq sequencing [55]. Sequencing was performed using a MiSeq v3 600 cycle kit in 2 × 300 bp mode.

2.6. Data Analysis

Sequence data from the Illumina libraries lib02645 and lib02679 were mapped against the ASFV Georgia 2007/1 genome sequence (FR682468.1) using the Newbler 3.0 Software (Roche, Basel, Switzerland), and mapped reads were assembled using SPAdes 3.11.1 [53] in the mode of read error correction prior to assembly.
Basecalled fastq data from Nanopore sequencing were assembled in a hybrid assembly with the Illumina reads mapped to FR682468.1 (see above) using SPAdes v3.10.0. Contigs were used as a further reference for an iterative mapping approach using the generic mapper in Geneious v11.1.5 (Biomatters, Auckland, New Zealand). Regions with a coverage lower than 2.5 times standard deviations from the mean were excluded from consensus generation and further analysis. Nanopore fastq reads were used in a mapping approach in parallel using KMA [56], applying pre-sets optimised for Nanopore reads sensitive for indels. Coding sequences (CDS) were predicted with Glimmer3 [52] using FR682468.1 (ASFV Georgia 2007/1) for model generation. CDS were manually annotated and used to update INSDC entry FR682468.1, now available under FR6824681.2.
Single nucleotide variants (SNVs) were determined using the generic SNP finder of the Geneious software suite, applying parameters of a maximum p-value of 10−5 and filter for strand bias. The threshold for SNP identification was set at 10% and a minimum coverage of 100 was defined. Variants were checked manually for accuracy. Due to the low read quality at the genome ends, variants within the ITR regions were not included in the analysis. Read duplicates from all data used for SNV analysis and coverage estimation were removed using Samtools v. 1.9 [57] prior to analysis.

2.7. PCR and Sanger Sequencing

For conventional PCR, Phire Green Hot Start II PCR Master Mix (Invitrogen, Carlsbad CA USA) and the primers 457F CGGGCCAGACAAAATTGACC and 1223R AACAGGAAATACAAGGCGGC were used in a peqSTAR Thermocycler (VWR, Darmstadt, Germany) with the following cycling conditions: 30 s at 98 °C, 35 cycles each of 10 s at 98 °C, 30 s at 64 °C, 2 min at 72 °C and final elongation 10 min at 72 °C. Corresponding bands (751 bp) were purified from a 1.5% agarose gel with the QIAquick Gel Extraction Kit (Qiagen) and used for sequencing reactions (BigDye Terminator v3.1 Cycle Sequencing Kit, Applied Biosystems, Darmstadt, Germany). The reaction products were purified using NucleoSEQ columns (Macherey–Nagel) and sequenced on an ABI PRISM® 3100 Genetic Analyzer (Life Technologies, Darmstadt, Germany). Obtained sequences were assembled using the Geneious software, v11.1.5.

2.8. Data Availability

Whole-genome sequences are available from the International Nucleotide Sequence Database Collaboration (INSDC) databases under the study accession number PRJEB33279, and under the GenBank accession number FR682468.2.

3. Results

3.1. ASFV-Specific Target Enrichment Prior to Illumina Sequencing Provides High Amounts of Target Reads

For the resequencing of ASFV Georgia 2007/1, we compared a shotgun and a target enrichment approach. For shotgun sequencing, two libraries were prepared and sequenced; Table 2 summarises the key figures of these sequencing efforts. Library lib02645 was sequenced twice on the Illumina MiSeq platform. The first run resulted in 1,764,078 reads, from which 8438 reads (0.48%) were identified as ASFV-specific by mapping and resulted in a whole-genome sequence. However, the mean coverage of 13 was too low to assemble a reliable whole-genome sequence, especially in repeat and homopolymer regions. Therefore, in a second run, we sequenced 7,316,924 total reads, from which 36,670 reads were identified as ASFV-specific by mapping (0.5%), resulting in a whole-genome sequence with a mean coverage of 57 (Table 1) and leading to a better resolution.
Using target enrichment on the shotgun lib02679 prior to sequencing, we were able to increase the amount of ASFV reads identified by mapping to 44,862 from 60,174 total reads (66.78%), and whole-genome assembly was successful with a mean coverage of 57.2 (Table 1), also leading to a good resolution combined with less sequencing capacity.
However, the assembly of the resulting data was complicated by the extensive repeat regions in the ASFV genome as shown by the drop in coverage in the specific regions (Figure 2).

3.2. Hybrid Assembly of Nanopore and Illumina Data Provide an Improved ASFV Georgia 2007/1 Whole-Genome Sequence with Novel Information about the Genome Length

To overcome the limitations of short reads and to analyse the sequence of the inverted terminal repeat regions, long reads were generated by sequencing a library (lib02923) prepared from ASFV Georgia 2007/1 infected PBMCs using Nanopore sequencing. A total of 136,627 reads was generated with an estimated base calling accuracy of around 88.8%, and 1648 ASFV-specific Nanopore reads (1.20%) were identified by mapping and used for further analysis. The mean length of the used reads was 7881.4 nt, the median read length was 5517 nt, the N50 read length was 13,165 nt, with the longest read spanning was 43,474 nt (Figure S1). ASFV-specific Nanopore reads were used in a hybrid assembly with the ASFV-specific Illumina reads (lib02645, lib02679). The resulting contig was used as reference for further mapping and assembly, resulting in a new ASFV Georgia 2007/1 (FR682468.2) sequence with an improved overall coverage, particularly in repeat regions (Figure 2), and a length of 190,594 bp.
With a mean coverage of 67.5-fold, Nanopore reads were found to cover 100% of the ASFV Georgia 2007/1 reference sequence (Figure 2 and Figure S1) using a method not adapted to the Nanopore read characteristics. A mapping approach using only Nanopore reads and applying an algorithm with parameters adapted to the specific features of these data (KMA see methods) gained a consensus with only four differences (99.9% identical) to the reference sequence, with differences occurring in only one large homopolymer region (17xC in reference; 14xC in gained consensus), underlining the potential of this technique to support the assembly of large virus genomes and pointing to its constraints to dissolve homopolymers.

3.3. Sequence Alignment with the Improved ASFV Genome Reveals 71 Differences in Homopolymer Regions and Open Reading Frames

After gene predictions using Glimmer 3 in Geneious, the improved ASFV Georgia 2007/1 sequence was annotated on the basis of the available ASFV Georgia 2007/1 sequence. The annotation was manually checked and modified to accommodate changes in ORFs; e.g., split, fusion or truncations. Modified annotations include the ORFs (including small ORFs and ORFs within larger ORFs): MGF360 1L, G ACD 00070, G ACD 00160, 285L, MGF 110-10L, MGF110-14L, MGF 110-13L, G ACD 00290, MGF 300-2R, G ACD 01760, DP63R, MGF100-3L and MGF 360-19R (Table S1).
Through the pairwise alignment of the original ASFV Georgia 2007/1 sequence (FR682468.1) with the improved ASFV Georgia 2007/1 sequence (FR682468.2), we identified 71 differences in the updated genome sequence leading to an overall sequence identity of 99.851%. From these differences, 58 are insertions or deletions (indels) in A/T homopolymer regions, seven are indels in G/C homopolymer regions, four are single nucleotide insertions and two are transversions. From these 71 differences, 22 are located in predicted and annotated open reading frames (ORF) causing frame shifts, while 49 do not affect ORFs. From the 22 frameshifts, eight lead to a shortened ORF, four to enlargement, three to the fusion of two ORFs, two to the split into two ORFs and five to the variation of codons in an ORF (Table S1). Altogether, 13 changes affect ORFs previously shown to be expressed including the CP204L and KP177R ORFs coding for the structural proteins P30 (CP204L) and P22 (KP177R), A151R, C84L, QP383R, three members of the MGFs 110 and 360 and the small ORFs ASFV G ACD 00290, 350, 01760 (Table S1). These differences to ASFV Georgia 2007/1 (FR682468.1) can also be identified in the published whole-genome sequences of ASFV Estonia 2014 [43], ASFV-POL/2015/Podlaskie [41] and ASFV Belgium 2018/1 [53] and are in agreement with differences to ASFV Georgia 2007/1 described in ASFV China [52], ASFV Georgia 2008 [45] and ASFV Poland [51].
As another important result, we were able to expand the 5’-ITR region of ASFV Georgia 2007/1 by 956 bp and the 3’-ITR region by 231 bp leading to a total length of 1378 bp of each ITR region, the annotation of the ORFs DP60R and ASFV G ACD 01990 at both ITR regions and a total genome length of 190,594 bp (Table S1).
For validation of the in silico assemblies, the junction at the 5’ ITR was confirmed by PCR and Sanger sequencing.

3.4. Application of Target Enrichment Prior to Illumina Sequencing Enabled Whole-Genome Assembly of ASFV Moldova from Organ Samples Using the New Improved Sequence as Reference

Sequencing spleen samples from an ASFV infected pig from Moldova using a shotgun sequencing approach, we generated 8,232,518 reads, from which only 0.05% were ASFV-specific. However, the total amount of 4042 reads was not sufficient to generate a whole-genome sequence, neither by mapping nor by assembly (Table 1). By using the established target enrichment protocol on a shotgun library prior to sequencing, we were able to increase the amount of ASFV-specific reads to 690,206 (83.89%), from which 207,763 (25.0%) were unique. Mapping these data against the improved reference sequence (ASFV Georgia 2007/1—FR682468.2) and using SPAdes for the assembly, we were able to generate a whole-genome sequence of ASFV Moldova 2017/1 with a mean unique coverage of 317 (Table 1).

3.5. Variant Analysis Reveals Possible Single Nucleotide Polymorphisms

By mapping all data from ASFV Georgia 2007/1 libraries against the improved ASFV Georgia 2007/1 consensus sequence, we identified 49 locations with nucleotide variant frequencies of more than 10%. From these, 30 are located within 21 annotated ASFV-ORFs, leading to 26 amino acid substitutions, one truncation and one extension, while two do not change the corresponding amino acid. While the observed variant frequency is below 15% at 44 positions, we observed frequencies between 17.6% and 29.3% at positions 2963, 12,525 and 176,519. In addition, we observed frequencies in two poly G/C homopolymer regions of 40.8% and 47.7% (Table S3).
For ASFV Moldova 2017/1, we identified 12 locations with nucleotide variant frequencies of more than 10%. From these, eight are located within seven annotated ORFs leading to frameshifts in six, a substitution and corresponding amino acid change in one and a silent mutation in another case (Table S4).

3.6. Alignment of ASFV Moldova 2017/1 and ASFV Georgia 2007/1 Reveals Only Single Nucleotide Differences and a Tandem Repeat Variation

Comparing the ASFV Moldova 2017/1 genome with the improved ASFV Georgia 2007/1 sequence, we identified an overall sequence identity of 99.98% and 11 locations showing differences in single nucleotides (Table S2). From these 11 locations, five were insertions or deletions in A/T homopolymer regions, three nucleotide transversions and three transitions. While four were located in non-coding regions, five are non-synonymous, leading to a codon change in annotated genes. From these, one leads to a stop codon integration causing the truncation of the MGF-110 1L ORF by 18 codons from the 5’-end (Table S2) not affecting the predicted endoplasmic reticulum retention motif [58]. Furthermore, one leads to a frameshift, causing the variation of 23 and addition of five codons of the DP60R ORFs 5’-end as previously described in other whole-genome sequences from Eastern Europe [51] (Table S2). In addition, we identified the integration of a 10 bp repeat in the previously described variable tandem repeat region between I73L and I329L [59] as well as differences in four non-coding poly G/C regions (Table S2). Comparing these differences with all available Eastern European whole-genome sequences, only three of these differences were found exclusively in the genome sequence of ASFV Moldova 2017/1 (Table S2).

3.7. Alignment and Phylogenetic Reconstruction of All Available ASFV Whole-Genome Sequences Shows a Very High Overall Nucleotide Sequence Identity of All Eurasian Strains

The alignment of all available ASFV whole-genome sequences revealed an overall nucleotide sequence identity of all Asian, Western and Eastern European whole genome sequences of over 99.9%. This high degree of identity is also shown by phylogenetic analysis; e.g., the resulting phylogenetic tree (Figure 3). While ASFV Moldova 2017/1 showed the highest nucleotide sequence identity to ASFV Belgium 2018/1 (99.99%), ASFV Georgia 2007/1 was most similar to ASFV Moldova 2017/1 and ASFV_Pol16_29413 _o23_ MG939586 (99.98%).

4. Discussion

In this study, we used a combined workflow including target enrichment, Illumina and Nanopore sequencing (Figure 4) for the resequencing of the ASFV Georgia 2007/1 isolate representing the first ASFV introduced into Eastern Europe in 2007.
Since the ASFV strains circulating in Eastern Europe and Asia are highly similar (Figure 3), molecular epidemiology analyses will rely on small changes and even single nucleotide variations. Therefore, the availability of high-quality sequences and the sharing of raw data as well as data on alignments and mappings is highly desirable to evaluate the quality of a sequence prior to analysis and search for a distinct pattern that could help to draw spatial and temporal conclusions.
With the use of ASFV target-specific enrichment prior to sequencing, we were able to increase the amount of target-specific reads allowing for whole-genome assembly while reducing the total number of reads considerably. As shown by our sequencing results, target-specific enrichment is especially useful when sequencing samples with a low virus to host ratio—e.g., organ tissue (Figure 4)—while samples with higher virus to host ratio—e.g., cell culture supernatant—do not necessarily call for enrichment (Figure 4).
However, since target-specific enrichment is performed using a bait set designed on known sequences, the enrichment and therefore the sequencing is limited to these. Although a hybridisation of baits to target sequences showing up to 90% identity can be expected [58], regions with extensive genetic variations or longer insertions might be missed. Therefore, bait sets should be designed carefully for the corresponding scientific question; e.g., specific sets for the whole-genome assembly of known viral genomes or sets including different genome variants for a broader enrichment. Furthermore, hybridisation parameters—e.g., temperature—can be optimised to include sequence variants [60]. In addition, the bias introduced by the amplification needed to generate the size of the library required for sequencing must be considered, especially when working on samples with a very low amount of target-specific sequences, and corresponding datasets need to be de-duplicated prior to analysis.
Since the ASFV genome includes extensive repeat regions, ITRs of unknown length and homopolymer regions with up to 17 C and 16 G, even data from target enrichment approaches have proven difficult for the de novo assembly of a high-quality reference sequence (Figure 2).
Through the employment of single molecule Nanopore sequencing (Figure 4), we were able to generate long viral reads of up to 43.5 kb. However, the low accuracy of the Nanopore reads (see above) hampers precise consensus sequence generation. While Nanopore reads reach an accuracy of around 90%, Illumina data sets usually have accuracies around 99.9%. Hence, the higher accuracy of the Illumina data sets can be used for better precision, whereas the Nanopore reads can be taken for correct assembly. Therefore, using these in a hybrid assembly, we were able to further improve the sequence quality, especially in the ITR regions. Although we were able to enlarge the ITR regions considerably, because of the decreasing sequence coverage at the ends, we cannot exclude the possibility that they are even longer. Due to their proposed involvement in ASFV replication involving head to head concatemers [61] (similar to poxviruses [62]), a complete sequence of the ITR regions including the terminal hairpin loops [63] would be especially interesting.
In our opinion, the presented workflows offer the best possible quality of ASFV sequences using the technology available to date. The sequences presented here are therefore of a reference nature and improve the analysis of additional ASFV sequences by far. Nonetheless, even with the combination of Nanopore and Illumina reads in a hybrid assembly, the sequencing of the poly G/C homopolymer regions has proven to be difficult, and we can neither completely exclude the existence of viral variants nor sequencing or bioinformatics artefacts, respectively. Therefore, it remains speculative if these homopolymer stretches might be involved in differential gene expression (maybe in the arthropod vector) or other unknown mechanisms, as might be suggested by the observed fusion or split of ORFs depending on the length of the homopolymer regions.
For all sequencing platforms, their specific sequencing limitations always need to be considered. For the 454 technology used for the generation of the original ASFV Georgia 2007/1 sequence (FR682468.1), limitations especially in sequencing homopolymer regions are well known [64]. Therefore, it is not surprising that we discovered numerous differences between the original and the improved sequence in A/T homopolymer regions, as also identified previously [41]. Since even the most recently developed sequencing technologies have their limitations, as shown by the evaluation of the IonTorrent and Illumina systems especially for homopolymer sequencing [64,65,66], selecting the adequate sequencing platform is crucial for each scientific question and should be done carefully. Furthermore, sequence information from the database should be evaluated thoroughly using the available information on the employed sequencing platform and data analysis prior to their use in experimental designs or bioinformatic analyses.
Using the ASFV-specific target enrichment prior to sequencing and the now available improved ASFV Georgia 2007/1 sequence as a reference, we could easily provide a reliable whole-genome sequence for ASFV Moldova 2017/1 from organ material without prior need for virus cultivation (Figure 4B). While the overall nucleotide sequence identity is more than 99.9% similar to the circulating ASFV strains from Eastern Europe and Asia, only a few single differences in nucleotides as well as a previously described tandem repeat insertion could be observed. However, the observed nucleotide differences, together with the available whole-genome sequence information, unfortunately do not allow for any further conclusions regarding phylogenetic or geographic relationships.
Although some of the single nucleotide changes in ASFV Georgia 2007/1 and ASFV Moldova 2017/1 reported in this study affect annotated ORFs, due to the lack of expression data for most of the ASFV ORFs and the lack of observations of virus attenuation from the field, no conclusion regarding their influence on virus replication or pathogenicity can be given. Therefore, transcriptome and proteome studies and infection experiments are needed to validate and update annotations, evaluate the influence of these changes on the viral phenotype and identify essential viral genes that might serve as targets for vaccine development (e.g., epitopes for T- or B-cells).
However, analysing SNVs in the viral genomes, we detected a clear difference in the number of SNVs. While we identified 49 SNVs in ASFV Georgia 2007/1, only 12 were detected in ASFV Moldova 2017/1. Since ASFV Georgia 2017/1 was passaged and sequenced from cell culture while ASFV Moldova 2017/1 was sequenced directly from organ material, it might be hypothesised that SNVs accumulate in cell culture-raised ASFV as was suggested for other viruses [67,68,69]. However, this observation might be influenced by the differences in coverage at the SNV sites, and although ASFV Moldova 2017/1 has a higher coverage, which should lead to a higher number of SNVs and a low sample size, further experiments with cell culture and field samples will be essential to confirm this hypothesis. In addition, more data on SNV in ASFV populations might be used in fine-tuning the intra-genotype phylogeny [70].
In conclusion, we used a workflow including target enrichment and Nanopore sequencing to provide an improved genomic sequence of ASFV Georgia 2007/1 (FR682468.2). Using this sequence as a reference, we subsequently generated the first whole-genome sequence from ASFV Moldova 2017/1, and showed that it is highly similar to the known circulating strains, showing only a few specific differences in single nucleotides from which only three seem to be specific for ASFV Moldova 2017/1.
We believe that our deep-sequencing workflow for ASFV and the provided reference sequences will be highly valuable to generate further high-quality ASFV whole-genome sequences and can serve as a basis for variant, transcriptome and proteome analyses. We are convinced that these sequences, along with coordinated efforts to improve data sharing and harmonised protocols, might pave the way for the identification of the new genetic markers that are desperately needed to understand virus evolution and trace routes of ASFV to eliminate the burden of this devastating animal disease.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/11/9/846/s1, Figure S1: Read length distribution of ASFV-specific ONT reads and mapping positions on the ASFV Georgia 2007/1 genome. Table S1: Comparison of ASFV Georgia 20007/1 FR682468.1 with ASFV Georgia 2007/1 FR682468.2. Table S2: Comparison of ASFV Moldova2017/1 and ASFV Georgia 2007/1 FR682468.2. Table S3: Single nucleotide variant analysis of ASFV Georgia 2007/1. Table S4: Single nucleotide variant analysis of ASFV Moldova 2017/1.

Author Contributions

Conceptualisation, J.H.F., T.B.R., S.B., M.B.; formal analysis, J.H.F., L.F.F., J.K., D.H., T.B.R.; investigation, J.H.F., L.F.F., J.K., A.H., A.S.O., D.H., T.B.R; resources, O.G., L.K.D., C.L.N..; data curation, A.P. writing—original draft preparation, J.H.F, J.K., A.P.; writing—review and editing, J.H.F., A.P., A.S.O, T.B.R., C.L.N., L.K.D., D.H., S.B., M.B.; visualisation, J.H.F., A.P..; supervision, S.B., M.B., T.B.R.

Funding

This work was funded through the FLI’s ASFV-Research Network and in part by COMPARE under the European Union’s Horizon 2020 program under grant agreement No. 643476. This publication was supported by the European Virus Archive goes Global (EVAg) project that has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 653316. L.K.D. and C.L.N. were supported by Biotechnology and Biological Sciences Research Council (BBSRC; https://bbsrc.ukri.org/) under grants BBS/E/1/00007031 and BBS/E/1/00007034.

Acknowledgments

The authors are very thankful to Patrick Zitzow for excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goodwin, S.; McPherson, J.D.; McCombie, W.R. Coming of age: Ten years of next-generation sequencing technologies. Nat. Rev. Genet. 2016, 17, 333–351. [Google Scholar] [CrossRef] [PubMed]
  2. Levy, S.E.; Myers, R.M. Advancements in Next-Generation Sequencing. Annu. Rev. Genom. Hum. Genet. 2016, 17, 95–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bayliss, S.C.; Verner-Jeffreys, D.W.; Bartie, K.L.; Aanensen, D.M.; Sheppard, S.K.; Adams, A.; Feil, E.J. The Promise of Whole Genome Pathogen Sequencing for the Molecular Epidemiology of Emerging Aquaculture Pathogens. Front. Microbiol. 2017, 8, 121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tang, P.; Croxen, M.A.; Hasan, M.R.; Hsiao, W.W.; Hoang, L.M. Infection control in the new age of genomic epidemiology. Am. J. Infect. Control 2017, 45, 170–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Alonso, C.; Borca, M.; Dixon, L.; Revilla, Y.; Rodriguez, F.; Escribano, J.M.; Ictv Report, C. ICTV Virus Taxonomy Profile: Asfarviridae. J. Gen. Virol. 2018, 99, 613–614. [Google Scholar] [CrossRef]
  6. Montgomery, R.E. On A Form of Swine Fever Occurring in British East Africa (Kenya Colony). J. Comp. Pathol. 1921, 34, 159–191. [Google Scholar] [CrossRef] [Green Version]
  7. Burrage, T.G. African swine fever virus infection in Ornithodoros ticks. Virus Res. 2013, 173, 131–139. [Google Scholar] [CrossRef]
  8. Blome, S.; Gabriel, C.; Dietze, K.; Breithaupt, A.; Beer, M. High virulence of African swine fever virus caucasus isolate in European wild boars of all ages. Emerg. Infect. Dis. 2012, 18, 708. [Google Scholar] [CrossRef]
  9. Blome, S.; Gabriel, C.; Beer, M. Pathogenesis of African swine fever in domestic pigs and European wild boar. Virus Res. 2013, 173, 122–130. [Google Scholar] [CrossRef]
  10. Arias, M.; de la Torre, A.; Dixon, L.; Gallardo, C.; Jori, F.; Laddomada, A.; Martins, C.; Parkhouse, R.M.; Revilla, Y.; Rodriguez, F.A.J.; et al. Approaches and Perspectives for Development of African Swine Fever Virus Vaccines. Vaccines 2017, 5, 35. [Google Scholar] [CrossRef]
  11. Arias, M.; Sánchez-Vizcaíno, J.M. African Swine Fever Eradication: The Spanish Model. In Trends in Emerging Viral Infections of Swine, 1st ed.; Morilla, A., Yoon, K.J., Zimmerman, J.J., Eds.; Iowa State Press: Iowa City, IA, USA, 2002; pp. 133–139. [Google Scholar]
  12. Sanchez-Cordon, P.J.; Montoya, M.; Reis, A.L.; Dixon, L.K. African swine fever: A re-emerging viral disease threatening the global pig industry. Vet. J. 2018, 233, 41–48. [Google Scholar] [CrossRef] [PubMed]
  13. Sanchez-Vizcaino, J.M.; Mur, L.; Gomez-Villamandos, J.C.; Carrasco, L. An update on the epidemiology and pathology of African swine fever. J. Comp. Pathol. 2015, 152, 9–21. [Google Scholar] [CrossRef] [PubMed]
  14. Sanchez-Vizcaino, J.M.; Mur, L.; Martinez-Lopez, B. African swine fever: An epidemiological update. Transbound. Emerg. Dis. 2012, 59, 27–35. [Google Scholar] [CrossRef] [PubMed]
  15. Beltran-Alcrudo, D.; Lubroth, J.; Depner, K.; De La Rocque, S. African swine fever in the Caucasus. FAO EMPRES Watch 2008, 1, 1–8. [Google Scholar]
  16. Gogin, A.; Gerasimov, V.; Malogolovkin, A.; Kolbasov, D. African swine fever in the North Caucasus region and the Russian Federation in years 2007–2012. Virus Res. 2013, 173, 198–203. [Google Scholar] [CrossRef] [PubMed]
  17. Kolbasov, D.; Titov, I.; Tsybanov, S.; Gogin, A.; Malogolovkin, A. African Swine Fever Virus, Siberia, Russia, 2017. Emerg. Infect. Dis. 2018, 24, 796–798. [Google Scholar] [CrossRef] [PubMed]
  18. Chenais, E.; Stahl, K.; Guberti, V.; Depner, K. Identification of Wild Boar-Habitat Epidemiologic Cycle in African Swine Fever Epizootic. Emerg. Infect. Dis. 2018, 24, 810–812. [Google Scholar] [CrossRef]
  19. Garigliany, M.; Desmecht, D.; Tignon, M.; Cassart, D.; Lesenfant, C.; Paternostre, J.; Volpe, R.; Cay, A.B.; van den Berg, T.; Linden, A. Phylogeographic Analysis of African Swine Fever Virus, Western Europe, 2018. Emerg. Infect. Dis. 2019, 25, 184–186. [Google Scholar] [CrossRef]
  20. Zhou, X.; Li, N.; Luo, Y.; Liu, Y.; Miao, F.; Chen, T.; Zhang, S.; Cao, P.; Li, X.; Tian, K.; et al. Emergence of African Swine Fever in China, 2018. Transbound. Emerg. Dis. 2018, 65, 1482–1484. [Google Scholar] [CrossRef] [Green Version]
  21. Luka, P.D.; Achenbach, J.E.; Mwiine, F.N.; Lamien, C.E.; Shamaki, D.; Unger, H.; Erume, J. Genetic Characterization of Circulating African Swine Fever Viruses in Nigeria (2007–2015). Transbound. Emerg. Dis. 2017, 64, 1598–1609. [Google Scholar] [CrossRef]
  22. Luka, P.D.; Erume, J.; Mwiine, F.N.; Shamaki, D.; Yakubu, B. Comparative Sequence Analysis of Different Strains of African Swine Fever Virus Outer Proteins Encoding Genes from Nigeria, 2009–2014. J. Virol. Mol. Biol. 2016, 5, 16–26. [Google Scholar]
  23. Alkhamis, M.A.; Gallardo, C.; Jurado, C.; Soler, A.; Arias, M.; Sanchez-Vizcaino, J.M. Phylodynamics and evolutionary epidemiology of African swine fever p72-CVR genes in Eurasia and Africa. PLoS ONE 2018, 13, e0192565. [Google Scholar] [CrossRef] [PubMed]
  24. Gallardo, C.; Fernandez-Pinero, J.; Pelayo, V.; Gazaev, I.; Markowska-Daniel, I.; Pridotkas, G.; Nieto, R.; Fernandez-Pacheco, P.; Bokhan, S.; Nevolko, O.; et al. Genetic variation among African swine fever genotype II viruses, eastern and central Europe. Emerg. Infect. Dis. 2014, 20, 1544–1547. [Google Scholar] [CrossRef] [PubMed]
  25. Gallardo, C.; Mwaengo, D.M.; Macharia, J.M.; Arias, M.; Taracha, E.A.; Soler, A.; Okoth, E.; Martin, E.; Kasiti, J.; Bishop, R.P. Enhanced discrimination of African swine fever virus isolates through nucleotide sequencing of the p54, p72, and pB602L (CVR) genes. Virus Genes 2009, 38, 85–95. [Google Scholar] [CrossRef] [PubMed]
  26. Nix, R.J.; Gallardo, C.; Hutchings, G.; Blanco, E.; Dixon, L.K. Molecular epidemiology of African swine fever virus studied by analysis of four variable genome regions. Arch. Virol. 2006, 151, 2475–2494. [Google Scholar] [CrossRef]
  27. Onzere, C.K.; Bastos, A.D.; Okoth, E.A.; Lichoti, J.K.; Bochere, E.N.; Owido, M.G.; Ndambuki, G.; Bronsvoort, M.; Bishop, R.P. Multi-locus sequence typing of African swine fever viruses from endemic regions of Kenya and Eastern Uganda (2011–2013) reveals rapid B602L central variable region evolution. Virus Genes 2017, 54, 111–123. [Google Scholar] [CrossRef] [PubMed]
  28. Dixon, L.K.; Chapman, D.A.; Netherton, C.L.; Upton, C. African swine fever virus replication and genomics. Virus Res. 2013, 173, 3–14. [Google Scholar] [CrossRef]
  29. Michaud, V.; Randriamparany, T.; Albina, E. Comprehensive phylogenetic reconstructions of African swine fever virus: Proposal for a new classification and molecular dating of the virus. PLoS ONE 2013, 8, e69662. [Google Scholar] [CrossRef]
  30. Dixon, L.K.; Wilkinson, P.J. Genetic Diversity of African Swine Fever Virus Isolates from Soft Ticks (Ornithodoros moubata) Inhabiting Warthog Burrows in Zambia. J. Gen. Virol. 1988, 69, 2981–2993. [Google Scholar] [CrossRef]
  31. Qin, L.; Evans, D.H. Genome scale patterns of recombination between coinfecting vaccinia viruses. J. Virol. 2014, 88, 5277–5286. [Google Scholar] [CrossRef]
  32. Yanez, R.J.; Rodriguez, J.M.; Nogal, M.L.; Yuste, L.; Enriquez, C.; Rodriguez, J.F.; Vinuela, E. Analysis of the complete nucleotide sequence of African swine fever virus. Virology 1995, 208, 249–278. [Google Scholar] [CrossRef] [PubMed]
  33. Sanger, F.; Coulson, A.R. A rapid method for determining sequences in DNA by primed synthesis with DNA polymerase. J. Mol. Biol. 1975, 94, 441–448. [Google Scholar] [CrossRef]
  34. Haresnape, J.M.; Wilkinson, P.J. A study of African swine fever virus infected ticks (Ornithodoros moubata) collected from three villages in the ASF enzootic area of Malawi following an outbreak of the disease in domestic pigs. Epidemiol. Infect. 2009, 102, 507. [Google Scholar] [CrossRef] [PubMed]
  35. Chapman, D.A.; Tcherepanov, V.; Upton, C.; Dixon, L.K. Comparison of the genome sequences of non-pathogenic and pathogenic African swine fever virus isolates. J. Gen. Virol. 2008, 89, 397–408. [Google Scholar] [CrossRef] [PubMed]
  36. de Villiers, E.P.; Gallardo, C.; Arias, M.; da Silva, M.; Upton, C.; Martin, R.; Bishop, R.P. Phylogenomic analysis of 11 complete African swine fever virus genome sequences. Virology 2010, 400, 128–136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Chapman, D.A.; Darby, A.C.; Da Silva, M.; Upton, C.; Radford, A.D.; Dixon, L.K. Genomic analysis of highly virulent Georgia 2007/1 isolate of African swine fever virus. Emerg. Infect. Dis. 2011, 17, 599–605. [Google Scholar] [CrossRef] [PubMed]
  38. Portugal, R.; Coelho, J.; Höper, D.; Little, N.S.; Smithson, C.; Upton, C.; Martins, C.; Leitao, A.; Keil, G.M. Related strains of African swine fever virus with different virulence: Genome comparison and analysis. J. Gen. Virol. 2015, 96, 408–419. [Google Scholar] [CrossRef] [PubMed]
  39. Bishop, R.P.; Fleischauer, C.; de Villiers, E.P.; Okoth, E.A.; Arias, M.; Gallardo, C.; Upton, C. Comparative analysis of the complete genome sequences of Kenyan African swine fever virus isolates within p72 genotypes IX and X. Virus Genes 2015, 50, 303–309. [Google Scholar] [CrossRef]
  40. Granberg, F.; Torresi, C.; Oggiano, A.; Malmberg, M.; Iscaro, C.; De Mia, G.M.; Belak, S. Complete Genome Sequence of an African Swine Fever Virus Isolate from Sardinia, Italy. Genome Announc. 2016, 4, e01220-16. [Google Scholar] [CrossRef] [Green Version]
  41. Olesen, A.S.; Lohse, L.; Dalgaard, M.D.; Wozniakowski, G.; Belsham, G.J.; Botner, A.; Rasmussen, T.B. Complete genome sequence of an African swine fever virus (ASFV POL/2015/Podlaskie) determined directly from pig erythrocyte-associated nucleic acid. J. Virol. Methods 2018, 261, 14–16. [Google Scholar] [CrossRef] [Green Version]
  42. Masembe, C.; Sreenu, V.B.; Da Silva Filipe, A.; Wilkie, G.S.; Ogweng, P.; Mayega, F.J.; Muwanika, V.B.; Biek, R.; Palmarini, M.; Davison, A.J. Genome Sequences of Five African Swine Fever Virus Genotype IX Isolates from Domestic Pigs in Uganda. Microbiol. Resour. Announc. 2018, 7, e01018-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zani, L.; Forth, J.H.; Forth, L.; Nurmoja, I.; Leidenberger, S.; Henke, J.; Carlson, J.; Breidenstein, C.; Viltrop, A.; Höper, D.; et al. Deletion at the 5′-end of Estonian ASFV strains associated with an attenuated phenotype. Sci. Rep. 2018, 8, 6510. [Google Scholar] [CrossRef] [PubMed]
  44. Rodriguez, J.M.; Moreno, L.T.; Alejo, A.; Lacasta, A.; Rodriguez, F.; Salas, M.L. Genome Sequence of African Swine Fever Virus BA71, the Virulent Parental Strain of the Nonpathogenic and Tissue-Culture Adapted BA71V. PLoS ONE 2015, 10, e0142889. [Google Scholar] [CrossRef] [PubMed]
  45. Farlow, J.; Donduashvili, M.; Kokhreidze, M.; Kotorashvili, A.; Vepkhvadze, N.G.; Kotaria, N.; Gulbani, A. Intra-epidemic genome variation in highly pathogenic African swine fever virus (ASFV) from the country of Georgia. Virol. J. 2018, 15, 190. [Google Scholar] [CrossRef] [PubMed]
  46. Bacciu, D.; Deligios, M.; Sanna, G.; Madrau, M.P.; Sanna, M.L.; Dei Giudici, S.; Oggiano, A. Genomic analysis of Sardinian 26544/OG10 isolate of African swine fever virus. Virol. Rep. 2016, 6, 81–89. [Google Scholar] [CrossRef] [Green Version]
  47. Quince, C.; Lanzen, A.; Curtis, T.P.; Davenport, R.J.; Hall, N.; Head, I.M.; Read, L.F.; Sloan, W.T. Accurate determination of microbial diversity from 454 pyrosequencing data. Nat. Methods 2009, 6, 639–641. [Google Scholar] [CrossRef] [PubMed]
  48. Luo, C.; Tsementzi, D.; Kyrpides, N.; Read, T.; Konstantinidis, K.T. Direct comparisons of Illumina vs. Roche 454 sequencing technologies on the same microbial community DNA sample. PLoS ONE 2012, 7, e30087. [Google Scholar]
  49. Borca, M.V.; O’Donnell, V.; Holinka, L.G.; Ramirez-Medina, E.; Clark, B.A.; Vuono, E.A.; Berggren, K.; Alfano, M.; Carey, L.B.; Richt, J.A.; et al. The L83L ORF of African swine fever virus strain Georgia encodes for a non-essential gene that interacts with the host protein IL-1beta. Virus Res. 2018, 249, 116–123. [Google Scholar] [CrossRef] [PubMed]
  50. Borca, M.V.; O’Donnell, V.; Holinka, L.G.; Sanford, B.; Azzinaro, P.A.; Risatti, G.R.; Gladue, D.P. Development of a fluorescent ASFV strain that retains the ability to cause disease in swine. Sci. Rep. 2017, 7, 46747. [Google Scholar] [CrossRef] [PubMed]
  51. Mazur-Panasiuk, N.; Wozniakowski, G.; Niemczuk, K. The first complete genomic sequences of African swine fever virus isolated in Poland. Sci. Rep. 2019, 9, 4556. [Google Scholar] [CrossRef]
  52. Bao, J.; Wang, Q.; Lin, P.; Liu, C.; Li, L.; Wu, X.; Chi, T.; Xu, T.; Ge, S.; Liu, Y.; et al. Genome comparison of African swine fever virus China/2018/AnhuiXCGQ strain and related European p72 Genotype II strains. Transbound. Emerg. Dis. 2019, 66, 1167–1176. [Google Scholar] [CrossRef] [PubMed]
  53. Forth, J.H.; Tignon, M.; Cay, A.B.; Forth, L.F.; Höper, D.; Blome, S.; Beer, M. Comparative Analysis of Whole-Genome Sequence of African Swine Fever Virus Belgium 2018/1. Emerg. Infect. Dis. 2019, 25, 1249. [Google Scholar] [CrossRef] [PubMed]
  54. Carrascosa, A.L.; Bustos, M.J.; de Leon, P. Methods for growing and titrating African swine fever virus: Field and laboratory samples. Curr. Protoc. Cell Biol. 2011, 53, 1–25. [Google Scholar]
  55. Wylezich, C.; Papa, A.; Beer, M.; Höper, D. A Versatile Sample Processing Workflow for Metagenomic Pathogen Detection. Sci. Rep. 2018, 8, 13108. [Google Scholar] [CrossRef] [PubMed]
  56. Clausen, P.; Aarestrup, F.M.; Lund, O. Rapid and precise alignment of raw reads against redundant databases with KMA. BMC Bioinform. 2018, 19, 307. [Google Scholar] [CrossRef] [PubMed]
  57. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Netherton, C.; Rouiller, I.; Wileman, T. The Subcellular Distribution of Multigene Family 110 Proteins of African Swine Fever Virus Is Determined by Differences in C-Terminal KDEL Endoplasmic Reticulum Retention Motifs. J. Virol. 2004, 78, 3710–3721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Goller, K.V.; Malogolovkin, A.S.; Katorkin, S.; Kolbasov, D.; Titov, I.; Höper, D.; Beer, M.; Keil, G.M.; Portugal, R.; Blome, S. Tandem Repeat Insertion in African Swine Fever Virus, Russia, 2012. Emerg. Infect. Dis. 2015, 21, 731–732. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Arbor_Bioscience Arbor Bioscience myBaits® Custom Target Capture Kits Homepage. Available online: https://arborbiosci.com/products/custom-target-capture/ (accessed on 8 April 2019).
  61. Rojo, G.; García-Beato, R.; Viñuela, E.; Salas, M.L.; Salas, J. Replication of African Swine Fever Virus DNA in Infected Cells. Virology 1999, 257, 524–536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Du, S.; Traktman, P. Vaccinia virus DNA replication: Two hundred base pairs of telomeric sequence confer optimal replication efficiency on minichromosome templates. Proc. Natl. Acad. Sci. USA 1996, 93, 9693–9698. [Google Scholar] [CrossRef]
  63. Gonzalez, A.; Talavera, A.; Almendral, J.M.; Vinuela, E. Hairpin Loop Structure of African Swine Fever Virus-DNA. Nucleic Acids Res. 1986, 14, 6835–6844. [Google Scholar] [CrossRef] [PubMed]
  64. Loman, N.J.; Misra, R.V.; Dallman, T.J.; Constantinidou, C.; Gharbia, S.E.; Wain, J.; Pallen, M.J. Performance comparison of benchtop high-throughput sequencing platforms. Nat. Biotechnol. 2012, 30, 434–439. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Schirmer, M.; D'Amore, R.; Ijaz, U.Z.; Hall, N.; Quince, C. Illumina error profiles: Resolving fine-scale variation in metagenomic sequencing data. BMC Bioinform. 2016, 17, 125. [Google Scholar] [CrossRef] [PubMed]
  66. Minoche, A.E.; Dohm, J.C.; Himmelbauer, H. Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and genome analyzer systems. Genome Biol. 2011, 12, R112. [Google Scholar] [CrossRef]
  67. Bonnaud, E.M.; Troupin, C.; Dacheux, L.; Holmes, E.C.; Monchatre-Leroy, E.; Tanguy, M.; Bouchier, C.; Cliquet, F.; Barrat, J.; Bourhy, H. Comparison of intra- and inter-host genetic diversity in rabies virus during experimental cross-species transmission. PLoS Pathog. 2019, 15, e1007799. [Google Scholar] [CrossRef]
  68. Renner, D.W.; Szpara, M.L. Impacts of Genome-Wide Analyses on Our Understanding of Human Herpesvirus Diversity and Evolution. J. Virol. 2018, 92, e00908-17. [Google Scholar] [CrossRef]
  69. Gong, Y.N.; Tsao, K.C.; Chen, G.W.; Wu, C.J.; Chen, Y.H.; Liu, Y.C.; Yang, S.L.; Huang, Y.C.; Shih, S.R. Population dynamics at neuraminidase position 151 of influenza A (H1N1)pdm09 virus in clinical specimens. J. Gen. Virol. 2019, 100, 752–759. [Google Scholar] [CrossRef]
  70. Höper, D.; Freuling, C.M.; Müller, T.; Hanke, D.; von Messling, V.; Duchow, K.; Beer, M.; Mettenleiter, T.C. High definition viral vaccine strain identity and stability testing using full-genome population data—The next generation of vaccine quality control. Vaccine 2015, 33, 5829–5837. [Google Scholar] [CrossRef] [PubMed]
Viruses 11 00846 g001 550
Figure 1. Workflow for the enrichment of ASFV-specific target sequences using MyBaits® (Arbor Bioscience) prior to Illumina MiSeq sequencing. Extracted DNA including ASFV target sequences (red) is fragmented and library is prepared by adding Illumina specific sequencing adapters (yellow) (A). After initial denaturation (95 °C), the library is cooled to hybridisation temperature (65 °C) and adapter specific blocking oligos (green) are added to prevent the re-hybridisation of adapter sequences during bait hybridisation (B). Biotinylated ASFV-specific RNA baits are added to the blocked library, and the reaction is incubated at hybridisation temperature for 16 h (B). Streptavidin-coated magnetic beads are added and bind to biotinylated RNA baits hybridised to ASFV target sequences (C). After magnetic separation and washing (C), ASFV target sequences are eluted from the baits and amplified (D). ASFV read-enriched libraries are sequenced on an Illumina MiSeq, and the resulting data are analysed in silico (E).
Figure 1. Workflow for the enrichment of ASFV-specific target sequences using MyBaits® (Arbor Bioscience) prior to Illumina MiSeq sequencing. Extracted DNA including ASFV target sequences (red) is fragmented and library is prepared by adding Illumina specific sequencing adapters (yellow) (A). After initial denaturation (95 °C), the library is cooled to hybridisation temperature (65 °C) and adapter specific blocking oligos (green) are added to prevent the re-hybridisation of adapter sequences during bait hybridisation (B). Biotinylated ASFV-specific RNA baits are added to the blocked library, and the reaction is incubated at hybridisation temperature for 16 h (B). Streptavidin-coated magnetic beads are added and bind to biotinylated RNA baits hybridised to ASFV target sequences (C). After magnetic separation and washing (C), ASFV target sequences are eluted from the baits and amplified (D). ASFV read-enriched libraries are sequenced on an Illumina MiSeq, and the resulting data are analysed in silico (E).
Viruses 11 00846 g001
Viruses 11 00846 g002 550
Figure 2. Coverage of ASFV Georgia 2007/1 with ASFV-specific reads from the different sequencing strategies. Blue line, coverage Illumina (left axis); red line, coverage Nanopore (right axis); green line, coverage total (Illumina + Nanopore) (left axis). Indicated below the figure are regions showing low coverage due to inverted terminal repeats (ITR) or G/C homopolymer stretches (red).
Figure 2. Coverage of ASFV Georgia 2007/1 with ASFV-specific reads from the different sequencing strategies. Blue line, coverage Illumina (left axis); red line, coverage Nanopore (right axis); green line, coverage total (Illumina + Nanopore) (left axis). Indicated below the figure are regions showing low coverage due to inverted terminal repeats (ITR) or G/C homopolymer stretches (red).
Viruses 11 00846 g002
Viruses 11 00846 g003 550
Figure 3. Phylogenetic tree showing all available ASFV whole-genome sequences. The maximum-likelihood (ML) tree was constructed using IQ-TREE v1.6.5 based on MAFFT v7.388 aligned ASFV whole-genome sequences. Standard model selection was used, resulting in the best-fit model K3Pu + F + R3 (three substitution types model + empirical base frequencies + FreeRate model with 3 categories). Statistical support of 10,000 ultrafast bootstraps is indicated at the nodes. Taxon names include, where available, ASFV-designation and INSDC accession number.
Figure 3. Phylogenetic tree showing all available ASFV whole-genome sequences. The maximum-likelihood (ML) tree was constructed using IQ-TREE v1.6.5 based on MAFFT v7.388 aligned ASFV whole-genome sequences. Standard model selection was used, resulting in the best-fit model K3Pu + F + R3 (three substitution types model + empirical base frequencies + FreeRate model with 3 categories). Statistical support of 10,000 ultrafast bootstraps is indicated at the nodes. Taxon names include, where available, ASFV-designation and INSDC accession number.
Viruses 11 00846 g003
Viruses 11 00846 g004 550
Figure 4. Comparison of sequencing workflows used in this study. (A) Overview of the possible strategies tested here for sequencing. (1) Illumina shotgun library preparation plus target enrichment prior to Illumina sequencing; (2) pure Illumina shotgun sequencing; (3) Nanopore sequencing. (B) Decision tree derived from the experiences in this study. Depending on the known or expected virus/host ratio and genome complexity, different workflows should be chosen. In this study, ASFV Georgia 2007/1 was treated as being a low virus/host ratio and high genome complexity sample and hence was sequenced according to all strategies (1/2/3), while ASFV Moldova 2017/1 had a low virus/host ratio but a suitable reference sequence was available, and hence target enrichment was sufficient, so only workflow (1) had to be applied.
Figure 4. Comparison of sequencing workflows used in this study. (A) Overview of the possible strategies tested here for sequencing. (1) Illumina shotgun library preparation plus target enrichment prior to Illumina sequencing; (2) pure Illumina shotgun sequencing; (3) Nanopore sequencing. (B) Decision tree derived from the experiences in this study. Depending on the known or expected virus/host ratio and genome complexity, different workflows should be chosen. In this study, ASFV Georgia 2007/1 was treated as being a low virus/host ratio and high genome complexity sample and hence was sequenced according to all strategies (1/2/3), while ASFV Moldova 2017/1 had a low virus/host ratio but a suitable reference sequence was available, and hence target enrichment was sufficient, so only workflow (1) had to be applied.
Viruses 11 00846 g004
Table
Table 1. African swine fever virus (ASFV) whole-genome sequences available from public databases (status 11.03.2019).
Table 1. African swine fever virus (ASFV) whole-genome sequences available from public databases (status 11.03.2019).
NumberAccession NumberASFV IsolateCountry of OriginSubmission DateCollection DateHostP72 GenotypeWGS PublicationMethodCoverage
1NC_001659.2BA71VSpain19951967Vero cellsI[32]Sanger sequencingN/A
2AY261360.1Kenya 1950Kenya20031950Domestic pigXN/AN/AN/A
3AY261362.1Mkuzi 1979South Africa20031979TickXIIN/AN/AN/A
4AY261365WarmbathsSouth Africa2003N/ATickIII / IN/AN/AN/A
5AY261363.1Pretorisuskop/96/4South Africa20031996TickXX / IN/AN/AN/A
6AY261361.1Malawi Lil-20/1Malawi20031983TickVIII[34]N/AN/A
7AY261366.1WarthogNamibia20031980WarthogIVN/AN/AN/A
8AY261364.1Tengani 62Malawi20031962Domestic pigV / IN/AN/AN/A
9AM712239.1Benin 97/1Benin20071997Domestic pigI[35]Sanger sequencingN/A
10AM712240.1OURT 88/3Portugal20071988Domestic pigI[35]Sanger sequencingN/A
11FN557520.1E75Spain20091975Domestic pigI[36] Roche 454 GS FLX, Sanger sequencingN/A
12FR682468.1Georgia 2007/1Georgia20102007Domestic pigII[37] Roche 454 GS FLXN/A
13KM102979.126544/OG10Italy (Sardinia)20142010Domestic pigI[46]Illumina HiScanSQ, Sanger sequencing20
14KJ747406.1Kashino 04/13Russia20142013Wild boarIIN/ASanger sequencingN/A
15KM111295.1Ken06.BusKenya20142006Domestic pigX[39]Illumina HiSeqN/A
16KM262844.1L60Portugal20141960Domestic pigI[38]Amplicon sequencing on Roche 454 GS FLX, Sanger sequencingN/A
17KP055815.1BA71Spain20141971Domestic pigI[44]Sanger sequencingN/A
18KM262845.1NHVSpain20141968Domestic pigI[38]Amplicon sequencing on Roche 454 GS FLX, Sanger sequencingN/A
19KM111294.1Ken05/Tk1Kenya20152005TickIX[39]Illumina HiSeqN/A
20KP843857.1Odintsovo_02/14Russia20152014Wild boarIIN/A Roche 454 GS FLXN/A
21LP643842.1Patent WO2015091322N/A2015N/AN/AN/AN/AN/AN/A
22KX354450.147/Ss/2008Italy (Sardinia)20162008Domestic pigI[40]Illumina MiSeq; PacBioN/A
23MG939585.1Pol16_20540_o10Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
24MG939587.1Pol17_03029_C201Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
25MG939583.1Pol16_20186_o7Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
26MG939588.1Pol17_04461_C210Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
27MG939584.1Pol16_20538_o9Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
28MG939586.1Pol16_29413_o23Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
29MG939589.1Pol17_05838_C220Poland20182016/2017Sus scrofaII[51]Illumina MiSeq20-40
30MH681419.1ASFV/POL/2015/PodlaskiePoland20182015Wild boarII[41]Illumina MiSeq103
31MH766894.1ASFV-SY18China20182018Domestic pigIIN/AN/AN/A
32MH025918.1R25Uganda20182015Domestic pigIX[42]Illumina NextSeq 500869
33MH025920.1R35Uganda20182015Domestic pigIX[42]Illumina NextSeq 5001487
34MH025917.1R7Uganda20182015Domestic pigIX[42]Illumina NextSeq 500439
35MH025916.1R8Uganda20182015Domestic pigIX[42]Illumina NextSeq 500309
36MH025919.1N10Uganda20182015Domestic pigIX[42]Illumina NextSeq 50023
37LS478113.1Estonia 2014Estonia20182014Domestic pigII[43]Illumina MiSeq100
38MH910495.1Georgia 2008/1Georgia20182008Domestic pigII[45]Illumina MiSeq8.5
39MH910496.1Georgia 2008/2Georgia20182008Domestic pigII[45]Illumina MiSeq118
40MK128995.1China/2018/AnhuiXCGQChina20192018Domestic pigII[52]BGISEQ-500271
41LR536725.1Belgium 2018/1Belgium20192018Wild BoarII[53]Illumina MiSeq292
N/A: data not available; WGS: whole-genome sequence.
Table
Table 2. Comparison of shotgun and target enrichment sequencing approaches for the generation of ASFV whole-genome sequences.
Table 2. Comparison of shotgun and target enrichment sequencing approaches for the generation of ASFV whole-genome sequences.
ASFVSample TypeLibrary NumberSequencing ModeTotal ReadsTotal ASFV Reads% ASFV ReadsMean Coverage
Georgia 2007/1Cell culture supernatantlib02645shotgun1,764,0788309 (8150)0.47 (0.46)12.7 (12.5)
lib02645shotgun7,317,74436,268 (33,454)0.5 (0.46)56.5 (52.1)
lib02679myBaits67,17444,862 (40,234)66.78 (59.89)57.2 (51.7)
Moldova 2017/1Spleenlib02487myBaits829,408690,206 (207,763)83.89 (25.0)1055 (317)
lib02577shotgun8,232,5184042 (3986)0.05 (0.048)N/A
() Number/percentage of unique reads and coverage of the corresponding ASFV whole-genome sequence are given in brackets.

Share and Cite

MDPI and ACS Style

Forth, J.H.; Forth, L.F.; King, J.; Groza, O.; Hübner, A.; Olesen, A.S.; Höper, D.; Dixon, L.K.; Netherton, C.L.; Rasmussen, T.B.; et al. A Deep-Sequencing Workflow for the Fast and Efficient Generation of High-Quality African Swine Fever Virus Whole-Genome Sequences. Viruses 2019, 11, 846. https://doi.org/10.3390/v11090846

AMA Style

Forth JH, Forth LF, King J, Groza O, Hübner A, Olesen AS, Höper D, Dixon LK, Netherton CL, Rasmussen TB, et al. A Deep-Sequencing Workflow for the Fast and Efficient Generation of High-Quality African Swine Fever Virus Whole-Genome Sequences. Viruses. 2019; 11(9):846. https://doi.org/10.3390/v11090846

Chicago/Turabian Style

Forth, Jan H., Leonie F. Forth, Jacqueline King, Oxana Groza, Alexandra Hübner, Ann Sofie Olesen, Dirk Höper, Linda K. Dixon, Christopher L. Netherton, Thomas Bruun Rasmussen, and et al. 2019. "A Deep-Sequencing Workflow for the Fast and Efficient Generation of High-Quality African Swine Fever Virus Whole-Genome Sequences" Viruses 11, no. 9: 846. https://doi.org/10.3390/v11090846

APA Style

Forth, J. H., Forth, L. F., King, J., Groza, O., Hübner, A., Olesen, A. S., Höper, D., Dixon, L. K., Netherton, C. L., Rasmussen, T. B., Blome, S., Pohlmann, A., & Beer, M. (2019). A Deep-Sequencing Workflow for the Fast and Efficient Generation of High-Quality African Swine Fever Virus Whole-Genome Sequences. Viruses, 11(9), 846. https://doi.org/10.3390/v11090846

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop