Development and Application of a Duplex Droplet Digital Polymerase Chain Reaction Assay for Detection and Differentiation of EP402R-Deleted and Wild-Type African Swine Fever Virus.

African swine fever (ASF) is a highly fatal porcine disease caused by the African swine fever virus (ASFV), and resulting in huge economic losses across the globe. ASF has been raging in China for 3 years, and recently EP402R-deleted ASFV strains emerged, showing sub-acute or chronic symptoms in pigs and providing novel difficulties to monitor and control the disease as EP402R-deleted strains possess no hemadsorption (HAD) ability. In addition, the gene deletion virus with low viral load is prone to results retest or false negative due to the high cycle threshold (Ct) value under the current real-time polymerase chain reaction (PCR) detection method. Thus, a new method is needed to detect and distinguish wild strains and gene-deleted viruses. In this study, a duplex droplet digital polymerase chain reaction (ddPCR) assay based on the ASFV B646L and EP402R genes was established and showed good linearity (R2 > 0.99). The limit of detection for duplex ddPCR was 52 copies per reaction and 8.6 copies per reaction for B646L and EP402R, respectively. No cross-reaction with other porcine viruses [classical swine fever virus (CSFV), porcine reproductive and respiratory syndrome virus (PRRSV), porcine epidemic diarrhea virus (PEDV), porcine parvovirus (PPV), Japanese encephalitis virus (JEV), and porcine circovirus type 2 (PCV2)] was identified by this assay. In addition, 44 ASFV-suspicious clinical samples as well as EP402R-deleted ASFV were tested in parallel by duplex real-time PCR and ddPCR, indicative of a higher sensitivity which belonged to the duplex ddPCR assay. In summary, this is the first time that duplex ddPCR assay has been successfully developed to provide an efficient method to detect and differentiate ASFV wild-type and gene-deleted strains.

Introduction

African swine fever (ASF) is a highly fatal porcine disease caused by the African swine fever virus (ASFV), leading to huge economic losses and food security threats worldwide. The ASFV genome is between 171 and 193 kbp, encoding more than 150 proteins and forming a mature ASFV virion with a multilayered structure and overall icosahedral morphology (1). The first ASF outbreak in China was reported in August 2018, and sequencing result revealed a genotype II ASFV (2, 3). Currently, ASF has raged in Asia and spread over 16 countries including Korea (4), Indonesia (5), and Vietnam (6). Pigs intramuscularly inoculated with epidemic strain Pig/HLJ/18 demonstrated acute clinical symptoms including fever, loss of appetite, depression, lethargy, respiratory distress, vomiting, and death in all of the pigs (3). However, EP402R gene (encoding CD2v protein)-deletion strains were isolated from China in 2020, which caused sub-acute or chronic disease, and persistent infection in pigs (7). It is worth mentioning that EP402R gene-deletion strains lost hemadsorption (HAD) ability, and hence conventional HAD assay could not detect EP402R gene-deletion virus infection. EP402R-deleted vaccine candidates, on the other side, also caused sub-acute or chronic symptoms in pigs (8, 9). No ASFV vaccine candidate has been licensed, but once the gene-deleted vaccine is successfully developed and commercialized in the future, it would be widely used for the control of ASF worldwide. Therefore, a novel method is needed to distinguish EP402R-deleted ASFV strains or vaccines developed in the future.

On account of the complex structure of the ASFV virion and the limited knowledge on ASFV infection mechanisms and immunity, no commercial vaccines and therapeutic drugs have been developed yet. Therefore, early phase monitoring of diseases, culling of infected animals, and strict biosecurity measures are dominating methods of disease prevention and control at this stage. A mass of detection methods for ASFV, such as virus isolation, HAD assay, fluorescent antibody detection, polymerase chain reaction (PCR), and enzyme-linked immunosorbent assay (ELISA) has been established and applied (10–14). Among them, real-time polymerase chain reaction (real-time PCR) is widely applied for early detection of pathogens in clinical samples due to rapid, highly sensitive, convenient and specific features. However, real-time PCR methods rely on external standards, standard curves, and the cycle within the amplification process that the reporter dye signal exceeds a threshold [cycle threshold (Ct) value].

