key: cord-0003671-t6cnjwpd
authors: Sung, Ming-Hua; Lin, Chao-Nan; Chiou, Ming-Tang; Cheng, I-Ju; Thanh, Quang-Hien; Chao, Day-Yu; Lan, Yu-Ching
title: Phylogeographic investigation of 2014 porcine epidemic diarrhea virus (PEDV) transmission in Taiwan
date: 2019-03-06
journal: PLoS One
DOI: 10.1371/journal.pone.0213153
sha: f226ecb8f7484ca7943245de562c6f2f08be80a4
doc_id: 3671
cord_uid: t6cnjwpd

The porcine epidemic diarrhea virus (PEDV) that emerged and spread throughout Taiwan in 2014 triggered significant concern in the country’s swine industry. Acknowledging the absence of a thorough investigation at the geographic level, we used 2014 outbreak sequence information from the Taiwan government’s open access databases plus GenBank records to analyze PEDV dissemination among Taiwanese pig farms. Genetic sequences, locations, and dates of identified PEDV-positive cases were used to assess spatial, temporal, clustering, GIS, and phylogeographic factors affecting PEDV dissemination. Our conclusion is that S gene sequences from 2014 PEDV-positive clinical samples collected in Taiwan were part of the same Genogroup 2 identified in the US in 2013. According to phylogenetic and phylogeographic data, viral strains collected in different areas were generally independent of each other, with certain clusters identified across different communities. Data from GIS and multiple potential infection factors were used to pinpoint cluster dissemination in areas with large numbers of swine farms in southern Taiwan. The data indicate that the 2014 Taiwan PEDV epidemic resulted from the spread of multiple strains, with strong correlations identified with pig farm numbers and sizes (measured as animal concentrations), feed mill numbers, and the number of slaughterhouses in a specifically defined geographic area.

