Impact of a Live Attenuated Classical Swine Fever Virus Introduced to Jeju Island, a CSF-Free Area

Here, we examine the effects of LOM(Low virulence of Miyagi) strains isolated from pigs (Jeju LOM strains) of Jeju Island, where vaccination with a live attenuated classical swine fever (CSF) LOM vaccine strain was stopped. The circulation of the Jeju LOM strains was mainly caused by a commercial swine erysipelas (Erysipelothrix rhusiopathiae) vaccine mixed with a LOM vaccine strain, which was inoculated into pregnant sows of 20 pig farms in 2014. The Jeju LOM strain was transmitted to 91 pig farms from 2015 to 2018. A histopathogenic investigation was performed for 25 farms among 111 farms affected by the Jeju LOM strain and revealed pigs infected with the Jeju LOM strain in combination with other pathogens, which resulted in the abortion of fetuses and mortality in suckling piglets. Histopathologic examination and immunohistochemical staining identified CSF-like lesions. Our results also confirm that the main transmission factor for the Jeju LOM strain circulation is the vehicles entering/exiting farms and slaughterhouses. Probability estimates of transmission between cohabiting pigs and pigs harboring the Jeju LOM strain JJ16LOM-YJK08 revealed that immunocompromised pigs showed horizontal transmission (r = 1.22). In a full genome analysis, we did not find genetic mutation on the site that is known to relate to pathogenicity between Jeju LOM strains (2014–2018) and the commercial LOM vaccine strain. However, we were not able to determine whether the Jeju LOM strain (2014–2018) is genetically the same virus as those of the commercial LOM vaccine due to several genetic variations in structure and non-structure proteins. Therefore, further studies are needed to evaluate the pathogenicity of the Jeju LOM strain in pregnant sow and SPF pigs and to clarify the characteristics of Jeju LOM and commercial LOM vaccine strains.


Introduction
Classical swine fever virus (CSFV), a member of the genus Pestivirus within the family Flaviviridae, comprises a single positive-stranded RNA genome of approximately 12.3 kb, which encodes a polyprotein of 3898 amino acids [1]. The viral genome comprises a 5' untranslated region (5' UTR), an N-terminal protease (N pro ), a capsid (C) protein, envelope (E) proteins (E rns , E1, E2, p7), non-structural (NS) proteins (NS2, NS3, NS4A, NS4B, NS5A, NS5B), and a 3' UTR [2]. CSFV is categorized into three genotypes, 1, 2, and 3, each of which can be subdivided into three subgenotypes (1.1-1.3, 2.1-2.3, and 3.1-3.4) [3]. Classical swine fever (CSF) is a highly contagious multisystemic hemorrhagic viral disease of domestic pigs and wild boar, which can manifest as acute, subacute, chronic, or late onset disease [4]. Vaccination is used to prevent and/or reduce the number of outbreaks of CSF and, together with other control measures, was an important factor in eradicating CSF from Holland in 1985 and from the Rivas region of Nicaragua [5]. Historically, CSF vaccines have been based on attenuated strains, i.e., lapinized Chinese and tissue culture-adapted strains. Modified live vaccines (MLVs) based on several attenuated virus strains (e.g., C-strain, Thiverval, PAV-250, GPE-, and K-strains) are used most widely. The advantages and disadvantages of MLV vaccination during an outbreak have been described in previous studies [5]. Briefly, the advantages include ease of use, low cost, and an induction of life-long immunity by a single dose. However, in the 1980s and early 1990s (the time when a lapinized Chinese strain was used in Mexico), vaccinated pigs occasionally exhibited adverse reactions; some even died [5]. The CSF vaccine virus multiplies in actively replicating cells such as fetal cells and reticuloendothelial cells; some reports show that animals vaccinated with a Chinese strain virus had more respiratory infections due to a concurrent Pasteurella multocida infection [6]. In addition, the CSF vaccine virus may induce embryonic death and myoclonia congenita when administered to pregnant sows [5]. When herds that include pregnant sows are vaccinated for the first time, there can be a fall in sow fertility in the 6 weeks following vaccination, after which it returns to normal; therefore, researchers concluded that vaccination reduces herd productivity [5]. The OIE requirements for the safety of MLVs in young animals states that vaccination should not induce a high body temperature or leukopenia and should not allow horizontal transmission. For pregnant sows, safe MLV vaccination requires no transplacental transmission and no evidence of reversion-to-virulence after passage in piglets. The CSF MLV (LOM vaccine strain) used to vaccinate pigs in South Korea since 1974 has caused abortion in some pregnant sows [7], and prolonged virus shedding by immunosuppressed pigs after vaccination [8].
The South Korean government maintained a CSF (LOM) vaccine policy to control the disease on the mainland; however, vaccination on Jeju Island (located at the southern end of Korea) was stopped in 2000 (Jeju Island declared itself a "CSF-free region"). However, pig farms of Jeju Island have suffered continuous outbreaks of the LOM vaccine strain via various routes from the mainland, raising suspicion about the safety and reversion-to-pathogenicity of commercial LOM vaccine strains. Five outbreaks of the LOM vaccine strain infection of pigs have occurred over 19 years. The first (from 2004 to 2007) occurred through feed contaminated with the commercial LOM vaccine strain present in an animal plasma protein supplement [9]. The second was caused by the delivery of incorrect vaccine material from the mainland; in 2010, one farm accidently inoculated pigs with a vaccine mixed with the LOM vaccine strain. The third and fourth occurred on one farm (2012) and two farms (2013); the route of LOM vaccine strain exposure was not revealed, but feed contaminated with the LOM vaccine strain, or an unintentional injection of an incorrect CSF (LOM) vaccine from the mainland was suspected. In the fifth case (in 2014), pregnant sows on 20 pig farms were inoculated with a commercial swine erysipelas vaccine mixed with the LOM vaccine strain. Vaccination was stopped immediately, but a total of 111 pig farms (20 in 2014, 22 in 2015, 32 in 2016, 26 in 2017, and 11 in 2018) were exposed to the Jeju LOM strains (LOM strains isolated from Jeju pigs). Farm-to-farm transmission patterns showed mainly between high-density pig farms in the Hanlim region (located in northwestern Jeju) from 2015 to 2018 and the Daejeong region (located in southwestern Jeju) from 2016 to 2018.
Here, we investigated pig-to-pig and/or farm-to-farm transmissions, and the possibility of reversion-to-pathogenicity by comparing genetic mutations in the Jeju LOM strain viruses.

