A Method for the Analysis of African Swine Fever by Viral Metagenomic Sequencing

In 2018, there was an outbreak of African swine fever (ASF) in China, which spread to other provinces in the following 3 years and severely damaged China's pig industry. ASF is caused by the African swine fever virus (ASFV). Given that the genome of the African swine fever virus is very complex and whole genome information is currently inadequate, it is important to efficiently obtain virus genome sequences for genomic and epidemiological studies. The prevalent ASFV strains have low genetic variability; therefore, whole genome sequencing analysis provides a basis for the study of ASFV. We provide a method for the efficient sequencing of whole genomes, which requires only a small number of tissues. The database construction method was selected according to the genomic types of ASFV, and the whole ASFV genome was obtained through data filtering, host sequence removal, virus classification, data assembly, virus sequence identification, statistical analysis, gene prediction, and functional analysis. Our proposed method will facilitate ASFV genome sequencing and novel virus discovery.

Introduction

ASF is a highly contagious and fatal swine disease. The pathogen of ASF is the African swine fever virus (ASFV), which is the only member of the family Asfarviridae and the only known DNA arbovirus (1, 2). Depending on the strain, ASFV has a large (170–193 kbp) double-stranded DNA genome containing 151–167 genes, which are involved in viral replication and assembly as well as in modulating host cellular functions and immune evasion (3). The virus can be transmitted through direct contact with infected swine, their products, and soft ticks of the genus Ornithodoros (4). ASF was first described in Kenya in 1921; it then spread to other African, European, Caribbean, and South American countries (3, 5, 6). The disease was introduced into Georgia in 2007 and then spread throughout Eastern Europe, including Russia, Belarus, Ukraine, Estonia, Lithuania, Latvia, Romania, Moldova, Czech Republic, and Poland (7–9). In August 2018, China reported its first ASF outbreak. Within 1 year, ASFV spread rapidly into all provinces in mainland China. The spread of ASF has resulted in huge losses to the Chinese pig industry (3).

ASFV has a high genetic and antigenic diversity. Based on the p72 protein (B646L), 24 genotypes have been identified, while at least eight serotypes are recognized based on hemadsorption inhibition (10, 11). The spread of the African swine fever virus in China is a serious threat to the diversity and survival of pigs. To facilitate much needed epidemiological investigations, advance research, and further vaccine development, it would be expedient to have a simple and reproducible method for full genome sequencing of the ASFV (12). In the early stage, we used the first-generation sequencing technology to sequence the whole genome of ASFV, which is very time-consuming and with heavy workload. This demands for a faster method for rapid sequencing of the whole ASFV genome, and second-generation sequencing is an important tool for sequencing large genomes, which is essential for effective emergency management in the event of disease outbreaks. Current methods of virus enrichment in second-generation sequencing are inefficient and time-consuming. By improving the enrichment method, we can effectively increase the proportion of virus samples and provide more effective data for subsequent analysis. As the ASFV genome contains a wide range of homopolymers and repeat regions, the short-read data generated by second generation sequencing platforms need to be processed carefully. ASFV genome recombination (such as inversion or duplication) may be missed when comparing reference sequences, and the quality of the consistent sequence is heavily dependent on the reference sequence and therefore, may be misassembled. After sequencing, data filtration, host sequence removal, virus classification, data assembly, virus sequence identification, statistical analysis of virus abundance, gene prediction, and functional analysis were performed. Through the analysis of these steps, an accurate sequence is finally obtained. Our results confirm the feasibility of sequencing an ASFV genome directly from positive clinical tissues, and provide a basis for further epidemiological research and evolutionary analysis (13).

Materials and Methods

Experimental Materials

