Detection and Genetic Diversity of Porcine Coronavirus Involved in Diarrhea Outbreaks in Spain
Héctor Puente
Héctor Argüello 
Óscar Mencía-Ares
Manuel Gómez-García
Ana Carvajal- Department of Animal Health, Faculty of Veterinary Medicine, Universidad de León, León, Spain
Porcine enteric coronaviruses include some of the most relevant viral pathogens to the swine industry such as porcine epidemic diarrhea virus (PEDV) or porcine transmissible gastroenteritis virus (TGEV) as well as several recently identified virus such as swine enteric coronavirus (SeCoV), porcine deltacoronavirus (PDCoV) or swine enteric alphacoronavirus (SeACoV). The aim of this study is the identification and characterization of enteric coronaviruses on Spanish pig farms between 2017 and 2019. The study was carried out on 106 swine farms with diarrhea outbreaks where a viral etiology was suspected by using two duplex RT-PCRs developed for the detection of porcine enteric coronaviruses. PEDV was the only coronavirus detected in our research (38.7% positive outbreaks, 41 out of 106) and neither TGEV, SeCoV, PDCoV nor SeACoV were detected in any of the samples. The complete S-gene of all the PEDV isolates recovered were obtained and compared to PEDV and SeCoV sequences available in GenBank. The phylogenetic tree showed that only PEDV of the INDEL 2 or G1b genogroup has circulated in Spain between 2017 and 2019. Three different variants were detected, the recombinant PEDV-SeCoV being the most widespread. These results show that PEDV is a relevant cause of enteric disorders in pigs in Spain while new emerging coronavirus have not been detected so far. However, the monitoring of these virus is advisable to curtail their emergence and spread.
Introduction
Coronaviruses (CoVs) are found in a wide variety of animals including both mammals and birds in which they cause a variety of respiratory, enteric or even hepatic and neurological disorders (1). They belong to the Coronaviridae family which recognizes four genera based on phylogenetic clustering: Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus. The CoVs are enveloped viruses, and their genome is composed of a non-segmented, positive sense and single-stranded RNA with a size of ~30 kb. From the 5′ end to the 3′ end, their genomic structure is made up of at least six open reading frames (ORFs) named ORF1a, ORF1b, spike (S), envelope (E), membrane (M), and nucleocapsid (N). ORF1a and ORF1b encode non-structural polyproteins, whereas S, E, M, and N genes encode the corresponding structural proteins (2).
Two species of the Alphacoronavirus genus, transmissible gastroenteritis virus (TGEV) and the porcine epidemic diarrhea virus (PEDV) have long been recognized as the cause of acute diarrhea outbreaks on swine farms affecting pigs of all ages and causing high mortality in lactating piglets. The relevance of TGEV on farms is scarce (1), mainly due to the widespread distribution of a respiratory mutant of TGEV, the porcine respiratory coronavirus (PRCV), which partially protects animals against the enteric disease. PEDV was recognized for the first time in Europe and Asia during the seventies and the eighties, respectively (3). In Europe its incidence markedly decreases in the nineties and subsequent years while in Asia PEDV has remained as a major cause of diarrhea outbreaks until now. Moreover, highly pathogenic variants of PEDV have been described in Asia since 2010. This virus emerged in America in 2013–2014 causing substantial economic losses (4). PEDV also re-emerged in Europe soon after its first description in the USA and PEDV outbreaks have been reported in several European countries since 2014 (5–7). Two main PEDV genogroups, named INDEL or G1 and non-INDEL or G2 have been described and differentiated by insertions-deletions in the S1 subunit of the S-gene (8, 9). Non-INDEL isolates have been associated with a higher virulence and better horizontal transmission (10). Both genotypes are reported on infected farms in Asia and America, while in Europe there is no evidence of the presence of the non-INDEL genogroup with the only exception of an Ukrainian isolate (11).
New coronaviruses affecting pigs have been unveiled in recent years. A chimeric virus produced by the recombination of TGEV/PRCV (backbone sequence) with PEDV (S gene), called swine enteric coronavirus (SeCoV), was reported in several European countries including Spain (12–15). A porcine deltacoronavirus (PDCoV) was detected in 2012 in Hong Kong (16) and subsequently on pig farms from the USA, Canada and several Asian countries (17). And finally, a new Alphacoronavirus named swine enteric alphacoronavirus (SeACoV), also known as swine acute diarrhea syndrome coronavirus (SADS-CoV) or porcine enteric alphacoronavirus (PEAV), was identified as the cause of severe diarrhea in neonatal piglets in China (18–22).
The emergence of these new CoVs and re-emergence of PEDV in Europe, immediately after its disruptive appearance in America, requires studies characterizing these CoVs on pig farms so as to allow for an accurate differential diagnosis of viral diarrhea. This study aims at disclosing the current situation of enteric CoVs in Spain, the largest pig producer in Europe, through the identification and characterization of CoVs in diarrhea outbreaks on Spanish pig farms between 2017 and 2019.
Materials and Methods
Samples Collection and Preparation
The study was conducted between January 2017 and March 2019 on 106 swine farms (105 from Spain and one from Portugal) with diarrhea outbreaks in which a viral etiology was suspected. The outbreaks affected nursing piglets (<21 days) (28 farms), postweaning-growing pigs (21–70 days) (17 farms), or fattening pigs (>70 days) (61 farms). Location of the farms include 22 provinces in the northeast, center and northwest of Spain (Figure 1). Fecal samples were submitted for diagnostic purposes to the Infectious Diseases Unit of the Animal Health Department of the University of León (Spain). From each farm, two to six individual fecal samples were submitted. Individual fecal samples were mixed to prepare one pooled sample per farm.
Figure 1. Map showing the distribution of the viral suspected diarrhea outbreaks investigated in this research (shaded area). The number of farms sampled (blue circle) and PEDV positive farms per province (orange circle) are given.
Pooled samples were diluted 1:2 (v/v) in sterile phosphate buffered saline (PBS), homogenized by vortex mixing and centrifuged for 10 min at 20,000 g. The RNA was extracted from 140 μl of the supernatant using QIAMP Viral RNA Mini Kit (QIAGEN), following the manufacturer's instructions.
Molecular Diagnosis of Porcine Enteric Coronaviruses
Two duplex RT-PCRs were developed for the detection of SeACoV (ORF1ab), TGEV/SeCoV (N gene), PEDV/SeCoV (S gene) and PDCoV (N gene) (Table 1). Two conventional RT-PCRs were also developed to confirm SeCoV by excluding a potential PEDV-TGEV co-infection, by detecting the M gene of PEDV and the S gene of TGEV. The RT-PCR was carried out with Verso 1-Step RT-PCR ReddyMix kit (Thermo Scientific). The reaction was conducted under the following conditions: 50°C for 30 min, 95°C for 2 min, 40 cycles at 95°C for 20 s, 50°C for 30 s, and 72°C for 1 min, followed by a final extension step at 72°C for 10 min.
Table 1. Primer sets used for the detection of porcine enteric coronaviruses by multiplex RT-PCR and for the amplification of the S gene of PEDV.
The RT-PCR products were visualized on a 1.5% agarose gel containing RedSafe Nucleic Acid Staining Solution (iNtRON Biotechnology, Inc.). The length of the PCR fragment generated for each CoV is shown in Table 1.
PEDV S-gene Sequencing
The S-gene of PEDV positive samples was amplified using four overlapping fragments with the primers described in Table 1 and the Verso 1-Step RT-PCR ReddyMix kit (Thermo Scientific). The reaction was conducted under the following conditions: 50°C for 30 min, 95°C for 2 min, 45 cycles at 95°C for 20 s, 50°C for 30 s, and 72°C for 2 min, followed by a final extension step at 72°C for 10 min. The RT-PCR products were purified using the GeneMATRIX Basic DNA Purification Kit (EurX). The complete sequences of the S-gene were obtained by using forward and reverse Sanger sequencing. The complete sequence of the S gene of 36 PEDV isolates from different farms can be accessed at the NCBI GenBank with the accession numbers MW251343-MW251378 (Table 2). The remaining five PEDV isolates included in this research were previously sequenced by using a RNA virus-specific tailor-made NGS protocol (15).
Table 2. List of porcine epidemic diarrhea virus (PEDV) isolates recovered in this study.
Phylogenetic Analysis
PEDV and SeCoV genome sequences available in the GenBank database were aligned together with S-gene sequences obtained in this study using CLUSTALW. After the alignment, the evolutionary relationships among sequences were analyzed with a phylogenetic analysis, using the neighbor joining method and the maximum composite likelihood method with MEGAX software (28).
Statistical Analysis
The Fisher exact test was used to compare the frequency of occurrence of PEDV positive outbreaks among age groups and provinces (only those with five or more submitted samples were included in the analysis). ANOVA test was used to compare the number of investigated outbreaks as well as the percentage of PEDV positive outbreaks among the different trimesters of the year. Epi Info™ was used for data analysis at α = 0.05.
Results
Prevalence of Enteric Coronaviruses in Porcine Diarrhea Outbreaks
Neither TGEV, SeCoV, PDCoV nor SeACoV were detected in any of the samples, while PEDV was the only coronavirus detected in 41 out of the 106 investigated outbreaks (38.7%). Most of these outbreaks occurred in fattening pigs (24 positive farms out of 61, 39.3%) or postweaning-growing pigs (nine positive farms out of 17, 52.9%). PEDV was involved to a lesser extent in diarrhea outbreaks affecting nursing piglets (eight positive farms out of 28, 28.6%), although no significant differences in the number of PEDV confirmed outbreaks between age groups arose when compared using the Fisher exact test (p = 0.094).
PEDV positive farms were detected throughout the sampled area in northeast, center, and northwest of the country with at least one positive outbreak in 10 out of 22 sampled provinces (Figure 1). No significant differences were found in the number of PEDV confirmed outbreaks between provinces (p = 0.286).
In addition, it was evidenced that most of the investigated outbreaks occurred during the first trimester of the year (p = 0.041) (Figure 2). Although most of the PEDV positive outbreaks also occurred during the first 3 months of the year, no significant differences arose when compared using ANOVA test (p = 0.097).
Figure 2. Distribution of the viral suspected diarrhea outbreaks and PEDV positive outbreaks from January, 2017 to March, 2019. The number of sampled farms and PEDV-positive farms is shown above the columns. The percentage of PEDV-positive outbreaks in each trimester studied. The percentage of PEDV-positive outbreaks in each trimester studied is showed.
Phylogenetic Analysis Based on the Nucleotide and Amino Acid Sequences of PEDV S Gene
The full-length S-gene sequences of all PEDV positive samples (41) were compared to previous and recent sequences of PEDV and SeCoV available in GenBank (Figure 3, Supplementary Data 1). The phylogenetic tree showed that all PEDV isolates recovered in Spain between 2017 and 2019 were allocated within the INDEL 2 or G1b genogroup together with several recent European PEDV isolates and isolates from the USA and Asia. They clustered in a branch clearly separate from the non-INDEL or G2 genogroup as well as from the original European or Asian PEDV isolates included in the INDEL 1 or G1a genogroup and SeCoV isolates.
Figure 3. Phylogenetic tree based on the complete S-gene sequences including enteric porcine coronaviruses available in GenBank. Numbers along the tree represent the confidence value for a given internal branch based on 500 Bootstrap replicates (only values >70 are shown). The symbols above the strains highlight the Spanish PEDV isolates of this study. The filled circles identified the isolates sequenced in this research while the filled triangles identified isolates previously sequenced by de Nova et al. (15). GenBank accession number, country and year of the outbreak are also shown below the strains. The genogroups and subgroups referred in the text are included on the right of the tree. Scale bars indicate nucleotide substitutions per site.
Three subgroups or clusters were identified from Spanish PEDV isolates identified as INDEL 2.1, 2.2, and 2.3 (Figures 3, 4). The first was formed using two Spanish isolates recovered in 2018 and 2019 together with other Spanish and European PEDV isolates from 2014 to 2016. The INDEL 2.2 cluster included Spanish isolates from 2014 to 2019 and several Asian and American PEDV strains from the same time period. It is worth mentioning that all the isolates identified in this research included in the INDEL 2.2 subgroup correspond to isolates recovered from farms located within the same region and belong to the same pig producing company. Finally, the INDEL 2.3 cluster include recent Spanish isolates from 2017 to 2019 distributed throughout the sampling area together with Hungarian, Slovenian, Dutch, German, and French coetaneous strains. This clade corresponds to a PEDV-SeCoV recombinant variant previously described in Hungary and Slovenia (29, 30) as well as in Spain (15).
Figure 4. Phylogenetic tree based on the complete S-gene sequences including all Spanish enteric porcine coronaviruses available in GenBank. Numbers along the tree represents the confidence value for a given internal branch based on 500 Bootstrap replicates (only values >70 are shown). The symbols above the strains stand out the Spanish PEDV isolates of this study. The filled circles identified the isolates sequenced in this research and the filled triangles identified the isolates previously sequenced by de Nova et al. (15). GenBank accession number, province and year of the outbreak are also shown below the strains. The genogroups and subgroups referred in the text are included on the right of the tree. Scale bars indicate nucleotide substitutions per site. Black, pre-2000 isolates; Red, 2014 isolates; Gray, 2015 isolates; Purple, 2016 isolates; Yellow, 2017 isolates; Green, 2018 isolates; Blue, 2019 isolates.
Three major regions of the S1 gene were further analyzed and characterized at the amino acid level. Compared with three different strains of PEDV (CV777 accession no. AF353511, OH851 accession no. KJ399978 and SLOreBAS-1 accession no. KY019623), four amino acid mutations were found in four out of 41 Spanish PEDV isolates (Supplementary Data 1).
Discussion
The recent emergence of several novel porcine enteric coronaviruses such as PDCoV or SeACoV together with the re-emergence of PEDV, a classical coronavirus, and their spread throughout the main pig producing countries over the last years have highlight porcine enteric coronavirus. They have caused significant losses to the pig-farming industry associated to high morbidity acute diarrhea in pigs of all ages and high mortality in neonatal pigs in naive farms. Enteric diseases caused by porcine enteric CoVs are clinically indistinguishable thus making differential diagnosis in the laboratory an essential tool.
A hundred and six viral suspected diarrhea outbreaks were investigated between 2017 and 2019 and confirmed that PEDV was the only enteric coronavirus detected in swine farms in Spain. PEDV was identified in about a 40% of the investigated outbreaks, confirming the re-emergence of PEDV in the Iberian Peninsula as it has been described in several European countries (31) after its emergence in 2013 in the United States. The difference in the percentage of PEDV positive outbreaks among the different trimesters of the year was near to statistical significance, with most of the PEDV positive outbreaks occurring in winter (68.3%, 28 out of 41 PEDV positive outbreaks), when low temperatures favor the environmental resistance of the virus facilitating its indirect transmission. This results confirms the seasonal distribution of the disease (3, 32). Besides, although no significant differences were shown in the proportion of PEDV outbreaks between age groups, most of the outbreaks were observed among postweaning-growing or fattening pigs. The fact that few PEDV outbreaks were detected in suckling piglets (28.6%, eight PEDV outbreaks out of 28 investigated) could be a consequence of maternal immunity in the sows or high biosecurity level in farrowing facilities. Clearance of maternal antibodies together with the mix of piglets after weaning could explain a higher percentage of positive outbreaks (52.9%, nine PEDV positive outbreaks out of 17 investigated) in postweaning pigs (21–70 days) (1).
Neither TGEV, SeCoV, PDCoV nor SeACoV were detected in any of the investigated diarrhea outbreaks. To our knowledge, this is the first study in Europe actively researching PDCoV or SeACoV on swine farms. While SeACoV has a limited geographic distribution and has only been detected in swine farms in China (18, 19), PDCoV has been reported in the USA, South Korea, Thailand and mainland China (33). Although the country of origin and transmission routes of this virus have not been elucidated, available sequence data suggests that PDCoV identified in South Korea was introduced from the USA (33) indicating its international spread. The ability of porcine CoVs to emerge and re-emerge showed in recent years together with the probable naive status of European pig population for these emerging CoVs allows us to conclude the need of monitoring programmes which allow for a prompt detection and alert to establish strict biosecurity measures which would curtail their spread between countries or continents.
In order to identify PEDV variants currently circulating in Spain, we obtained the complete sequences of the S-gene of all the isolates and compared them to a representative selection of 44 PEDV genome sequences from Europe, Asia and America available in the GenBank (Figure 3). Like in other European countries (5–7, 32), only PEDV isolates included in the INDEL 2 or G1b genogroup were identified in Spain between 2017 and 2019. This genogroup has been described as causing a less severe disease than the non-INDEL or G2 genogroup (3, 8) which could explained the limited economic impact of PEDV re-emergence in Europe as compared to its dramatic consequences in the United States or Asia (34).
Phylogenetic analysis allows us to classify recent Spanish PEDV isolates into three clusters with some geographical relationships. The INDEL 2.1 and 2.2 clusters have a limited geographical spread, with the first restricted to two isolates recovered from two farms from Catalonia in the northeast of the country and the second including a number of isolates from farms located in a single province in the northwest and belonging to the same pig-producing company. Nevertheless, the INDEL 2.3 cluster included isolates recovered throughout the country together with recombinant PEDV-SeCoV isolates described in Hungary, Slovenia, Italy, or Spain (14, 15, 29) as well as Dutch, German, and French PEDV recent isolates. These isolates harbor a recombinant fragment of ~400 nt in the 5′ end of S-gene with PEDV and SeCoV being the major and minor parents, respectively. Since S protein is a key target for PEDV neutralizing antibodies, it has been proposed that this recombination event might provide some advantages (35, 36). Our results confirm this recombinant PEDV-SeCoV variant as the most widespread in Spain between 2017 and 2019 as previously suggested by de Nova et al. (15) in a research with a limited number of recent PEDV isolates. A similar displacement of other PEDV subgroups by the PEDV-SeCoV cluster has also been reported in Italy (35).
To sum up, this research confirms that PEDV has become endemic in Spain being detected in almost 40% of the viral suspected diarrhea outbreaks between 2017 and 2019 with most of the outbreaks occurring in postweaning-growing or fattening pigs, which limits the severity of the disease. Only the PEDV INDEL 2 or G1b genogroup is circulating and the recombinant PEDV-SeCoV variant is the most widespread strain. In contrast, neither PEDV non-INDEL or G2 genogroup, nor TGEV, SeCoV, PDCoV or SeACoV were detected. The emergence of new virus or variants due to spillover or through mutation or recombination events makes monitoring of porcine enteric CoVs of outmost importance in order to prevent their spread and allow for updated diagnostic tools.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Materials.
Author Contributions
HP, HA, PR, and AC conceived and designed the experiments, analyzed the data, and wrote and revised the manuscript. HP, MG-G, and ÓM-A performed the experiments. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the program from the National Institute of Agricultural and Food Research and Technology (INIA project E-RTA2015-0003-C02-02) of Spanish Government. HP, ÓM-A, and HA were supported by Spanish Government (FPU17/00466, FPU16/03485, and BEAGAL-18-106, respectively) and MG-G by Junta de Castilla y León (LE131-18).
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.
Acknowledgments
We acknowledge the excellent technique assistance provided by Diana Molina.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.651999/full#supplementary-material
Supplementary Data 1. Comparison of S1 amino acid sequences of PEDV isolates from this study with those of reference strains (GenBank accession nos. AF353511, KJ399978, and KY019623). Amino acid sequence variations in three major epitope regions (aa 220-233, aa 770-776, and aa 784-792) in the S1 genes of 41 Spanish PEDV strains.
References
1. Saif LJ, Wang Q, Vlasova AN, Jung K, Xiao S. Coronaviruses. In: Zimmerman JJ, Karriker LA, Ramirez A, Schwartz KJ, Stevenson GW, Zhang J, editors. Diseases of Swine. Hoboken: John Wiley and Sons. (2019) p. 488−523. doi: 10.1002/9781119350927.ch31
2. Fehr AR, Perlman S. Coronaviruses: an overview of their replication and pathogenesis. Methods Mol Biol. (2015) 1282:1–23. doi: 10.1007/978-1-4939-2438-7_1
3. Carvajal A, Argüello H, Martínez-Lobo FJ, Costillas S, Miranda R, de Nova PJG, et al. Porcine epidemic diarrhoea: new insights into an old disease. Porc Heal Manag. (2015) 1:1–8. doi: 10.1186/s40813-015-0007-9
4. Stevenson GW, Hoang H, Schwartz KJ, Burrough ER, Sun D, Madson D, et al. Emergence of Porcine epidemic diarrhea virus in the United States: clinical signs, lesions, and viral genomic sequences. J Vet Diagnostic Investig. (2013) 25:649–54. doi: 10.1177/1040638713501675
5. Grasland B, Bigault L, Bernard C, Quenault H, Toulouse O, Fablet C, et al. Complete genome sequence of a porcine epidemic diarrhea S gene indel strain isolated in France in December 2014. Genome Announc. (2015) 3:2014–5. doi: 10.1128/genomeA.00535-15
6. Hanke D, Jenckel M, Petrov A, Ritzmann M, Stadler J, Akimkin V, et al. Comparison of porcine epidemic diarrhea viruses from Germany and the United States, 2014. Emerg Infect Dis. (2015) 21:493–6. doi: 10.3201/eid2103.141165
7. Theuns S, Conceição-Neto N, Christiaens I, Zeller M, Desmarets LMB, Roukaerts IDM, et al. Complete genome sequence of a porcine epidemic diarrhea virus from a novel outbreak in Belgium, January 2015. Genome Announc. (2015) 3:1–2. doi: 10.1128/genomeA.00506-15
8. Vlasova AN, Marthaler D, Wang Q, Culhane MR, Rossow KD, Rovira A, et al. Distinct characteristics and complex evolution of pedv strains, North America, May 2013-February 2014. Emerg Infect Dis. (2014) 20:1620–8. doi: 10.3201/eid2010.140491
9. Lee S, Kim Y, Lee C. Isolation and characterization of a Korean porcine epidemic diarrhea virus strain KNU-141112. Virus Res. (2015) 208:215–24. doi: 10.1016/j.virusres.2015.07.010
10. Gallien S, Andraud M, Moro A, Lediguerher G, Morin N, Gauger PC, et al. Better horizontal transmission of a US non-InDel strain compared with a French InDel strain of porcine epidemic diarrhoea virus. Transbound Emerg Dis. (2018) 65:1720–32. doi: 10.1111/tbed.12945
11. Dastjerdi A, Carr J, Ellis RJ, Steinbach F, Williamson S. Porcine epidemic diarrhea virus among farmed pigs, Ukraine. Emerg Infect Dis. (2015) 21:2235–7. doi: 10.3201/eid2112.150272
12. Akimkin V, Beer M, Blome S, Hanke D, Höper D, Jenckel M, et al. New chimeric porcine coronavirus in Swine Feces, Germany, 2012. Emerg Infect Dis. (2016) 22:1314–5. doi: 10.3201/eid2207.160179
13. Belsham GJ, Rasmussen TB, Normann P, Vaclavek P, Strandbygaard B, Bøtner A. Characterization of a novel chimeric swine enteric coronavirus from diseased pigs in Central Eastern Europe in 2016. Transbound Emerg Dis. (2016) 63:595–601. doi: 10.1111/tbed.12579
14. Boniotti MB, Papetti A, Lavazza A, Alborali G, Sozzi E, Chiapponi C, et al. Porcine epidemic diarrhea virus and discovery of a recombinant swine enteric coronavirus, Italy. Emerg Infect Dis. (2016) 22:83–7. doi: 10.3201/eid2201.150544
15. de Nova PJG, Cortey M, Díaz I, Puente H, Rubio P, Martín M, et al. A retrospective study of porcine epidemic diarrhoea virus (PEDV) reveals the presence of swine enteric coronavirus (SeCoV) since 1993 and the recent introduction of a recombinant PEDV-SeCoV in Spain. Transbound Emerg Dis. (2020) 67:2911–22. doi: 10.1111/tbed.13666
16. Woo PCY, Lau SKP, Lam CSF, Lau CCY, Tsang AKL, Lau JHN, et al. Discovery of seven novel mammalian and Avian coronaviruses in the genus Deltacoronavirus supports Bat coronaviruses as the gene source of Alphacoronavirus and Betacoronavirus and Avian coronaviruses as the gene source of Gammacoronavirus and Deltacoronavi. J Virol. (2012) 86:3995–4008. doi: 10.1128/JVI.06540-11
17. Wang L, Byrum B, Zhang Y. Detection and genetic characterization of deltacoronavirus in pigs, Ohio, USA, 2014. Emerg Infect Dis. (2014) 20:1227–30. doi: 10.3201/eid2007.140296
18. Gong L, Li J, Zhou Q, Xu Z, Chen L, Zhang Y, et al. A new bat-HKU2–like coronavirus in swine, China, 2017. Emerg Infect Dis. (2017) 23:1607–9. doi: 10.3201/eid2309.170915
19. Xu Z, Zhang Y, Gong L, Huang L, Lin Y, Qin J, et al. Isolation and characterization of a highly pathogenic strain of Porcine enteric alphacoronavirus causing watery diarrhoea and high mortality in newborn piglets. Transbound Emerg Dis. (2018) 66:119–30. doi: 10.1111/tbed.12992
20. Fu X, Fang B, Liu Y, Cai M, Jun J, Ma J, et al. Newly emerged porcine enteric alphacoronavirus in southern China: identification, origin and evolutionary history analysis. Infect Genet Evol. (2018) 62:179–87. doi: 10.1016/j.meegid.2018.04.031
21. Zhou P, Fan H, Lan T, Yang X-L, Shi WF, Zhang W, et al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. Nature. (2018) 556:255–9. doi: 10.1038/s41586-018-0010-9
22. Pan Y, Tian X, Qin P, Wang B, Zhao P, Yang Y-L, et al. Discovery of a novel swine enteric alphacoronavirus (SeACoV) in southern China. Vet Microbiol. (2017) 211:15–21. doi: 10.1016/j.vetmic.2017.09.020
23. Kim O, Choi C, Kim B, Chae C. Detection and differentiation of porcine epidemic diarrhoea virus and transmissible gastroenteritis virus in clinical samples by multiplex RT-PCR. Vet Rec. (2000) 146:637–43. doi: 10.1136/vr.146.22.637
24. Kim S, Song D, Park B. Differential detection of transmissible gastroenteritis virus and porcine epidemic diarrhea virus by duplex RT-PCR. J Vet Diagnostic Investig. (2001) 13:516–20. doi: 10.1177/104063870101300611
25. Ding G, Fu Y, Li B, Chen J, Wang J, Yin B, et al. Development of a multiplex RT-PCR for the detection of major diarrhoeal viruses in pig herds in China. Transbound Emerg Dis. (2019) 67:678–85. doi: 10.1111/tbed.13385
26. Huang YW, Dickerman AW, Piñeyro P, Li L, Fang L, Kiehne R, et al. Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the united states. MBio. (2013) 4:1–8. doi: 10.1128/mBio.00737-13
27. Chen Q, Li G, Stasko J, Thomas J, Stensland W, Pillatzki A, et al. Isolation and characterization of porcine epidemic diarrhea viruses associated with the 2013 disease outbreak among swine in the United States. J Clin Microbiol. (2014) 52:234–43. doi: 10.1128/JCM.02820-13
28. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:1547–9. doi: 10.1093/molbev/msy096
29. Valkó A, Biksi I, Cságola A, Tuboly T, Kiss K, Ursu K, et al. Porcine epidemic diarrhoea virus with a recombinant S gene detected in Hungary, 2016. Acta Vet Hung. (2017) 65:253–61. doi: 10.1556/004.2017.025
30. Valkó A, Albert E, Cságola A, Varga T, Kiss K, Farkas R, et al. Isolation and characterisation of porcine epidemic diarrhoea virus in Hungary – short communication. Acta Vet Hung. (2019) 67:307–13. doi: 10.1556/004.2019.031
31. Choudhury B, Dastjerdi A, Doyle N, Frossard J-P, Steinbach F. From the field to the lab — an European view on the global spread of PEDV. Virus Res. (2016) 226:40–9. doi: 10.1016/j.virusres.2016.09.003
32. Dortmans JCFM, Li W, van der Wolf PJ, Buter GJ, Franssen PJM, van Schaik G, et al. Porcine epidemic diarrhea virus (PEDV) introduction into a naive Dutch pig population in 2014. Vet Microbiol. (2018) 221:13–8. doi: 10.1016/j.vetmic.2018.05.014
33. Zhang J. Porcine deltacoronavirus: overview of infection dynamics, diagnostic methods, prevalence and genetic evolution. Virus Res. (2016) 226:71–84. doi: 10.1016/j.virusres.2016.05.028
34. Alvarez J, Sarradell J, Morrison R, Perez A. Impact of porcine epidemic diarrhea on performance of growing pigs. PLoS ONE. (2015) 10:e0120532. doi: 10.1371/journal.pone.0120532
35. Boniotti MB, Papetti A, Bertasio C, Giacomini E, Lazzaro M, Cerioli M, et al. Porcine epidemic diarrhoea virus in Italy: disease spread and the role of transportation. Transbound Emerg Dis. (2018) 65:1935–42. doi: 10.1111/tbed.12974
Keywords: swine coronavirus, pig, PEDV, S or Spike gene, INDEL strain
Citation: Puente H, Argüello H, Mencía-Ares Ó, Gómez-García M, Rubio P and Carvajal A (2021) Detection and Genetic Diversity of Porcine Coronavirus Involved in Diarrhea Outbreaks in Spain. Front. Vet. Sci. 8:651999. doi: 10.3389/fvets.2021.651999
Received: 11 January 2021; Accepted: 05 February 2021;
Published: 25 February 2021.
Edited by:
Valentina Stefanetti, University of Perugia, ItalyCopyright © 2021 Puente, Argüello, Mencía-Ares, Gómez-García, Rubio and Carvajal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Héctor Puente, hpuef@unileon.es