Pig-to-Pig Transmission and Reproduction Rate (R)
Six pigs were inoculated with the Jeju LOM strain JJ16LOM-YJK08 and, after 24 h, were placed in a pen with six non-inoculated pigs. For Group 1, we used pigs (unhealthy) with PRRSV or PCV2, and Group 2 used pigs (healthy) without specific wasting diseases. In Group 1, two non-inoculated pigs were positive for anti-CSF (Jeju LOM) antibodies on Day 21 and Day 28. However, no non-inoculated pigs in Group 2 had detectable anti-CSF (Jeju LOM) antibodies on Day 45 ( Table 5). The transmission possibility estimate for group 1 was R0 = 1.22 (95% confidence interval (CI), 0.980-1.765), whereas that for Group 2 was R0 = 0.00 (95% CI, not applicable) ( Table 5).

Root-to-Top Divergence and Positive Selection Analyses
Application of the heuristic residual mean squared method to all strains (the commercial LOM vaccine and Jeju LOM strains) using the TempEst program revealed a slope of 3.93 × 10 −5 (rate), an X-intercept (TMRCA) of 1797.24, a correlation coefficient of 0.1272, an R squared value of 1.6173 × 10 −2 , and a residual mean squared value of 6.9076 × 10 −6 . Root-to-top divergence for Jeju LOM strain JJ07LOM-JSG02 was >0.0130. Mainly, the Jeju LOM strains isolated from the field showed a high divergence value of between 0.0070 and 0.0130 ( Figure 2). However, seven LOM vaccine strains (excluding 88LOM-Suri) and two Jeju LOM strains (JJ14LOM-WSH01 and JJ17LOM-LHH10) showed a low divergence value (0.0040-0.0060) ( Figure 2). The omega value (dN/dS) for the commercial LOM vaccine strains and the Jeju LOM strains isolated from the field (2004-2007 and 2014-2018) showed high homology with respect to the E rns , E1, E2, NS2, NS3, NS4A, NS4B, NS5A, and NS5B proteins ( Figure 3A). However, the omega value for the C and P7 proteins of commercial LOM vaccine strains was a little higher (0.34224 and 0.43781, respectively) than that of the five Jeju LOM strains isolated in 2004 to 2007. The omega values for the C, E1, and E2 structural proteins of the Jeju LOM strains isolated in the field in 2014 to 2018 were higher than for the NS proteins ( Figure 3A). Jeju LOM strains isolated from pigs from 2004 to 2007 were 98.8% to 99.2% identical at the nucleotide level and 99.2% to 99.4% identical at the amino acid level, with an omega value of 0.14088. Jeju LOM strains isolated from 2014 to 2018 showed 98.2% to 99.4% identity at the nucleotide level and 99.1% to 99.6% identity at the amino acid level, with an omega value of 0.16196 (Supplemental Table S3 Table S3).