Droplet digital polymerase chain reaction (ddPCR) is a novel pathogen detection technology that enables absolute quantification of nucleic acids without using the standard samples. In a ddPCR assay, the amplification reaction is divided into thousands of independent partitions containing either a single target molecule or none at all in each individual reaction mixture. After amplification in partitioned reactions to the endpoint with positive and negative partitions is counted, the number of DNA targets per partition and the copy number in the original sample are calculated by Poisson's Law. On comparing with real-time PCR, one of the golden standard methods for ASFV detection, ddPCR has advantages of absolute quantification, less affected by sample inhibitors, and better detection of low-copy-number samples (15). Methods of ASFV detection by ddPCR with a high degree of linearity and specificity have been established (16, 17). However, current ASFV ddPCR strategies only detect single viral gene and possess no ability to distinguish gene-deleted virus from wild strains.

In this study, based on the B646L gene (conserved and essential gene encoding the major capsid protein p72) and the EP402R gene (encoding the viral hemagglutinin protein CD2v), a probe-based duplex ddPCR was developed to distinguish wild-type or gene-deleted strains ASFV infection. The duplex ddPCR assay in this study showed superior characteristics, such as highly sensitivity, low detection limit, and gene-deleted ASFV strains differentiation capability, and thus possessed the potential to monitor both wild-type ASFV and EP402R-deleted strains infection.

Materials and Methods

Recombinant Standard Plasmid and Virus Genome Preparation

This sentence has been rewritten as “The full length of B646L gene (1,941 base pairs) and EP402R gene (1,083 base pairs) were obtained from ASFV strain GZ201801 (GenBank accession number MT496893.1) by PCR and cloned into pMD18-T vector, which named pB646L and pEP402R. Primers were showed in Table 1. GZ201801-ASFV genome was extracted from inactivated cell-cultured ASFV supernatant by the AxyPrep Body Fluid Viral DNA/RNA Miniprep Kit (Axygen; Corning Inc., Corning, NY, United States). Genomes of classical swine fever virus (CSFV), porcine reproductive and respiratory syndrome virus (PRRSV), porcine epidemic diarrhea virus (PEDV), porcine parvovirus (PPV), Japanese encephalitis virus (JEV), porcine circovirus type 2 (PCV2), ΔEP402R/ΔMGF360-505-ASFV, ΔEP402R-ASFV, and ΔMGF360-505-ASFV are preserved in the South China Agricultural University (18).

Table 1

Primers used for the recombinant standard plasmids construction.

PrimerQ11	Sequence(5′-3′)	Position
pB646L-F	TTAGGTACTGTAACGCA	103607–103623
pB646L-R	ATGGCATCAGGAGGA	105547–105533
pEP402R-F	ATGATAATACTTATTTTTTTAATAT	73383–73407
pEP402R-R	TTAAATAATTCTATCTACGT	74446–74465

Primer and Probe Design

To design the generalized primers and probe of EP402R genes, 37 ASFV genomes were obtained from GenBank and used for multiple sequence alignments by DNAMAN. The upstream primer sequences (5′-ATGTTGAAGAAATAGAAAGTC-3′) were fully matched to the target gene of all viral genomes, except for individual gene mutations in the downstream primer (5′-GACTGTAAGGCTTAGGAA-3′) and probe [5′-(HEX)-TGACACCACTTCCATACATGAACCA-(BHQ1)-3′] (Table 2). The TaqMan real-time PCR detection of the ASFV B646L gene was carried out according to the Office International des Epizooties (OIE)-recommended procedure described previously (19). Moreover, the designed primers and probes were fully matched to genotype II (dominating genotype in China) and genotype I (emerging genotype in China) ASFV. In duplex detection, the probe for the B646L gene was labeled with the 5′-reported dye FAM and the 3′-quencher MGB, while the probe for the EP402R gene was labeled with the 5′-reported dye HEX and the 3′-quencher BHQ1. All primers and probes were purchased from Sangon Biotech (Shanghai, China).