Porcine epidemic diarrhea virus (PEDV) causes acute diarrhea, vomiting, and dehydration, resulting in high mortality rates for suckling piglets [1] . Since 2010, a new PEDV variant belonging to Genogroup 2 has spread throughout the United States and across multiple Asian countries, including China and Taiwan [2, 3] . An initial identification in the US was made in April of 2013 [4] . A 40.5% premises-level incidence of PEDV caused the deaths of more than 8 a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 million newborn piglets in a single year-an event that significantly affected the American swine industry [5] [6] [7] [8] [9] . Subsequent Genogroup 2 epidemics have been reported in major swineproducing countries such as Canada, Mexico, Taiwan, Korea and Japan [8, [10] [11] [12] [13] . One of the most notable PEDV outbreaks occurred in South Korea in late 2013 [14] , more than a half million pigs died from PEDV infections in Japan between 2013 and 2016 [15] , and a significant increase in PEDV outbreaks occurred in Taiwan around the same time [16] . PEDV is now considered the world's most catastrophic swine disease, with major financial impacts noted throughout the global pork industry.

The PEDV genome is comprised of at least seven open reading frames (ORF1a, ORF1b and ORF2-6) encoding four structural (S), envelope (E), membrane (M), and nucleocapsid (N) proteins [17] . A high degree of genetic diversity has been observed in the S glycoprotein gene [18] [19] [20] . Partial spike (S) polyprotein genes located in the virus envelope are central to PEDV biological properties such as interactions with cellular receptors during virus entry, the neutralizing of antibody induction in natural hosts, growth adaptation in vitro, and virulence attenuation in vivo [19] . The PEDV spike (S) protein is a type 1 transmembrane envelope glycoprotein with a 4,158 nucleotide sequence divided into S1 and S2 domains. The S1 region (aa 26-734) is responsible for viral binding, while the S2 domain (aa 735-1383) serves as an anchor for viral membrane and fusion activity [21, 22] . Thus, the S glycoprotein is considered a primary target for PEDV vaccine development. As the major envelope glycoprotein found in virion, S serves as an important viral component for studying genetic relationships among PEDV isolates, and for determining PEDV epidemiological status [23, 24] .

PEDV is believed to infect pigs by both direct and indirect fecal-oral routes. Due to the scales and complexities of modern swine production systems, PEDV is likely transmitted between farms via diarrheic feces or vomitus; contaminated environmental sources involving clinically or sub-clinically infected pigs; trailers used to transport livestock, manure, or food sources; farmers or visitors wearing contaminated clothes; or wild animals and birds [5, 25] . Other potential sources include contaminated fomites (e.g., raw food, feed, sow milk), food ingredients or additives, and environmental features such as wind direction, farm altitude, terrain slope, and tree coverage [26] . After an initial outbreak, PEDV may spread at an increasingly rapid rate due to inadequate farm hygiene management procedures such as improper disinfection and poor biosecurity. The virus can remain dormant in weaning pigs or growth finishing units, eventually triggering mild symptoms and resulting in low mortality rates [6] .

Although researchers believe that PEDV infections primarily result via fecal-oral routes, the rapid regional spread of the disease raises the possibility of airborne transmission [4] . Support for this hypothesis includes an identified correlation between disease-spread direction and prevailing wind direction [7] , with environmental features such as land coverage, altitude, and slope possibly influencing airborne disease dissemination [26] . To determine specific temporal and geographic relationships associated with PEDV strain transmission, we used phylogenetic, phylodynamic and phylogeographic methods to systematically evaluate potential temporal and spatial transmission routes among Taiwanese swine farms during the 2014 outbreak.

Epidemiological and geographic data were collected from 92 animal lab reports of PEDV viral infections involving 8,557 pig farms, 37 pig feed mills, 57 slaughterhouses, and 5,806,237 animals. These reports are available from an open database maintained by the Taiwan government. Additional epidemiological and genetic information was gathered for purposes of determining details of the disease spread. A total of 48 global PEDV whole genome sequences (S1 Table) and 49 Taiwan partial S1 gene sequences (648 nucleotide of PEDV S1 gene position 1468-2115) (S2 Table) were downloaded from GenBank, and information for the TW4 whole genome sequence was collected from a previous study [27] . Information datasets focused on pig feeding and disease were collected from the Taiwan open data website (http:// data.gov.tw/), the Animal Health Research Institute of the Taiwan Executive Yuan's Agricultural Council (http://eng.nvri.gov.tw/fmodule/Default.aspx), and GenBank (https://www.ncbi. nlm.nih.gov/genbank/).

In the first stage of this study, PEDV phylogenetic and phylogeographic data analyses were performed for purposes of organizing viral transmission evidence and tracking possible transmission routes. In the second, data for variables of interest associated with the pig feed industry were collected and combined with geographic information system (GIS) data for investigation using open source Quantum GIS (QGIS v3.2.3) software [28] . An SAS MIXED procedure was used to create linear regression models, with livestock breeding variables employed to predict PEDV infections.

Gene sequences for the 49 global PEDV whole genomes and 49 Taiwan PEDV partial S1 gene sequences were downloaded from GenBank and aligned with a single Taiwan PEDV complete genome from a previous study [27] using the ClustalW multiple alignment feature in BioEdit [29] . The best-fit model, as determined using jModelTest 2.1.7 [30] , had a gamma distribution (GTR+G+I). To ensure topological consistency, phylogenetic trees were constructed using maximum likelihood (ML) and Bayesian methods (MEGA version 6.0 and BEAST v1.8.2, respectively). Branch support was evaluated using bootstrap analyses based on 1,000 ML tree replications. Bootstrap values >75% were considered as belonging to the same monophyletic group.

To identify the specific locations of migration events, we grouped the 49 PEDV isolates into different counties and used various PEDV viral strains as references. Spatial location reconstruction and viral migration activity were estimated using discrete coalescent tree and Bayesian phylogeographic methods. Bayesian Markov Chain Monte Carlo (MCMC) sampling using BEAST v1.8.2 [31] was employed to infer the time-scaled phylogenies of partial PEDV S genes. HKY+G and relaxed clock exponential models were used prior to setting coalescent population constants in the MCMC simulations. Estimated convergence and effective sampling sizes were visually assessed using Tracer v1.6. Multiple chains were combined based on a 10% burnin using the version of LogCombiner (v1.8.2) included in the BEAST package. Maximum clade credibility trees with temporal and spatial annotations were summarized with the 10% burn-in removed using TreeAnnotator (v1.8.2, also in the BEAST package). FigTree (v1.4.2) was used to generate presentation figures.

Bayes factor (BF) tests were conducted to build statistical support for transmission routes among geographic locations using SPREAD3 (v0.9.6; BF cutoff = 3) [32] . BF values were used to indicate differences between posterior and prior probabilities so that rates between any two locations were non-zero. Routes with high BF values were considered as having greater potential for viral strain migration.

To create animations of viral dispersion over time, annotated MCC trees were converted into a keyhole markup language (KML) file using SPREAD3 (v0.9.6). The KML file can be visualized with an open-access Earth map downloaded from Natural Earth (http://www. naturalearthdata.com) as a QGIS software base layer.

The maximum-likelihood phylogenetic tree (Fig 1) was constructed from 49 global PEDV whole genome sequences. Sequences found in Taiwan clustered with 100% bootstraps containing viruses from other Asian countries and the US between 2011 and 2015. A correlation was determined between this cluster and the Genogroup 2 originally identified in the US. Among the Taiwan sequences, some of the viruses collected and identified in 2014 were strongly correlated with one another (bootstrap values >75%), while others were mixed with viral strains collected in other countries at different times, suggesting multiple transmission events (Fig 2 &  Table 1 ).

To further evaluate transmission periods and to verify viral strain origins, all 49 partial S gene sequences determined from the Taiwan PEDV samples were used for phylodynamic data analyses (Fig 3) . Four clusters (A, B, C and D) exhibited statistically significant posterior probabilities (>0.8), including 9 viral strains on 4 monophyletic branches. The most recent common ancestor (MRCA) for the 9 viral strains was traced to July 2013. All other coalescent tree sequences were determined as independent. Viral branches had MRCAs in 2013 or the spring of 2014, with none identified as statistically significant. Thus, the phylodynamic coalescent tree established for this study indicates the involvement of multiple virus strains in the 2014 PEDV outbreak.

Results from our phylogeographic analysis were matched to a map of the country to specify the geographic boundaries of the Taiwan outbreak (Fig 4) . Next, molecular sequence data were combined with isolation time data and geographic coordinates to determine the spatiotemporal distribution of Taiwan PEDV strains. Lineages were identified in several agricultural communities in the far south, with additional virus strains found in central Taiwan at approximately the same time. The southern infections were spatially closer to each other. The data indicate an absence of natural and artificial barriers restricting the spread of the virus.

The pig feed industry risk map shown was created to assist in the identification of significant risk factors associated with infections ( Fig 5) . Results from our GIS system analysis of positive cases indicate correlations between PEDV infection transmission and both pig farm size (number of pigs, not physical size) (0.000232, 95% CI [0.000102, 0.000362]) and slaughterhouse distribution (0.8043, 95% CI [0.4351, 1.1735]) ( Table 2 ). Our data also indicate higher infection rates in counties with fewer pig feed mills (-0.7857, 95% CI [-1.2298, -0.3415]).

Phylogeographic inferences are a potential tool for identifying the transmission and dissemination routes of PEDV and other potentially much deadlier infectious diseases. However, to date very few research efforts in Asia have utilized full genome sequencing for determining geographic structures due to the high costs and enormous amounts of computational time Phylogeographic investigation of 2014 porcine epidemic diarrhea virus transmission in Taiwan required for analyses [33, 34] . Several researchers have suggested using partial S genes for phylogenetic tree construction and for phylodynamic analyses specifically aimed at studying the genetic relatedness of PEDV strains [3, 10, 19] . In this study, we investigated temporal and geographic relationships among PEDV strains identified as having been transmitted among farms in Taiwan in 2014, using partial sequences from S genes extracted from porcine fetus samples and GenBank sequences as reference panels.

According to the phylogenetic tree we constructed based on these partial S gene sequences, the primary PEDV strain in Taiwan is related to the Genogroup 2 strain identified in samples collected in the US in 2013 [35] . As previously suggested, the most recently identified Taiwan PEDV strains have greater similarity with US strains than with Chinese or earlier Taiwanese strains [3, 36] . Transmission may have occurred as early as December 2013 [36] , but the crossborder route to Taiwan remains unknown. Results from our phylogenetic analysis of PEDV viruses associated with the Taiwan outbreak confirm independence in multiple counties. Although these viruses share a common ancestry with the US Genogroup 2 PEDV, our coalescent tree data indicate that only 9 of the Taiwanese viruses were significantly clustered (Fig 3) . Further, no recombination involving Taiwanese strains was observed in the present study.

Five independent PEDV strains were identified in Taiwan on or before September 16, 2013 in Wanda (KM246707 and KM246708) and Jhutian (KM246673, KM246674 and KM246675) townships. These are considered starting points for PEDV complex dissemination throughout sections of Taiwan up to May 15, 2014. According to a mix of phylogenetic and coalescent data ( Table 2) , clustering was limited to 3 counties in southern Taiwan (Wanluan, Jiouru and Daliao) (Fig 3) , with most of the identified viruses existing independently. Combined, data for PEDV divergence during the Taiwan outbreak suggest a common ancestry shared by multiple virus lineages. Since that time period, no evidence has been found indicating outbreaks involving multiple virus strains in Taiwan. According to our phylogenetic analysis, most of the identified PEDV strains from the 2014 outbreak were independent, despite sharing a common ancestor. Multiple PEDV invasions from abroad were also identified in Japan during the 2014 outbreak [15, 37] . The widespread dissemination that followed presumably resulted from the movement of pigs, agricultural vehicles, farmers, farm visitors, commercial feed products, and other materials [15, 37] . Reproducing novel approaches used in the molecular and spatial surveillance of the porcine reproductive and respiratory syndrome virus (PRRSV) in the US and in PEDV studies in Japan [26, 38] , we utilized a combination of phylogeographic and GIS approaches in our effort to profile the Taiwan PEDV outbreak. GIS has been used to investigate correlations between diseases and factors that include pig farm size, number of feed mills, and number of pig slaughterhouses in a specified geographic zone (Fig 5) . We used a mixed linear regression to identify factors associated with the number of PEDV cases in Taiwan. Results indicate positive correlations between the number of cases and both slaughterhouse number and pig farm size, and a negative correlation with number of feed mills ( Table 1 ). The highest concentration of PEDV cases was identified in a multi-county section of southern Taiwan characterized by a large number of pig farms with high animal densities. Our results are in agreement with findings from previous studies suggesting that excessive farm capacity (measured as the total number of pigs on a farm) is a risk factor for the spread of PEDV [5] . Further, aerosol transmission is considered a viable dissemination route in environments marked by high pig densities and close animal proximities [4] . Regarding the negative correlation between PEDV cases and number of feed mills, our data indicate that pig farms at the end of feed routes likely have higher probabilities of infections. Some reports suggest that vehicles used for the dual purposes of transporting swine to slaughterhouses and delivering feed to farms may increase the potential for PEDV due to feedbag contamination [5, 25, 39, 40] . Our finding of a high correlation between the number of PEDV cases and the number of pig slaughterhouses suggests that transport trucks may have been a factor in the 2014 Taiwan outbreak.

Other risk factors requiring further research include spray-dried porcine plasma (SDPP, an important blood-based component of nursery pig diets) and improper disposal procedures when pig corpses are collected and sent to rendering plants. There is a clear need to collect more data on feed truck routes, SDPP supplement distribution, rendering plant procedures (especially delivery), and slaughterhouse processes when trying to identify the sources of various strains in Taiwan. Due to the potential for significant financial losses, there is a strong need to act on these and other possible factors before detailed studies can be designed, funded, and completed. Farms, slaughterhouses, and feed suppliers in counties with high pig densities need to immediately enhance their biosecurity measures to prevent future PEDV outbreaks, and greater effort is required to monitor potential transmission routes. Table. https://doi.org/10.1371/journal.pone.0213153.g004

Supporting information S1 

Conceptualization: Yu-Ching Lan. 

Porcine Epidemic Diarrhea, Diagnosis, and Elimination

Introduction of PEDV outbreak and prevention method

Evolutionary characterization of the emerging porcine epidemic diarrhea virus worldwide and 2014 epidemic in Taiwan. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases

Evidence of infectivity of airborne porcine epidemic diarrhea virus and detection of airborne viral RNA at long distances from infected herds

Epidemiological factors associated to spread of porcine epidemic diarrhea in Japan. Preventive veterinary medicine

Porcine epidemic diarrhea virus: An emerging and re-emerging epizootic swine virus

A Porcine Epidemic Diarrhea Virus Outbreak in One Geographic Region of the United States: Descriptive Epidemiology and Investigation of the Possibility of Airborne Virus Spread

Deadly pig virus slips through US borders

Risk factors for porcine reproductive and respiratory syndrome virus infection and resulting challenges for effective disease surveillance

Isolation and characterization of porcine epidemic diarrhea viruses associated with the 2013 disease outbreak among swine in the United States

Fighting a deadly pig disease. Industry, veterinarians trying to contain PED virus, new to the US

Comparison of porcine epidemic diarrhea viruses from Germany and the United States

Novel porcine epidemic diarrhea virus variant with large genomic deletion, South Korea. Emerging infectious diseases

Isolation of Porcine Epidemic Diarrhea Virus during Outbreaks in South Korea

Phylogeographic investigation of 2014 porcine epidemic diarrhea virus transmission in Taiwan PLOS ONE

Isolation and molecular characterization of porcine epidemic diarrhea viruses collected in Japan in 2014

US-like strain of porcine epidemic diarrhea virus outbreaks in Taiwan

Porcine epidemic diarrhoea virus: a comprehensive review of molecular epidemiology, diagnosis, and vaccines

Novel Porcine Epidemic Diarrhea Virus (PEDV) Variants with Large Deletions in the Spike (S) Gene Coexist with PEDV Strains Possessing an Intact S Gene in Domestic Pigs in Japan: A New Disease Situation. PloS one

Heterogeneity in spike protein genes of porcine epidemic diarrhea viruses isolated in Korea

Detection and molecular diversity of spike gene of porcine epidemic diarrhea virus in China

The S2 glycoprotein subunit of porcine epidemic diarrhea virus contains immunodominant neutralizing epitopes

Cellular entry of the porcine epidemic diarrhea virus

Sequence analysis of the partial spike glycoprotein gene of porcine epidemic diarrhea viruses isolated in Korea

Coronaviruses: structure and genome expression

Role of transportation in spread of porcine epidemic diarrhea virus infection, United States. Emerging infectious diseases

Land altitude, slope, and coverage as risk factors for Porcine Reproductive and Respiratory Syndrome (PRRS) outbreaks in the United States

Molecular Characterization of the Porcine Epidemic Diarrhea Virus TW4/2014 in Taiwan. Austin Virology and Retro Virology

QGIS Geographic Information System. Open Source Geospatial Foundation

BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT

jModelTest 2: more models, new heuristics and parallel computing

Bayesian phylogeography finds its roots

SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics

Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States

Porcine epidemic diarrhea virus infection: Etiology, epidemiology, pathogenesis and immunoprophylaxis

Distinct Characteristics and Complex Evolution of PEDV Strains

US-Like Strain of Porcine Epidemic Diarrhea Virus Outbreaks in Taiwan

Molecular characterization of pig epidemic diarrhoea viruses isolated in Japan from 2013 to 2014. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases

Spatial dynamics of porcine epidemic diarrhea (PED) spread in the southern Kyushu, Japan. Preventive veterinary medicine

An evaluation of contaminated complete feed as a vehicle for porcine epidemic diarrhea virus infection of naive pigs following consumption via natural feeding behavior: proof of concept

An evaluation of a liquid antimicrobial (Sal CURB(R)) for reducing the risk of porcine epidemic diarrhea virus infection of naive pigs during consumption of contaminated feed