Geographic Distance and MCC Tree Analysis
A correlation between geographic distance and genetic p-distance of the Jeju LOM strains contaminating the commercial LOM vaccine strain 16LOM-KM00 from 2014 to 2018 was confirmed,  Table S3).

Geographic Distance and MCC Tree Analysis
A correlation between geographic distance and genetic p-distance of the Jeju LOM strains contaminating the commercial LOM vaccine strain 16LOM-KM00 from 2014 to 2018 was confirmed,

Geographic Distance and MCC Tree Analysis
A correlation between geographic distance and genetic p-distance of the Jeju LOM strains contaminating the commercial LOM vaccine strain 16LOM-KM00 from 2014 to 2018 was confirmed, as shown in Figure 3B. Jeju LOM strain JJ14LOM-WSH01, isolated at a pig farm in the Aeweal region in 2014, showed a genetic p-distance from commercial LOM vaccine strain 16LOM-KM00 of 0.0018. Eleven Jeju LOM strains (Hanlim region: geographic distance, 13.5-20.1 km) and four Jeju LOM strains (Daejeong region: geographic distance, 27.4-28.1 km) isolated from pig farms showed a genetic p-distance of 0.0045-0.0139 and 0.0011-0.0086, respectively ( Figures 3B and 4). Among genotypes 1, 2, and 3 in the beast tree constructed after global analysis of the complete E2 protein of CSFV, all LOM strains (including the Jeju LOM strains) belonged to independent groups within subgenotype 1.1 ( Figure 5). From the mid-1980s, LOM strains were divided into two clusters: a lower cluster comprising commercial LOM vaccine strains and an upper cluster comprising mainly Jeju LOM strains isolated from pigs ( Figure 5). In the above cluster, Jeju LOM strains isolated from pigs on Jeju Island were as shown in Figure 3B. Jeju LOM strain JJ14LOM-WSH01, isolated at a pig farm in the Aeweal region in 2014, showed a genetic p-distance from commercial LOM vaccine strain 16LOM-KM00 of 0.0018. Eleven Jeju LOM strains (Hanlim region: geographic distance, 13.5-20.1 km) and four Jeju LOM strains (Daejeong region: geographic distance, 27.4-28.1 km) isolated from pig farms showed a genetic p-distance of 0.0045-0.0139 and 0.0011-0.0086, respectively ( Figures 3B and 4). Among genotypes 1, 2, and 3 in the beast tree constructed after global analysis of the complete E2 protein of CSFV, all LOM strains (including the Jeju LOM strains) belonged to independent groups within subgenotype 1.1 ( Figure 5). From the mid-1980s, LOM strains were divided into two clusters: a lower cluster comprising commercial LOM vaccine strains and an upper cluster comprising mainly Jeju LOM strains isolated from pigs ( Figure 5). In the above cluster, Jeju LOM strains isolated from pigs on Jeju Island were divided according to the

Discussion
The main characteristics of the LOM strain are attenuated; however, inoculation into pregnant sows lacking anti-CSF antibodies triggers abortion [7]. In addition, LOM vaccination of immunocompetent pigs leads to an antigen release period of more than 4 weeks, which results in the worsening of lesions if an animal is co-infected with other pathogens [8]. A previous study tested the effect of a live attenuated BVDV vaccine on pig-to-pig transmission and found that the R value was <1 [10]. Here, we found that transmission of the Jeju LOM strain differed significantly depending on the health status of the pig (Group 1 (unhealthy): r = 1.22; and Group 2 (healthy): r = 0). This is because, as for other wasting diseases, the period of virus release from infected pigs is prolonged and the amount of virus shed in the feces increases [8]. A commercial LOM vaccine strain showed 10% of the average transmission rate from inoculated pigs to cohabitating noninoculated pigs [11]. Although a small number of pigs were used for the transmission experiment in