We obtained 0.45 μm and 0.22 μm filter (PVDF) membranes from Merck Display Materials Co., Ltd. (Shanghai, China), OptiPrep Density Gradient Medium purchasedfrom Guangzhou Fan-si Biotechnology Co., Ltd. (Guangzhou, ChinaGuangzhou, ChinaGuangzhou, ChinaGuangzhou, ChinaGuangzhou, China), DNase I (Promega Corporation, Madison, USA), BSA (configured as 1% BSA-SM solution, filtered by 0.22 μm membrane), and gelatin from porcine skin purchased from Sigma-Aldrich Trading Co., Ltd. (Shanghai, China). Furthermore, we used an overspeed centrifuge (Beckman Coulter, USA), HiPure Viral DNA Kit (D3191) purchased from Majorbio (Shanghai, China), and 10 × SM buffer (pH 7.5, 1 M NaCl, 100 mM MgSO4, 500 mM Tris, 0.1% gelatin; working concentration was diluted with ultra-pure water to 1 × SM). Op density gradient solution was made using Optiprep original solution, mixed with 10 × SM buffer (9:1).

Enrichment and Purification of Virus-Like Particles (VLPs)

Cases were identified from the Ministry of Agriculture and Rural Affairs of Huangpu District, Guangzhou (information released by the Ministry of Agriculture and Rural Affairs: http://www.moa.gov.cn/gk/yjgl_1/yqfb/201812/t20181223_6165395.htm). The spleen tissue was completely fixed with 10% formalin solution for 72 h. Two formalin fixed pig spleens 1 g each were taking and cut into small pieces with sterile scissors, and then poured into 2 × sucrose -Triton washing solution for washing, 1,000 r/min centrifugation for 15 min, precipitation was beaten evenly with TE Buffer (pH 9.0), adding SDS to the final concentration of 1%, protease K to 200 μg/mL, and kept in water bath at 48°C for 48 h. And 5–10 mL of pre-cooled 1 × SM buffer (SM buffer filtered by 0.22 μm) was added. The spleen was evenly homogenized and placed in liquid nitrogen and a 37°C water bath alternately three times. Spleens were centrifuged at 4°C for 5 min at 1,000 rpm and 12,000 rpm in succession, and the supernatant was filtered using a 0.45 μm filter. Before using the 0.45 μm filter, 1% BSA solution was filtered through the wetting membrane. Then, qPCR was performed on each liquid layer using specific ASFV primers to determine the virus Cq value in each liquid layer. The liquid layer with the low Cq value was selected, each two liquid layers were mixed, and 1 × μSM buffer solution was added to fully mix. The solution was centrifuged at 160,000 g for 1 h, and the supernatant was discarded. According to the amount of precipitation, 100–500 μL 1 × SM buffer solution was added for resuspension (with repeated pipette mixing to avoid violent shock). DNaseI was added according to the manufacturer's instructions, with treatment at 37°C for 1 h. Single- and double-stranded DNA were excised simultaneously to fragment the DNA for library construction. EDTA (Promega Corporation, Madison, USA) was added according to the manufacturer's instructions to terminate the reaction. The HiPure Viral DNA Kit (D3191) (Majorbio, Shanghai, China) was used to extract viral DNA, and the concentration was measured. Sequencing libraries were generated using NEB Next® Ultra™ DNA Library Prep Kit for Illumina® (New England Biolabs, MA, USA) following manufacturer's recommendations and index codes were added. Finally, the library was sequenced on an Illumina Novaseq 6,000 and 150 bp paired-end reads were generated. Each sample added 507,333,357 reads. (The extracted DNA was sent to China Guangdong Magigene Technology Co., Ltd [Guangdong, China]) for sequencing. All the experiments involving the ASF virus were carried out in a biosafety level (BSL)-3 laboratory at South China Agricultural University (Guangzhou, ChinaGuangzhou, ChinaGuangzhou, ChinaGuangzhou, ChinaGuangzhou, China).

Data Filtering

After obtaining the metagenomic sequencing data of the sample, it was necessary to evaluate the quality of the sequencing data and remove low-quality data to ensure the credibility of the subsequent analysis results. High-quality sequences obtained by quality control were used for downstream data analysis. The quality control process used the software SOAPnuke (14), with the specific processing steps as follows: (1) Removal of Adapter Paired reads; (2) Removal of single-ended reads with N's (N, uncertain base information) and >5% of paired reads; (3) When the single-ended sequencing read was low quality (sQ ≤ 20) and when the number of bases was >20% of the total number of read bases, these paired reads were removed; (4) Replicated reads produced by PCR amplification were removed; (5) Removal of polyX (ATCG) sequences.

Removal of Host Sequences

SOAPaligner (15) and BWA (16) software were used to compare clean reads to the specified host genome and to remove host sequences.

Virus Classification

The comparison software SOAPaligner (15) and BWA (16) were used to compare clean reads to the virus reference database in order to quickly obtain virus classification information in the samples.

Data Assembly

The assembly software IDBA (17), SPAdes (18), metaSPAdes, MEGAHIT (19), and Trinity (20) were used to assemble high-quality reads of each sample to obtain a longer contig sequence. Specific software information is provided in the Results section. Then, the number, length, and N50 statistic of the assembly sequence were counted. BWA software was used to compare high-quality reads to the assembly sequence to calculate the utilization rate of assembled reads, and the assembly effect was evaluated using these statistical data.

Virus Sequence Identification

A variety of methods [including BLAST, HMMSearch (21), and Metagenemark (22)] and databases [including NT, NR, VPFS (23), VFam (24), PFAM (25), and KEGG] were used to identify viral sequences. Annotation was based on the reference database; the corresponding virus sequences were isolated from the NT database, and BLASTN was used to compare contigs with the constructed virus database for species annotation. Using novel virus identification methods to find candidate virus sequences, contigs were compared with multiple databases, as long as one of the following three conditions were satisfied: (1) Comparison between BLASTN (v2.9.0+) and the virus database isolated from NT (virus-NT, including phages) was used to screen the comparison results with e ≤ 1 × 10⁻⁵ (e: exponent); (2) Comparison between BLASTX (v2.9.0+) and the virus-NR database isolated from NR (including bacteriophages) was used to screen the comparison results with e ≤ 1 × 10⁻³; (3) Metagenemark (v3.38) was used to predict the genes, and then HMMSearch (v3.2.1) software was used to compare the protein sequences with the HMM database (VPFS and VFAM), and the comparison results were screened with e ≤ 1 × 10⁻⁵.

Elimination of False Positives

The candidate virus sequences obtained above were compared with the NT database BLAST (v2.9.0+) and screened at e ≤ 1 × 10⁻¹⁰. The sequences not aligned in the previous step were compared with the NR database Diamond (v0.9.10) and screened with e ≤ 1 × 10⁻³. NCBI taxonomy data was used to annotate the above-mentioned alignment results. If more than 20% of the alignment results in the first 50 alignment results were non-viral sequences (annotated results were Eukaryota, Bacteria, and Archaea), the sequences were considered to be non-viral sequences, and the rest were considered to be viral sequences. Virus contigs were annotated according to the best hit comparison results of virus contigs and virus-NT (e ≤ 1 × 10⁻⁵).

Virus Abundance Statistics

Reads were compared with identified virus contigs, and the reads per kilo bases per million reads (RPKM) values of each contig were calculated for comparative analysis between samples.

$\begin{array}{l}\text{RPKM}=\frac{\text{Contig reads}}{\text{Total mapped reads}(\text{millions})*\text{ Contig length}(\text{KB})}\end{array}$

Note: (1) Contig reads: number of reads in a contig; (2) Total mapped reads: number of reads in millions; (3) Contig length: contig length, in kbp.

Gene Prediction

Metagenemark (22) was used to predict the gene sequences of virus contigs, and the number and length of the predicted genes were evaluated.

Functional Analysis

The predicted gene protein sequence and UniProtKB/Swiss-Prot database sequence of the virus [ViralZone (26), reviewed proteins, https://viralzone.expasy.org/] were used for functional annotation information.

ASFV Statistical Analysis

The contig on ASFV was identified and compared according to the final virus sequence for analysis. MAFFT (27) software was used to perform multiple comparisons between the ASFV genome and the ASFV reference strain genome.

Results

DNA Extraction

After centrifugation, the products were layered into Eppendorf tubes, which were weighed before use and after adding each liquid layer. The weights of the liquid layers were recorded and the density of each liquid layer was calculated (Figure 1A) to assess the stratification. The DNA of each liquid layer was extracted with the OMEGA nucleic acid extraction kit, and then qPCR was performed on each liquid layer with specific African swine fever virus primers to determine the virus Cq value in each liquid layer (Figure 1B). The results showed that the Cq value of layer 16–19 was the lowest, and the virions were enriched in layer 16–19. Then, the two liquid layers adjacent to layer 16–19 were mixed, 1 × SM buffer was added, and the virions were fully mixed.

FIGURE 1

Figure 1. Density and Cq value stratified graphs. (A) Density diagram of each liquid layer after overspeed centrifugation; (B) Cq value of each liquid layer detected by qPCR after ultracentrifugation. (A,B) Are the same two samples of ultracentrifugation.

Quality Control of Sequencing Data

To ensure the accuracy of subsequent analysis, SOAPnuke (v2.0.5) software was first used to process the raw data from the machine, and high-quality clean reads were obtained. The quality distribution is shown in Figure 2. The results showed that the average base mass (green line) at all positions was above 30, and the quality of the data was high enough to be used in the following analysis.

FIGURE 2

Figure 2. Base quality values of sequencing data from the Illumina sequencing platforms. Quality values are represented by Qphred, where Qphred = −10 × log₁₀(e), with e representing the base sequencing error rate. The corresponding relationship between base recognition error rate and Phred score in Illumina Base Calling software is as follows: when Phred score is 10, the base correct recognition rate is 90%; when the Q-score is Q10 and Phred score is 20, the correct base identification rate is 99%; when the Q-score is Q20 and Phred score is 30, the correct base recognition rate is 99.9%; when the Q-score is Q30 and Phred score is 40, the correct base recognition rate is 99.99%; and when the Q-score is Q40, the correct base identification rate is 99.99%. The abscissa represents the position of the base, the ordinate represents the quality of the base at each position, the left graph represents the data before quality control, and the right graph represents after quality control.

Removal of Host Contamination

To avoid the influence of the host sequence on subsequent analysis, BWA (v0.7.17, default parameter: mem-k30) software was used to compare clean reads with the host database. Sus scrofa was used as the host reference information (accession: NC_010443.5). The comparison results where the length of the comparison was >80% of the total read length were filtered, and then the corresponding sequence was removed. The results showed that 16.5% of clean reads (PE: paired-end) were obtained after host removal, as shown in Table 1.

TABLE 1

Table 1. Removal of host sequence and virus sequence statistics.

Virus Composition Analysis

To quickly obtain virus composition information in samples, BWA (v0.7.17, default parameter: mem-k30) software was used to compare clean reads to virus reference data (isolated from NT data); the comparison results were filtered if the length of the comparison was >80% of the total read length. Statistical analysis showed that virus reads (PE) accounted for 3.19% of clean reads (PE), as shown in Table 1. According to the annotation information of the NCBI Taxonomy Database, the virus classification information was counted. To improve the accuracy of the results, comparison results with >5 reads covered were filtered during species annotation. Reads (PE) of the African swine fever virus family (Asfarviridae, red column in the figure) accounted for 87.12% of the total reads (PE). Asfarviridae had only one ASFV member, so it was not necessary to make annotations at the genus level for the display. The statistical results are shown in Figure 3.

FIGURE 3

Figure 3. Statistics of the results of horizontal annotation of the viral family. Note: (1) Species with proportion of reads >1% were selected for display; proportion of virus = number of single virus reads/number of total virus reads × 100.

Data Assembly

MEGAHIT (v1.1.2, default parameters: —presets meta-large —min-contig-len 300) software was used to assemble clean data, and BWA (v0.7.17) software was used to compare clean reads with assembly results. The calculated utilization rate of the reads was 75.30%. At the same time, BLAST (v2.9.0+) software was used to align the assembled contigs with the host sequence and to remove the host sequence. The statistical results showed that the maximum number of 400 bp contigs was 3,896 (Figure 4), and the proportion of host sequences was 0.57% (Table 2).

FIGURE 4

Figure 4. Contig length distribution. (1) The abscissa represents the length of contigs, and the ordinate represents the number of contigs; (2) The red number represents the number of contigs in the length segment.

TABLE 2

Table 2. Removal of host sequence statistics.

Virus Sequence Identification

Identification Based on Reference Sequence

The obtained contigs were compared to the virus database (isolated from the NT database) using BLAST (v2.9.0+) software. Comparison and screening criteria were as follows. (1) Confirmative: comparison of similarity ≥80%, compare length ≥500 bp, e ≤ 1 × 10⁻⁵, and high reliability; (2) Suspected: did not meet the conditions in (1), compare length ≥100 bp, e ≤ 1 × 10⁻⁵, and low credibility, requiring in-depth analysis and verification; (3) According to NCBI Taxonomy annotation information, the number of phages and other viruses, as well as the types of RNA and DNA, were counted (Tables 3–5, and Figure 5).

TABLE 3

Table 3. Virus sequence statistics.

TABLE 4

Table 4. Virus category statistics.

TABLE 5

Table 5. Virus category statistics.

FIGURE 5

Figure 5. Virus annotation statistics (family level). (1) abscissa: number of contigs; ordinate: comment information.

Identification of De novo Virus Based on Sequence

Based on the comparison of reference data, there were two defects: (1) Only known viruses could be identified; (2) There were false positives in the results. To reduce false positives and identify unknown viruses, new methods were used to further identify virus sequences by combining multiple databases (Table 6). The analysis process is described in Section 2.7.

TABLE 6

Table 6. De novo virus identification result statistics.

Comparison of the Two Methods

The virus sequences obtained by the two methods were compared, and Venn diagram were drawn. Then, the confirmed and novel virus sequences (Section 3.6.1 and 3.6.2) were combined as the final virus sequences (Figure 6 and Table 7).

FIGURE 6

Figure 6. Venn and annotated statistical diagrams of the virus sequence. (A) Comparison of the number of contigs obtained by different methods. (1) Novel method to obtain virus sequence; (2) Confirmed: confirmed virus sequence; (3) Suspected sequence of a virus. (B) Statistics of final viral sequence annotations (family level).

TABLE 7

Table 7. Virus contig statistics.

Virus Abundance Analysis

BWA (v0.7.17, default parameter: mem-k 30) software was used to compare the host-removed clean reads with the obtained virus contigs, and the proportion of virus reads was calculated when the filter comparison length was lower than 80% of the total read length. The RPKM value of each viral contig was calculated. The maximum RPKM value was 85,317.57 and the minimum RPKM value was 224.62 (Table 8).

TABLE 8

Table 8. Comparison results of clean reads and virus contigs.

Gene Prediction

Metagenemark (v3.38) software was used to predict the gene corresponding to virus contigs, and sequences with gene nucleic acid length >150 bp were filtered. The statistical results are presented in Table 9.

TABLE 9

Table 9. Gene prediction result statistics.

Functional Analysis

Based on the BLASTP (v2.9.0+) software, the gene sequences were compared with those collected from the UniProtKB/SWISS-PROT database (https://viralzone.expasy.org/). The best hit screened with e < 1 × 10⁻³ was compared to obtain the function information of the virus. The results showed that the percentage of genes in the database accounted for 37.14% of the total gene number (Table 10). Functional genes annotated by UniProtKB category statistics included biological processes, cell components, and molecular functions (Figure 7). The polymerase gene number statistics are shown in Table 11.

TABLE 10

Table 10. Gene annotation rate statistics.

FIGURE 7

Figure 7. UniProtKB category statistics.

TABLE 11

Table 11. Polymerase gene number statistics.

ASFV Analysis

Clean reads were compared to the reference sequence MT496893.1 by using BWA and SamTools software, and a consensus sequence was obtained. The average sequencing depth was 1,233. According to the final virus sequence in Section 3.6.3, the ASFV contig was found and compared with the reference sequence MT496893.1. Sequences were filtered using similarities <90% and lengths <500 bp; the results showed that there were five matched contigs, and the longest was 114,708 bp. After overlapping, correcting, and assembling with the MT496893.1 sequence, the whole genome sequence, named GZ201801-1, was obtained with a length of 188,035 bp. However, the assembled genome GZ201801-1 was 1,358 bp shorter than GZ201801. To construct the phylogenetic tree, representative ASFV sequences from GenBank were selected, including genotype II ASFV strain (n = 12), genotype I ASFV strain (n = 8), genotype II ASFV strain (n = 6), genotype I ASFV strain (n = 2), and genotype VIII, IV, XX, III, V, and XII ASFV strains (each one). Phylogenetic analysis of the entire viral genome was performed using the NJ method in Mega X, and the phylogenetic tree bootstrap value was set to 1,000 (Table 12 and Figure 8).

TABLE 12

Table 12. Statistical analysis of total length and similarity of ASFV.

FIGURE 8

Figure 8. Sequencing depth map and genome-wide phylogenetic tree. (A) Sequencing depth map of African swine fever virus. (B) Genome-wide phylogenetic tree of the GZ201801 strain.

Discussion

The ASFV genome is relatively conserved and its natural variation is very slow, but interaction with the host and stimulation by environmental factors can accelerate the transmission of ASFV (28). Given the current lack of genome-wide information, it is necessary to promote the whole-genome sequencing of ASFV. In recent years, the development of sequencing technology has provided technical support for whole genome sequencing, including the macro virus group method (29).

The macro virus group (virome) method directly uses the genetic material of all viruses in the samples for analysis. First, virus particles are enriched, the genome sequence information is obtained, and then the composition and relative abundance of all viruses are identified. This method is a powerful means for the discovery of new viruses and the early detection and control of virus infection. It is of great use in the study of the origin and evolutionary patterns of viruses, genetic diversity, geographical distribution, and the relationship between viruses and their hosts. In the early stage, we used first-generation sequencing to sequence the whole genome of the Guangdong outbreak ASFV, and submitted the sequence to GenBank (serial number: MT496893.1). Due to the time-consuming characteristic of first-generation sequencing, we improved the sequencing analysis method based on the macro virus group (virome) method. Previous studies used a density gradient buffer of Percoll, sucrose, and cesium chloride (30), whereas our steps of using Optiprep separation and improved centrifugation effectively increased the proportion of virus samples provided more effective data for subsequent analysis. Optiprep Separation liquid is characterized by high density, low viscosity, and low permeability. Continuous or discontinuous isotonic gradient solutions can be formed when appropriate concentrations of buffers or basic media are added to separate various cells, nuclei, organelles, and lipoproteins. Compared with the traditional separation solutions (such as sucrose, Percoll, and cesium chloride), it not only has good separation and purification effect, but also has no effect on the life activities of various cells and organelles. In this experiment, the virions in solution were purified and concentrated by overspeed centrifugation of spleens, and the enrichment layer of the virions was then determined by qPCR to obtain higher purity virions for viral DNA extraction. High-purity samples yielded high-quality data for downstream assembly and analysis. Data quality control, host sequence removal, virus classification, data assembly, virus sequence identification, statistical analysis of virus abundance, gene prediction, and functional analysis were performed. Due to the high error rate of second-generation sequencing methods, we used different virus identification methods to restore the diversity and abundance of viruses in the sample as much as possible.

The results of data quality control showed that the average base mass at all positions was above 30, indicating high data quality. A total of 1,351,282 reads (3.88% of the 34,783,974 reads) were mapped to the entire genome region of the MT496893.1 strain, with an average depth of 1,233. Finally, the MT496893.1 strain was used as a reference for genome correction and assembly, and a 188,035 bp genome was finally obtained. If tissue can be ground and supercentrifuged, the high purity and concentration of virus particles can be theoretically stratified for second-generation sequencing, and virus isolation is not required, which will shorten the time to obtain the results (31). Using clinical samples and overspeed centrifugation to obtain high purity of virus particles, we confirmed that we could get high quality sequencing data; this was advantageous to the subsequent assembly and analysis, due to more timely whole genome sequencing of the outbreak strain, facilitating the analysis of the origin and evolution of ASFV, genetic diversity, and geographical distribution. In summary, our results confirm the utility of virome sequencing as a non-culture direct sequencing method for ASFV genomes from PCR-positive clinical tissues.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Author Contributions

HW and GZ: conceived and designed the experiments. CJ, JJ, and YW: performed the experiments. YCh, YW, ZM, MC, HY, GLi, and LG: sample collection. CJ, CQ, YCh, and GLi: analyzed the data. CJ and HW: contributed to the writing. All authors have read and approved the final manuscript.

Funding

The Key-Area Research and Development Program of Guangdong Province (2019B020211003 and 2019B020211005). National Natural Science Foundation of China (31941004) and China Agriculture Research System of MOF and MARA.

Conflict of Interest

CJ, YW, YCh, GZ, and HW were employed by company Shenzhen Kingkey Smart Agriculture Times Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank Editage (www.editage.cn) for English language editing.

References

1. Wang A, Jia R, Liu Y, Zhou J, Qi Y, Chen Y, et al. Development of a novel quantitative real-time PCR assay with lyophilized powder reagent to detect African swine fever virus in blood samples of domestic pigs in China. Transboundary and Emerg Dis. (2020) 67:284–97. doi: 10.1111/tbed.13350

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Luo Y, Atim SA, Shao L, Ayebazibwe C, Sun Y, Liu Y, et al. Development of an updated PCR assay for detection of African swine fever virus. Arch Virol. (2017) 162:191–99. doi: 10.1007/s00705-016-3069-3

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zhao D, Liu R, Zhang X, Li F, Wang J, Zhang J, et al. Replication and virulence in pigs of the first African swine fever virus isolated in China. Emerg Microbes Infect. (2019) 8:438–47. doi: 10.1080/22221751.2019.1590128

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Alejo A, Matamoros T, Guerra M, Andrés G. A proteomic atlas of the African swine fever virus particle. J Virol. (2018) 92:e01293–18. doi: 10.1128/JVI.01293-18

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Aira C, Ruiz T, Dixon L, Blome S, Rueda P, Sastre P. Bead-based multiplex assay for the simultaneous detection of antibodies to African swine fever virus and classical swine fever virus. Front Vet Sci. (2019) 6:306. doi: 10.3389/fvets.2019.00306

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Barongo MB, Bishop RP, Fèvre EM, Knobel DL, Ssematimba A. A mathematical model that simulates control options for African swine fever virus (ASFV). PLoS ONE. (2016) 11:e0158658. doi: 10.1371/journal.pone.0158658

PubMed Abstract | CrossRef Full Text | Google Scholar

7. 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. doi: 10.1016/j.virusres.2012.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Lu G, Pan J, Zhang G. African swine fever virus in Asia: its rapid spread and potential threat to unaffected countries. J Infec. (2020) 80:350–71. doi: 10.1016/j.jinf.2019.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Revilla Y, Pérez-Núñez D, Richt JA. African swine fever virus biology and vaccine approaches. Adv Virus Res. (2018) 100:41–74. doi: 10.1016/bs.aivir.2017.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Wang T, Sun Y, Huang S, Qiu HJ. Multifaceted immune responses to African swine fever virus: Implications for vaccine development. Vet Microbiol. (2020) 249:108832. doi: 10.1016/j.vetmic.2020.108832

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ge S, Li J, Fan X, Liu F, Li L, Wang Q, et al. Molecular characterization of African swine fever virus, China, 2018. Emerg Infect Dis. (2018) 24:2131–3. doi: 10.3201/eid2411.181274

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Olasz F, Mészáros I, Marton S, Kaján GL, Tamás V, Locsmándi G, et al. A simple method for sample preparation to facilitate efficient whole-genome sequencing of African swine fever virus. Viruses. (2019) 11:1129. doi: 10.3390/v11121129

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Jia L, Jiang M, Wu K, Hu J, Wang Y, Quan W, et al. Nanopore sequencing of African swine fever virus. Sci China Life Sci. (2020) 63:160–4. doi: 10.1007/s11427-019-9828-1

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chen Y, Chen Y, Shi C, Huang Z, Zhang Y, Li S, et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience. (2018) 7:1–6. doi: 10.1093/gigascience/gix120

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Li R, Li Y, Kristiansen K, Wang J. SOAP: short oligonucleotide alignment program. Bioinformatics. (2008) 24:713–14. doi: 10.1093/bioinformatics/btn025

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. (2009) 25:1754–60. doi: 10.1093/bioinformatics/btp324

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Peng Y, Leung HC, Yiu SM, Chin FY. IDBA-UD a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. (2012) 28:1420–28. doi: 10.1093/bioinformatics/bts174

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. (2012) 19:455–77. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Li D, Luo R, Liu CM, Leung CM, Ting HF, Sadakane K, et al. MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. (2016) 102:3–11: 10.1016/j.ymeth.2016.02.020

Google Scholar

20. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. (2011) 29:644–52. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Eddy SR. A new generation of homology search tools based on probabilistic inference. Genome Inform. (2009) 23:205–11. doi: 10.1142/9781848165632_0019

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acid Res. (2010) 38:e132. doi: 10.1093/nar/gkq275

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Paez-Espino D, Eloe-Fadrosh EA, Pavlopoulos GA, Thomas AD, Huntemann M, Mikhailova N, et al. Uncovering earth's virome. Nature. (2016) 536:425–30. doi: 10.1038/nature19094

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Skewes-Cox P, Sharpton TJ, Pollard KS, DeRisi JL. Profile hidden Markov models for the detection of viruses within metagenomic sequence data. PLoS ONE. (2014) 9:e105067. doi: 10.1371/journal.pone.0105067

PubMed Abstract | CrossRef Full Text | Google Scholar

25. El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The Pfam protein families database in 2019. Nucleic Acid Res. (2019) 47:427–32. doi: 10.1093/nar/gky995

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hulo C, de Castro E, Masson P, Bougueleret L, Bairoch A, Xenarios I, et al. ViralZone: a knowledge resource to understand virus diversity. Nucleic Acid Res. (2011) 39:D576–582. doi: 10.1093/nar/gkq901

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Nakamura T, Yamada KD, Tomii K, Katoh K. Parallelization of MAFFT for large-scale multiple sequence alignments. Bioinformatics. (2018) 34:2490–2. doi: 10.1093/bioinformatics/bty121

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Dixon LK, Chapman DA, Netherton CL, Upton C. African swine fever virus replication and genomics. Virus Res. (2013) 173: 3–14. doi: 10.1016/j.virusres.2012.10.020

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Thurber RV, Haynes M, Breitbart M, Wegley L, Rohwer F. Laboratory procedures to generate viral metagenomes. Nat Protoc. (2009) 4: 470–83. doi: 10.1038/nprot.2009.10

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Carrascosa AL, Del VM, Santaren JF, Vinuela E. Purification and properties of African swine fever virus. J Virol. (1985) 54:337–44. doi: 10.1128/jvi.54.2.337-344.1985

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Forth JH, Forth LF, King J, Groza O, Hubner A, Olesen AS, 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. doi: 10.3390/v11090846

PubMed Abstract | CrossRef Full Text | Google Scholar