Table 2

Primer and probe of EP402R were aligned with 37 African swine fever virus (ASFV) strains.

ASFV isolate	Probe	Anti-sense primer
M78 (genotype V)	TGACACCACTTCTATACATGAACCA	GACTATATGGCTTAGGAA
MK-200 (genotype V)	TGACACCACTTCTATACATGAACCA	GACTATATGGCTTAGGAA
Ndjassi-77 (genotype I)	TGACACCACTTCCATACATGAACCA	GACTATAAGGCTTAGGAA
Spencer (genotype I)	TGACACCACTTCCATACATGAACCA	GACTATAAGGCTTAGGAA
L-50 (serotype I)	TGACACCACTTCCATACATGAACCA	GACTATAAGGCTTAGGAA
TS7 (genotype X)	TGACACCACTTCTATACATGAACCG	GACTGTAAGGCTTAGGAA
691/88 (serotype IV)	TGACACCACTTCTATACATGAACCG	GACTGTAAGGCTTAGGAA
Katanga (genotype I)	TGACACCACTTCCATACATGAACCA	GACTATAAGGCTTAGGAA
Krasnodar_2012/dom (serotype VIII)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Tver_2012/wb (serotype VIII)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Uganda (genotype X)	TGACACCACTTCTATACATGAACCG	GACTGTAAGGCTTAGGAA
Volgograd_2012/dom (serotype VIII)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
TSP80 (genotype X)	TGACACTACTTCTATACATGAGCCG	GACTGTAAGGCTTAGGAA
Estonia 2014 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
CN201801 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Wuhan 2019-1 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Wuhan 2019-2 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
ASFV-wbBS01 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
CADC_HN09 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
CN/2019/InnerMongolia-AES01 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
DB/LN/2018 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Pig/HLJ/2018 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
BA71V (genotype I)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
ASFV/Timor-Leste/2019/1 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
ASFV/pig/China/CAS19-01/2019 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
ASFV/Ulyanovsk 19/WB-5699 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
China/2018/AnhuiXCGQ (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
CN/2019/InnerMongolia-AES01 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Georgia_2007-1 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Georgia_2008-1 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
GZ201801 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
HLJ-2018 (genotype II)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
Ken05-Tkl (genotype X)	TGACACCACTTCTATACATGAACCA	GACTGTAAGGCTTAGGAA
Ken06-Bus (genotype X)	TGACACTACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
L60 (genotype I)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
NHV (genotype I)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA
OURT88/3 (genotype I)	TGACACCACTTCCATACATGAACCA	GACTGTAAGGCTTAGGAA

The bold letters indicate the mutated bases.

Duplex ddPCR and Real-Time PCR Amplification Conditions

For duplex ddPCR assays, 20 μl of solution contained 10 μl of ddPCR supermix (Bio-Rad, Hercules, CA, USA), 2 μl of template, 900 nM primer and 500 nM probe at final concentration. The mixtures were then emulsified automatically with 70 μl droplet generator oil (QX-200 droplet generator; Bio-Rad), and 40 μl mixed droplets were transferred to a 96-well-reaction plate to heat seal at 180°C for 5s. The ddPCR cycling conditions were as follows: 95°C predenaturation for 10 min, followed by 40 cycles of denaturation at 94°C for 30s and annealing at 53.9°C for 1 min, and 1 cycle of enzyme inactivation at 98°C for 10 min. Finally, the plate was placed into a droplet reader (QX-200 droplet reader; Bio-Rad). For duplex real-time PCR assays, a 20 μl reaction mixture was made containing 10 μl of 2×AceQ real-time PCR Probe Master Mix (Vazyme Nanjing, Jiangsu, China), 0.4 μl of each primer (final concentration of 400 nM), 0.2 μl of each probe (final concentration of 200 nM), and 2μl of DNA template. The duplex real-time PCR was performed on the LightCycler96 Real-Time PCR System with cycling conditions as follows: 95°C for 5 min, followed by 40 cycles of 95°C for 10 sand 60°C for 30 s.

Sensitivity and Specificity of ddPCR and Real-Time PCR

To evaluate the sensitivity of the primers and probes of the B646L gene and EP402R gene, gradient-diluted pB646L and pEP402R plasmid were used as templates and tested by ddPCR and real-time PCR. Briefly, copy numbers of the standard plasmids of pB646L and pEP402R were calculated based on the following formula:

These two methods were then performed in parallel using the same standard plasmid as the template after 10-fold serial dilutions. To evaluate the specificity of the primers and probes of the B646L gene and EP402R gene, genomes of ASFV and other six common swine viruses (CSFV, PRRSV, PEDV, JEV, PPV, and PCV2) were used as templates and tested by duplex ddPCR with ASFV-specific primers and probes.

Clinical Sample Detection

To compare the diagnostic specificity and identify the ability to distinguish between wild strains and gene-deletion viruses of duplex ddPCR and real-time PCR assay, a total of 44 pig blood samples collected from suspected ASFV-infected pigs were processed for DNA extraction and evaluated by duplex ddPCR and real-time PCR. EP402R-positive samples (inactivated GZ201801 and ΔMGF360-505-ASFV) and EP402R-negative samples (water, ΔEP402R-ASFV, and ΔEP402R/ΔMGF360-505-ASFV) were set at the same time.

Results

Optimization of the Duplex ddPCR Assay

To validate the optimum annealing temperature of the duplex ddPCR assay, pB646L and pEP402R were used as templates and a gradient annealing temperature was set at 60, 59.4, 58.3, 56.3, 53.9, 52, and 50.7°C. Blue droplets and green droplets were shown as positive droplets for B646L and EP402R, respectively, while black droplets were negative droplets (Figure 1). For the primers and probe of the B646L gene, positive droplets were detected under all annealing temperature conditions but were differentiated well above 53.9°C. As for the primers and probe of the EP402R gene, positive droplets were differentiated well below 53.9°C. Therefore, an annealing temperature of 53.9°C was determined for the duplex ddPCR assay.

Figure 1

Annealing temperature optimization of duplex droplet digital polymerase chain reaction (ddPCR) assay. The duplex ddPCR assay for pB646L (A) and pEP402R (B) was optimized with different annealing temperatures and the fluorescence amplitude of amplification determined the optimal reaction condition. The annealing temperature gradient (lanes 1–8) was 60, 59.4, 58.3, 56.3, 53.9, 52, 50.7, and 50°C, respectively.

Linearity and Limit Detection of Duplex ddPCR and Real-Time PCR Assay

To evaluate the linearity, sensitivity, and efficiency of the quantification of the ASFV duplex real-time PCR and ddPCR, 10-fold dilution standard plasmids of pB646L and pEP402R with concentrations ranging from 106 to 100 copies/μl were used as templates and tested by duplex real-time PCR and ddPCR. The standard curve of the duplex real-time PCR was plotted using the LightCycler96 Real-Time PCR System (ROCHE). As shown in Figures 2A,B, both duplex real-time PCR and ddPCR have a strong linear correlation (R2 > 0.99) between Ct values and the corresponding copy numbers of pB646L and pEP402R. The standard curves of pB646L and pEP402R were plotted with slopes of−3.9345 and−3.2950, respectively. The lowest detectable concentration of pB646L and pEP402R was 10 copies/μl in duplex real-time PCR assay but 1 copy/μl in duplex ddPCR assay (Figures 2C,D). In addition, the actual concentration value was revised as 52 copies/reaction and 8.6 copies/reaction by Poisson distribution for B646L and EP402R, respectively. These results indicated that the primers and probes used in the ddPCR assay possessed higher linearity, sensitivity, and efficiency, and lower limit detection compared to real-time PCR.

Figure 2

Sensitivity of the duplex real-time PCR and the ddPCR assay. The duplex real-time PCR assay for pB646L (A) and pEP402R (B) was evaluated at concentrations of 106-100 copies /μl template and nuclease-free water as negative control (NC). The duplex ddPCR assay for pB646L (C) and pEP402R (D) was evaluated at the same concentration condition and the template plasmid concentrations gradient (lanes 1–8) was 106, 105, 104, 103, 102, 101 and 100 copies /μl, and NC.

Specificity of the Duplex ddPCR

To evaluate the specificity of the duplex ddPCR assay, ASFV, CSFV, PRRSV, PEDV, JEV, PPV, and PCV2 were used as templates and tested by duplex ddPCR with primers and probes of B646L and EP402R. As shown in Figures 3A,C, no cross-reactivity was detected. Moreover, the positive events of ASFV B646L and EP402R molecules in other swine disease genome were zero (Figures 3B,D). These results demonstrated that both duplex ddPCR and real-time PCR assay possessed high specificity.

Figure 3

Specificity of the duplex ddPCR assay. The fluorescence amplitude of African swine fever virus (ASFV) (lane1) B646L (A) and EP402R (C) amplification was observed without cross-reactivity with CSFV, PRRSV, PEDV, JEV, PPV, PCV2 (lanes 2–7), and NC (lane 8). The ratio of B646L (B) and EP402R (D) for classical swine fever virus (CSFV), porcine reproductive and respiratory syndrome virus (PRRSV), porcine epidemic diarrhea virus (PEDV), porcine parvovirus (PPV), Japanese encephalitis virus (JEV), porcine circovirus type 2 (PCV2) (lanes 2–7), and NC (lane 8) positive events to total partitions was 0.

Application of the Duplex ddPCR Assay in Clinical Samples

Forty-four swine blood samples as well as EP402R-positive samples (inactivated GZ201801 and Δ MGF360-505R-ASFV) and EP402R-negative samples (water, ΔEP402R-ASFV, and ΔEP402R/Δ MGF360-505R-ASFV) were tested by the established duplex ddPCR and real-time PCR. Duplex ddPCR detected 44 FAM and HEX signals in all ASFV-suspicious swine blood samples (44/44, 100%) (Supplementary Table S1). However, duplex real-time PCR demonstrated 42 samples showing FAM fluorescence curves (42/44, 95%) and 41 samples showing HEX fluorescence curves (41/44, 93%). Moreover, no HEX fluorescence signal was detected in either duplex ddPCR assay or duplex real-time PCR assay. Clinical samples of B646L or EP402R showed negative results by real-time PCR, which were also detected as very low nucleic acid concentration in ddPCR, indicating low viral loads in these samples.

Discussion

African swine fever virus is a severe contagious swine disease. Strict biosecurity measures and slaughter are the main means of prevention and remediation before commercial vaccines and antiviral drugs were authorized. It caused huge economic losses, especially for a country like China that has a large-scale pig industry. CD2v, encoding by the EP402R gene, is required for protection against the challenge of the homologous ASFV strain. Several EP402R-deleted live attenuated vaccine candidates have demonstrated that CD2v is associated with virus replication and virulence in vivo, and EP402R-deleted strains can confer protection against homologous and heterologous ASFV strains (8, 9, 20, 21). In addition, ASFV strains, HLJ/HRB1/20 and HeB/Q3/20, demonstrated deletion in the EP402R gene, causing consequent loss of protein function (7). Thus, it is critical to develop an efficient detection method to differentiate not only EP402R-deleted but also mutated wild-type viruses from wild-type virus.

The PCR-based detection method is widely used in pathogen detection due to the swift, accurate, and efficient characteristics. Multifarious PCR-based detection methods have been established and applied to ASFV detection, such as real-time PCR, droplet digital PCR, loop-mediated isothermal amplification (LAMP) assay, recombinase polymerase amplification (RPA) assay, and CRISPR/Cas12a assay (16, 19, 22–25). Among them, ddPCR achieves theoretical single-molecule amplification through extreme dilution, and then uses PCR amplification and Poisson distribution to calculate the original concentration of the sample (26). Compared with real-time PCR, ddPCR possesses higher detection sensitivity without the requirement of standard products and standard curves. Meanwhile, end-point PCR signal count detection does not depend on the Ct value, and can effectively overcome the influence of PCR inhibitors. A ddPCR assay was developed based on the TaqMan probe targeting the ASFV K205R gene, showing a high degree of linearity and specificity and a 10-fold greater sensitivity compared with real-time PCR (16). Afterwards, a nanofluidic chip dPCR (cdPCR) was designed on the QuantStudio 3D digital PCR platform, which showed a 33-times higher limit detection compared with real-time PCR approved by OIE (17). Although several multiplex real-time PCR assays were developed to differentiate gene-deleted and wild-type viruses, existing ddPCR methods detected a single ASFV gene and possessed no ability to differentiate gene-deleted strains (27–29).

In this study, we developed and evaluated a duplex ddPCR assay targeting the B646L gene and EP402R gene. Two signals (FAM and HEX) were generated to represent the B646L and the EP402R gene, respectively. The B646L gene, which is a conserved region for wild-type and gene-deleted viruses, was selected as a reference gene to distinguish ASFV from other viruses. Detection of the FAM signal indicated the presence of ASFV nucleic acid, whether it was wild-type ASFV or EP402R-deleted ASFV. However, only FAM signal was detected implying that the EP402R was deleted, and further verification, sequencing for example, should be made to verify whether it is a vaccine candidate or mutant strain. The digital PCR result on 44 pig blood samples showed that no EP402R deletion virus was found, but samples with extreme low viral load demonstrated false negative results by real-time PCR assay. In addition, EP402R-deleted virus strains previously characterized (ΔEP402R/ΔMGF360-505-ASFV and ΔEP402R-ASFV) were detected as zero positive droplet by established duplex ddPCR assay, indicating a good application in EP402R-deleted and wild-type ASFV differentiation. A triplex real-time PCR assay targeting B646L gene, MGF_360-14L gene, and EP402R gene demonstrated a similar limit detection compared with the duplex ddPCR (27). Also, the limit detection of duplex ddPCR was similar with that of cdPCR assay established (30.1995 copies per reaction for B646L gene) (17). However, duplex ddPCR assay can only distinguish between wild virus strains and EP402R-deleted strains, which may lead to false test results when detecting other gene-deleted strains. Therefore, protocol-based multiple-cluster ddPCR assay is permissible to be developed in the future to differentiate other mutations associated with attenuated ASFV phenotypes from wild-type virus (30).

In summary, a specific, sensitive, and accurate duplex ddPCR assay to simultaneously detect ASFV B646L and EP402R was first established, which could distinguish EP402R-deleted viruses and wild-type viruses.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

JZ and WQ: conceptualization, WJ: methodology, YH, QG, FG, and HC: data curation. JZ: writing—original draft preparation. WQ: writing—review and editing. GZ, ML, and WQ: funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by the Key-Area Research and Development Program of Guangdong Province (No. 2019B020211003), National Natural Science Foundation of China (No. 31941014), Project of Swine Innovation Team in Guangdong Modern Agricultural Research System (No. 2020KJ126), and Guangdong Province Rural Revitalization Strategy Special Project (No. 200-2018-XMZC-0001-107-0145).

Conflict of Interest

The 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.