Discussion
The main characteristics of the LOM strain are attenuated; however, inoculation into pregnant sows lacking anti-CSF antibodies triggers abortion [7]. In addition, LOM vaccination of immunocompetent pigs leads to an antigen release period of more than 4 weeks, which results in the worsening of lesions if an animal is co-infected with other pathogens [8]. A previous study tested the effect of a live attenuated BVDV vaccine on pig-to-pig transmission and found that the R value was <1 [10]. Here, we found that transmission of the Jeju LOM strain differed significantly depending on the health status of the pig (Group 1 (unhealthy): r = 1.22; and Group 2 (healthy): r = 0). This is because, as for other wasting diseases, the period of virus release from infected pigs is prolonged and the amount of virus shed in the feces increases [8]. A commercial LOM vaccine strain showed 10% of the average transmission rate from inoculated pigs to cohabitating non-inoculated pigs [11]. Although a small number of pigs were used for the transmission experiment in this study, the transmission characteristics of the Jeju LOM strain were revealed to be similar to the commercial LOM vaccine. In a previous study, after oral inoculation of the LOM vaccine strain of 2 × 10 4.0 TCID 50 /mL, the LOM vaccine virus was detected in blood and feces in the 9 days post-inoculation [9]. We also detected 11 samples of live Jeju LOM viruses with more than 10 4.0 TCID 50 /mL in the slaughterhouse experiment, which support the possibility of farm-slaughterhouse-farm transmission.
From 2004 to 2005, LOM-infected pigs in Jeju showed varying symptoms. The clinical symptoms showed both the similarities and differences from those associated with CSF [10]. Similarities to more current Jeju LOM strain infections (2014-2018) included skin ulcers, ulceration of the tonsil, petechial hemorrhage in the kidney, button ulceration of the ileocecal region, infarction of the spleen, and pneumonia. Differences in histologic findings were a weakness in endothelial cell lesions of organs such as the spleen, lymph nodes, and brain. These lesions seem to be more related with PRRS not CSF. In addition, meningitis was observed in the central nervous system, CSF did not show meningitis, and lymphoid organs showed little increase in volume other than lymphocytic atrophy. The Jeju LOM strain that was prevalent from 2004 to 2007 was not related to CSF pathogenicity [10]. IHC (2004)(2005) was Jeju LOM virus positive in 52.1% (25/48) of suckling piglets and 28.6% (6/21) of weaned pigs [10]. Our Jeju LOM antigen results (2014-2018) were lower than the results of the Jeju LOM antigen (2004)(2005) test (58.3%) from a previous report [10]. Jeju LOM strains (2014-2018) were also suspected to have caused CSF-like lesions in 7 young piglets (average age of 10 days), which was proved by histopathological findings and IHC analysis. It may be deduced that some CSF lesions in suckling piglets, as well as abortion during pregnancy, are associated with LOM strains. However, in this case, it is unlikely that persistent infection will continue to circulate the virus among pigs.
Some experts suggest that Jeju LOM strains detected on Jeju pig farms from 2014 to 2018 have recovered pathogenicity and that it is now different from that of the commercial LOM vaccine strain on the current market. However, gene analysis of Jeju LOM strains revealed no evidence for restoration of pathogenicity. Previous studies have identified differences in pathogenic regions between the ALD and GPE strains as E2 (T830A), NS4B (V2475A, A2563V), and N pro (N136D) [12,13]. Other studies revealed that Shimen's pathogenicity-related gene locus is within E2 (T745I, M979K) [14] and E2 (N850S) [15]. Our analysis of the amino acid sequences of the commercial LOM vaccine strain and the Jeju LOM strains detected in Jeju pigs revealed no changes at these sites. In addition, although there are differences in pathogenicity linked to changes in N-linked glycosylation sites on the E2 protein of CSF, there was no change of N-linked glycosylation sites in the Jeju LOM strains [16]. A study comparing the complete nucleotide and amino acid sequences of the ALD and GPE strains revealed 98.2% and 98.8% identity, respectively [17]. Although the Jeju LOM strain (2014-2018) harbors several amino acid mutations, it is difficult to interpret this as conferring recovery of pathogenicity or additional attenuation. There was no significant difference in the gene mutation rate among Jeju LOM strains, nor was there a significant difference in omega values. Therefore, Jeju LOM strains (2014-2018) from Jeju cannot be considered a pathogenically-reversed LOM strain, which is virulent in pigs but shows a safety profile characteristic of the commercial LOM vaccine strain described in previous study [7]. More precisely, the presence of Jeju LOM strains with recovered pathogenicity should be evaluated in animals (i.e., pregnancy sows and SPF pigs).
In conclusion, the CSF-like histopathogenic lesions of Jeju pigs revealed to be more related to other viral pathogens rather than the Jeju LOM strains (2014-2018), despite the presence of the Jeju LOM strain in organs of piglets. We also confirmed that the pig-to-pig transmission of the Jeju LOM strain and farm-to-farm transmission may have been caused by vehicles visiting the slaughterhouse. Although we found some genetic differences between the Jeju LOM strains (2014-2018) and commercial LOM vaccine strains, more pathogenesis studies may be needed using animals such as pregnant sows and SPF pigs.

Detection of Histopathologic Lesions and Immunohistochemical Straining Analysis
From 2014 to 2018, tests in this study performed by the diagnostic services of a veterinary medicine college in Jeju national university to detect animal infectious disease identified 122 samples (100 piglets and 22 fetuses) from 25 pig farms that were positive for CSFV by RT-PCR, histopathology, and IHC. The amplification targets for RT-PCR were CSFV [9], porcine circovirus type 2 (PCV2) [18], porcine reproductive and respiratory syndrome virus (PRRSV), swine influenza virus (SIV), cytomegalovirus (PCMV), Salmonella spp, Streptococcus spp, Actinobacillus, pleuropneumoniae, and Haemophilus parasuis. Histopathologic lesions specific to CSFV were examined by staining tissues with hematoxylin and eosin. The following were examined: skin (hemorrhage and endothelial damage), lymph nodes (peripheral hemorrhage, lymphoid depletion, and reticular cell hyperplasia), kidney and bladder (hemorrhage, swelling and peripheral hemorrhage), spleen (hemorrhagic infarction, endothelial damage), cecum and colon (necrotic or ulcerative colitis, vascular congestion, and sub-serous hemorrhage), heart (myocardial hemorrhage), and brain and spinal cord (nonsuppurative encephalitis, proliferation of endothelial cells, perivascular cuffing, microgliosis, and focal necrosis). To confirm the presence of CSF antigen in tissues, IHC was performed using the EnVision TM peroxidase-conjugated polymer reagent (DAKO, Denmark). In brief, tissues were reacted with a primary antibody (mouse anti-CSFV) (WH303, Animal and Plant Health Agency, New Haw, Addlestone, UK), followed by EnVision TM /HRP, rabbit/mouse (EVN) regent (DAKO, USA), and 3,3'-diamino-benzidine tetrahydrochloride (DAB). The presence of CSF antigen is denoted by a dark brown deposit in the tissue section.

Detection of Anti-CSF Antibodies
Seven pig farms with the Jeju LOM infection were selected to estimate how Jeju LOM strains were introduced into the farms and how they spread within the farms. The anti-CSF antibody-positive rates for each pig farm and antibody titer levels for pigs within farms were examined according to age. Blood samples were collected from 20 pigs in each group: sows, suckling pigs (10-20 days old), weaning pigs (40-60 days old), growing pigs (90-120 days old), and finishing pigs (150-180 days old). To detect CSF-specific neutralizing antibodies, a neutralizing peroxidase-linked assay (NPLA) was performed according to the standards manual of the World Organization for Animal Health [19]. For PK-15 cell staining, the monoclonal antibody 3B6 (Median Diagnostics, Chuncheon, Korea) was used to detect the CSF E2 protein.

Environmental Samples Taken from a Slaughterhouse
In 2017, 242 samples were obtained from a slaughterhouse in Jeju to investigate CSFV (Jeju LOM strain). These included 71 samples from the driver's foot floor, 71 from vehicle wheels, 74 from the vehicle's pig-holding compartment, and 26 from other sites. Real-time quantitative PCR (qRT-PCR) was performed to detect the CSF antigen copy number in fecal and environmental samples. The VDx ® CSFV qRT-PCR (MEDIAN Diagnostic Co. Cat No. NS-CSF-31, Gangwon-do, Korea), which uses TaqMan probes, detects the CSFV 5' UTR with high specificity; it does not detect BVDV or border disease virus, which also belong to the Pestivirus genus. Briefly, the qRT-PCR program comprised the following steps: cDNA synthesis (50 • C, 30 min) and initial inactivation (95 • C, 15 min), followed by a two-step PCR comprising 42 cycles of denaturation (95 • C, 10 s) and extension (60 • C, 60 s). A peroxidase-linked immunosorbent assay (PLA) was used to confirm viability of Jeju LOM strains. Briefly, PK-15 cells (grown to 80% confluence in 24-well plates) were inoculated for 72 h at 37 • C with 10% tissue homogenates including the Jeju LOM strain in a minimum essential medium. The PK-15 cells were fixed in pre-chilled 80% acetone after 72 h and reacted with a 3B6 monoclonal antibody specific for CSFV E2. Subsequently, the PK-15 cells were reacted with biotinylated anti-mouse IgG (H + L) (VECTOR Laboratories, Cat No. BA-9200, Burlingame, CA, USA) and ABC solution (VECTOR Laboratories, Cat No. PK-4000, Burlingame, CA, USA). After staining with DAB peroxidase substrate (VECTOR Laboratories, Cat No. SK-4100, Burlingame, CA, USA) following the manufacturer's instruction, the PK-15 cells were observed under a microscope.

Horizontal Transmission between Pigs
To investigate the transmission probability of the Jeju LOM strain among pigs, 40-day-old CSF antigen-and antibody-negative pigs from two pig farms (n = 24, 12 per farm) were used. Group 1 comprised of pigs with a wasting disease (PRRSV or PCV2) and Group 2 comprised of pigs with specific non-disease. Six pigs from each group were inoculated intramuscularly with 2 mL (10 4.0 TCID 50 /mL) of the Jeju LOM strain JJ16LOM-YJK08, which is a representative strain circulating in Jeju between 2014 and 2018. After 24 h, the six pigs from each group were housed with non-injected pigs and observed for clinical signs and symptoms. Blood was collected every week for 45 days and tested in the NPLA assay and by qRT-PCR. The estimated reproduction number (R0) and confidence interval (CI) were calculated as described previously [20] using the following formula: R0 = −In((1−AR)/S0)/(AR−(1−S0)); CI = AR ± 1.96 ARx(1 − AR)/n).

Genomic Analysis of Commercial LOM Vaccine Strains and Jeju LOM Strains
The complete genome sequences of seven commercial LOM vaccine strains (16LOM-GC00, 16LOM-JY00, 16LOM-KM00 and 16LOM-KR00 (collected from a market in 2016), 02LOM-JY00 (collected from a market in 2002), 88LOM-suri (isolated from LOM-850 strain in 1988), and LOM-850 (original master seed in 1987)) were examined. In addition, five Jeju LOM strains isolated from pigs on Jeju Island between 2004 and 2007, and 16 Jeju LOM strains isolated from pigs on Jeju Island between 2014 and 2018, were analyzed. Positive selection analysis of the complete ORF (N pro , C, E rns , E1, E2, P7, NS2, NS3, NS4A, NS4B, NS5A, and NS5B) was conducted using several models available in the BASEML and CODEML modules of the PAML 4.6 software package [21]. Different values of the non-synonymous/synonymous (dN/dS) ratio (the omega parameter), were considered in accordance with the user manual. A dN/dS ratio of <1 indicates a purifying selection, a dN/dS ratio = 1 suggests an absence of selection (i.e., neutral evolution), and dN/dS >1 indicates a positive selection. The Bayes empirical Bayes (BEB) calculation of the posterior probabilities of site classes was used to calculate the probability that a site is under positive selection pressure [22]. TempEst (formerly known as "Path-O-gen") is a tool for investigating the temporal signals and "clocklikeness" of molecular phylogenies. The contemporaneous trees (in which all sequences were collected at the same time) and dated-tip trees (in which sequences are collected at different dates) were analyzed using the TempEst v1.5.1 program, which is temporal signal estimator tool [23]. The genome sequences of the Jeju LOM strains (n = 21) and commercial LOM vaccine strains (n = 8) in this study were submitted to the GenBank under accession numbers MN558862-MN558889.

Maximum Clade Credibility Tree
Complete sequences for the E2 gene of 145 CSFV isolates available in the GenBank Database, including complete E2 gene sequences of 8 commercial LOM vaccine strains and 21 Jeju LOM strains, were used to generate a BEAST input file using BEAUti within the BEAST package v1.8.1 [24]. Rates of nucleotide substitution per site and per year, and the most recent common ancestor (tMRCA), were estimated using a Bayesian MCMC approach. Each dataset was simulated with the following options: generation = 100,000,000, burn-in of 10%, and ESSs > 200. The exponential clock and expansion growth population model in the BEAST program was used to obtain the best-fit evolutionary model. The resulting convergence was analyzed using Tracer 1.5 [25]. Trees were summarized as a maximum clade credibility (MCC) tree using TreeAnnotator 1.7.4 [26] and visualized using Figtree 1.4 [27]. For each tree node, estimated divergence times and 95% HPD intervals, which summarize statistical uncertainties, were indicated.