Abstract
Swine acute diarrhea syndrome (SADS) is a highly contagious infectious disease characterized by acute vomiting and watery diarrhea in neonatal piglets. The causative agent for SADS is the swine acute diarrhea syndrome coronavirus (SADS-CoV), an alphacoronavirus in the family Coronaviridae. Currently, SADS-CoV was identified only in Guangdong and Fujian provinces of China, not in any other regions or countries in the world. To explore the genetic diversity of SADS-CoV isolates, herein we comparatively analyzed 44 full-length genomes of viruses isolated in Guangdong and Fujian provinces during 2017–2019. The spike glycoprotein gene of SADS-CoV strain CH/FJWT/2018 isolated in Fujian province is distinct from that of other viral isolates in either spike glycoprotein gene-based phylogenetic analysis or whole genome-based gene similarity analysis. Moreover, at least 7 predicted linear B cell epitopes in the spike glycoprotein of CH/FJWT/2018 would be affected by amino acid variations when compared with a representative virus isolated in Guangdong province. The spike glycoprotein of coronaviruses determines viral host range and tissue tropism during virus infection via specific interactions with the cellular receptor and also plays critical roles in eliciting the production of neutralizing antibodies. Since SADS-CoVs have a broad cell tropism, the results in this report further emphasize that the spike glycoprotein gene is a pivotal target in the surveillance of SADS-CoV.
Introduction
The swine acute diarrhea syndrome (SADS) was first reported in Guangdong province, China in 2017 (–). It is a highly contagious infectious disease that is characterized by acute vomiting and watery diarrhea in neonatal piglets, which leads to high mortality and significant economic losses for the local pork industry. The causative agent for SADS is the swine acute diarrhea syndrome coronavirus (SADS-CoV), which is also named swine enteric alphacoronavirus (SeACoV) () or porcine enteric alphacoronavirus (PEAV) (). The genome of SADS-CoV is closely related to that of the species Rhinolophus bat coronavirus HKU2, both of which are members of the genus Alphacoronavirus within the subfamily Orthocoronavirinae of family Coronaviridae of suborder Cornidovirineae of order Nidovirales. Rhinolophus bat coronavirus HKU2 was first reported from Guangdong province and Hong Kong in 2004 and 2006, respectively, in China (, ). The full-length genome of SADS-CoV is about 27 kb in length with a 5′ cap and a 3′ polyadenylated tail, and contains a 5′ untranslated region (UTR) followed by nine open reading frames (ORFs) and a 3′ UTR (). ORF1a and ORF1b encode polyprotein 1a (pp1a) and polyprotein 1b (pp1b), respectively, and occupy 5′ two-thirds of the genome. During translation, pp1a and pp1b will be processed into 16 non-structural proteins (Nsp1-16) responsible for the viral RNA synthesis (, ). The remaining part of virus genome contains ORFs that encode four structural proteins including spike (S) glycoprotein, envelope (E) protein, membrane (M) protein, and nucleocapsid (N) protein, one bigger accessory protein NS3 and two smaller accessory proteins (NS7a and NS7b) (, ). The M and E proteins of coronaviruses are required for the viral particle assembly (–). The N protein plays a critical role in the packaging of viral RNA genomes into virus particles via interactions with each other and the viral genomic RNA ().
The S glycoprotein of coronaviruses is responsible for the entry of virus into host cells by interacting with the cellular receptor. This protein can also elicit the production of neutralizing antibodies in the natural host (10, 11). Thus, the S glycoprotein serves as a primary target in the development of effective vaccines against coronaviruses and a major determinator for the efficacy of vaccines. The monomer S glycoprotein will be cleaved into S1 and S2 subunits during virus entry into host cells. S1 subunit is responsible for cell attachment and receptor binding, whereas S2 subunit mediates the membrane fusion between virus and cell host (10, 11). Thus far, the cellular receptor for S glycoprotein of SADS-CoV is still not identified. The cryo-EM structures of S glycoprotein trimers have been determined for both SADS-CoV and HKU2 (12). The overall structure of S glycoprotein trimer adopts a mushroom-like shape (~150 Å in height and ~115 Å in width), in which the β-sheets of the S1 subunit form a cap, the α-helices of the S2 subunit form a central stalk, and twisted β-sheets and loops of the S2 subunit form a root (12). The S1 subunit of monomer S glycoprotein can be further divided into several subdomains including NTD (N-terminal domain), CTD (C-terminal domain), SD1, and SD2 (12).
The first remarkable outbreak of SADS occurred in newborn piglets in February 2017 in commercial farms around the foot of a mountain located in Guangdong province of southern China, which killed 24,693 piglets (–). A retrospective investigation of SADS-CoV infection in pig farms showed that SADS-CoV had emerged in Guangdong province as early as August 2016 (). In February 2019, SADS-CoV infection re-emerged in pig herds in Guangdong, in which about 2000 pigs were killed at other pig farm near the original farms of the first outbreak (13). In 2018, a new SADS-CoV strain (CH/FJWT/2018) was identified in pig farms in Fujian province neighboring to the Guangdong province (14). In this study, we sought to explore the genomic characteristics in order to understand the genetic relatedness and genomic diversity of SADS-CoV field isolates in China.
Materials and Methods
Dataset, Phylogenetic Tree Construction, and Genomic Similarity Analysis
A total of 44 full-length genomic sequences of SADS-CoVs that were isolated in China from 2017 to 2019 were retrieved from GenBank. The sequences were aligned with Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/) (15). The phylogenetic trees were constructed using the neighbor-joining method in the MEGA-X software (16, 17) based on the whole genomic sequences or the full-length sequences of ORFs encoding pp1a and pp1b, or S, E, M, or N proteins. The viruses in this report were identified by their GenBank accession number, viral name, country, and year of collection (in a format GenBank accession number-virus name/country-year of collection). The phylogenetic inference was tested with the bootstrap method with 1,000 replications. Bootstrap values >50% were indicated. The genomic similarity plot was generated by using Simplot ver. 3.5.1 (18) to compare the whole genomes of SADS-CoV isolates.
Prediction of Linear B Cell Epitopes
The linear B cell epitopes in the S glycoprotein of SADS-CoV were predicted using BepiPred-2.0 server, which predicts B-cell epitopes using a Random Forest algorithm trained on epitopes and non-epitope amino acids determined from crystal structures (19). The residues with scores above the threshold value (set at 0.5) were predicted to be part of an epitope and shown in yellow on the graph (Y-axes depicts residue scores and X-axes indicates amino acid residue positions in the sequence). Only conserved fragments of more than 5 amino acid residues that were predicted as potential epitopes by BepiPred-2.0 were taken into consideration in this study.
Three-Dimensional Structure of Partial S Glycoprotein of SADS-CoV
Tertiary structure modeling of part of SADS-CoV S glycoprotein containing S1 subunit and the N-terminus of S2 subunit (aa 1-692 relative to the S glycoprotein of CH/FJWT/2018) was carried out using I-TASER (Iterative Threading ASSEmbly Refinement) server (20–22). I-TASER is a hierarchical approach to protein structure prediction and structure-based function annotation, which first identifies structural templates from the RCSB protein data bank (PDB) by multiple threading approach LOMETS, with full-length atomic models constructed by iterative template-based fragment assembly simulations.
Results
Comparative Analysis of Genomes of SADS-CoVs Isolated in China
Thus far, the SADS-CoV was identified only in Guangdong and Fujian provinces in China, not in any other region or countries in the world. The identified sequences in pigs are closely related to the bat HKU 2-like coronavirus sub-lineage (). To comparatively analyze genetic characteristics of SADS-CoVs, we performed a phylogenetic analysis based on the full-length genome sequence of viral strains isolated in China. The phylogenetic analysis was also done based on the complete sequence of individual ORFs encoding pp1a and pp1b, S, E, M, or N proteins. Results show that the virus CH/FJWT/2018 isolated in Fujian province (GenBank accession number MH615810.1) represents a lineage distinct from other viruses isolated in China in both whole genome-based (Figure 1A) and complete S glycoprotein ORF-based (Figure 1B) phylogenetic trees. However, CH/FJWT/2018 virus belongs to the same genogroup as most of viruses isolated in Guangdong province in the phylogenetic trees based on the complete ORFs encoding pp1a and pp1b (Figure 1C), E (Figure 1D), M (Figure 1E), or N (Figure 1F) proteins.
Figure 1
To further analyze genetic characteristics of SADS-CoV strains, the genomic similarity plot was generated using Simplot ver 3.5.1 (18) for all the non-identical full-length genomic sequences. As shown in Figure 2, the S glycoprotein ORF of CH/FJWT/2018 at nt 20,000–25,000 of viral genome shows the lowest similarity between all the viruses included. In addition, the second lowest similarity was found for the NS3 ORF of CH/FJWT/2018. The results demonstrate that S glycoprotein ORF may distinguish CH/FJWT/2018 from all other viruses isolated in China.
Figure 2
Linear B Cell Epitopes in S Glycoprotein
SADS-CoV GDWT-P18 strain (GenBank accession number: MK994935.1) was obtained by propagating in Vero cells a virulent virus CN/GDWT/2017 (GenBank accession number: MG557844) isolated in Guangdong province (23). To analyze potential antigenicity of S glycoprotein of CH/FJWT/2018, the linear B cell epitopes in the S glycoprotein were predicted using BepiPred-2.0 server (19). As shown in Figure 3, the S1 subunit shows the most obvious differences in the distribution pattern of predicted linear B cell epitopes between CH/FJWT/2018 and GDWT-P18. Furthermore, at least 7 epitopes would be affected by amino acid variations in the S glycoprotein, five of which are located in the S1 subunit and two in the N-terminus of S2 subunit (Figures 4A–C).
Figure 3
Figure 4
Discussion
SADS-CoV CH/FJWT/2018 was isolated from fecal and small intestinal samples of diarrheal piglets from pig farms in Fujian province, China, which represents the only virus identified outside of Guangdong province. The results in this report show that the S glycoprotein gene of CH/FJWT/2018 is distinct from other SADS-CoV isolates while the genome of those viruses shares over 99% nucleotide identities (14).
The S glycoprotein of coronaviruses mediates virus entry into host cells by binding host receptor (11, 24). Therefore, this protein determines the viral host range and tissue tropism. Several cellular proteins have been identified to be the receptors for coronaviruses. For example, angiotensin-converting enzyme 2 (ACE2) serves as the receptor for severe acute respiratory syndrome coronavirus (SARS-CoV or SARS-CoV-1), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and human coronavirus NL63 (HCoV-NL63) (25–27); dipeptidyl peptidase IV (DPP4) for Middle East respiratory syndrome coronavirus (MERS-CoV) (28, 29); aminopeptidase N (APN) for transmissible gastroenteritis virus (TGEV), porcine respiratory coronavirus (PRCV), and human coronavirus 229E (HCoV-229E) (30–32). For SADS-CoV, the specific receptor is still unknown and none of the known coronavirus protein receptors aforementioned are essential for the entry of SADS-CoV into the host cells (, 33). SADS-CoV has been reported to have a broad cell tropism (33). Luo et al. further confirmed this observation as SADS-CoV isolate CN/GDWT/2017 (GenBank accession number MG557844) can efficiently replicate in 20 cell lines derived from various tissues of bat, swine, human and other animal species (34). The broad cell tropism of SADS-CoV highlights the potential of cross-species infection risk. S glycoprotein of CH/FJWT/2018 shares the lowest similarity with that of all other viral isolates (Figure 2). Moreover, 7 predicted B cell epitopes in S glycoprotein of CH/FJWT/2018 have been found to contain amino acid mutations when compared with that of other viral isolates including GDWT-P18 (Figure 3), suggesting that a significant antigenic drift may occur. Collectively, these findings emphasize that S glycoprotein gene of SADS-CoV would be an informative target either in the surveillance of SADS-CoV or in the development of preventive strategies for disease control.
The finding in this report is also reminiscent of another alphacoronavirus, feline coronavirus (FCoV) that is the causative agent for feline infectious peritonitis in wild and domestic cats (35, 36). According to the antigenic and genetic differences in the S glycoprotein, FCoV exists in two distinct types, type I and type II (37–39). Feline aminopeptidase N (fAPN) is a membrane glycoprotein with metalloproteinase activity and expressed in a variety of cells (40). It serves as a cellular receptor for the binding of S glycoprotein of FCoV type II, but not for type I, during virus entry into host cells (41–43). Therefore, the primary receptor for type I FCoV remains unknown. It has yet to be explored whether the differences in the S glycoprotein of SADS-CoVs reported in this study would affect the usage of cellular receptor or the pathogenesis of disease.
Statements
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 at: https://www.ncbi.nlm.nih.gov/, MH615810.1.
Author contributions
P-HW, Y-YG, and LX: conceptualization. P-HW, Y-YG, R-ZS, Y-QL, FG, and Y-QP: data analysis. P-HW, Y-YG, R-ZS, and LX: visualization and writing. P-HW, Y-QP, Y-YG, and LX: manuscript revision. All authors contributed to the article and approved the submitted version.
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.
References
1.
PanYTianXQinPWangBZhaoPYangYLet al. Discovery of a novel swine enteric alphacoronavirus (SeACoV) in southern China. Vet Microbiol. (2017) 211:15–21. 10.1016/j.vetmic.2017.09.020
2.
GongLLiJZhouQXuZChenLZhangYet al. A new Bat-HKU2-like coronavirus in Swine, China, 2017. Emerg Infect Dis. (2017) 23:1607–9. 10.3201/eid2309.170915
3.
ZhouPFanHLanTYangXLShiWFZhangWet al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. Nature. (2018) 556:255–8. 10.1038/s41586-018-0010-9
4.
LauSKWooPCLiKSHuangYWangMLamCSet al. Complete genome sequence of bat coronavirus HKU2 from Chinese horseshoe bats revealed a much smaller spike gene with a different evolutionary lineage from the rest of the genome. Virology. (2007) 367:428–39. 10.1016/j.virol.2007.06.009
5.
YangYLLiangQZXuSYMazingEXuGHPengLet al. Characterization of a novel bat-HKU2-like swine enteric alphacoronavirus (SeACoV) infection in cultured cells and development of a SeACoV infectious clone. Virology. (2019) 536:110–8. 10.1016/j.virol.2019.08.006
6.
NeumanBWKissGKundingAHBhellaDBakshMFConnellySet al. A structural analysis of M protein in coronavirus assembly and morphology. J Struct Biol. (2011) 174:11–22. 10.1016/j.jsb.2010.11.021
7.
RuchTRMachamerCE. The coronavirus E protein: assembly and beyond. Viruses. (2012) 4:363–82. 10.3390/v4030363
8.
SchoemanDFieldingBC. Coronavirus envelope protein: current knowledge. Virol J. (2019) 16:69. 10.1186/s12985-019-1182-0
9.
McBrideRvanZyl MFieldingBC. The coronavirus nucleocapsid is a multifunctional protein. Viruses. (2014) 6:2991–3018. 10.3390/v6082991
10.
LiF. Receptor recognition mechanisms of coronaviruses: a decade of structural studies. J Virol. (2015) 89:1954–64. 10.1128/JVI.02615-14
11.
LiF. Structure, function, and evolution of coronavirus spike proteins. Annu Rev Virol. (2016) 3:237–61. 10.1146/annurev-virology-110615-042301
12.
YuJQiaoSGuoRWangX. Cryo-EM structures of HKU2 and SADS-CoV spike glycoproteins provide insights into coronavirus evolution. Nat Commun. (2020) 11:3070. 10.1038/s41467-020-16876-4
13.
ZhouLLiQNSuJNChenGHWuZXLuoYet al. The re-emerging of SADS-CoV infection in pig herds in Southern China. Transbound Emerg Dis. (2019) 66:2180–3. 10.1111/tbed.13270
14.
LiKLiHBiZGuJGongWLuoSet al. Complete genome sequence of a novel swine acute diarrhea syndrome coronavirus, CH/FJWT/2018, isolated in Fujian, China, in 2018. Microbiol Resour Announc. (2018) 7:e01259–18. 10.1128/MRA.01259-18
15.
SieversFWilmADineenDGibsonTJKarplusKLiWet al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. (2011) 7:539. 10.1038/msb.2011.75
16.
KumarSStecherGLiMKnyazCTamuraK. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:1547–9. 10.1093/molbev/msy096
17.
TamuraKNeiMKumarS. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA. (2004) 101:11030–5. 10.1073/pnas.0404206101
18.
LoleKSBollingerRCParanjapeRSGadkariDKulkarniSSNovakNGet al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. (1999) 73:152–60. 10.1128/JVI.73.1.152-160.1999
19.
JespersenMCPetersBNielsenMMarcatiliP. BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res. (2017) 45:W24–9. 10.1093/nar/gkx346
20.
YangJYanRRoyAXuDPoissonJZhangY. The I-TASSER Suite: protein structure and function prediction. Nat Methods. (2015) 12:7–8. 10.1038/nmeth.3213
21.
RoyAKucukuralAZhangY. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. (2010) 5:725–38. 10.1038/nprot.2010.5
22.
YangJZhangY. I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res. (2015) 43:W174–81. 10.1093/nar/gkv342
23.
SunYChengJLuoYYanXLWuZXHeLLet al. Attenuation of a virulent swine acute diarrhea syndrome coronavirus strain via cell culture passage. Virology. (2019) 538:61–70. 10.1016/j.virol.2019.09.009
24.
TortoriciMAVeeslerD. Structural insights into coronavirus entry. Adv Virus Res. (2019) 105:93–116. 10.1016/bs.aivir.2019.08.002
25.
WuKLiWPengGLiF. Crystal structure of NL63 respiratory coronavirus receptor-binding domain complexed with its human receptor. Proc Natl Acad Sci USA. (2009) 106:19970–4. 10.1073/pnas.0908837106
26.
LiFLiWFarzanMHarrisonSC. Structure of SARS coronavirus spike receptor-binding domain complexed with receptor. Science. (2005) 309:1864–8. 10.1126/science.1116480
27.
LanJGeJYuJShanSZhouHFanSet al. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature. (2020) 581:215–20. 10.1038/s41586-020-2180-5
28.
LuGHuYWangQQiJGaoFLiYet al. Molecular basis of binding between novel human coronavirus MERS-CoV and its receptor CD26. Nature. (2013) 500:227–31. 10.1038/nature12328
29.
WangNShiXJiangLZhangSWangDTongPet al. Structure of MERS-CoV spike receptor-binding domain complexed with human receptor DPP4. Cell Res. (2013) 23:986–93. 10.1038/cr.2013.92
30.
RegueraJSantiagoCMudgalGOrdonoDEnjuanesLCasasnovasJM. Structural bases of coronavirus attachment to host aminopeptidase N and its inhibition by neutralizing antibodies. PLoS Pathog. (2012) 8:e1002859. 10.1371/journal.ppat.1002859
31.
LiZTomlinsonACWongAHZhouDDesforgesMTalbotPJet al. The human coronavirus HCoV-229E S-protein structure and receptor binding. Elife. (2019) 8:51230. 10.7554/eLife.51230
32.
WongAHMTomlinsonACAZhouDSatkunarajahMChenKSharonCet al. Receptor-binding loops in alphacoronavirus adaptation and evolution. Nat Commun. (2017) 8:1735. 10.1038/s41467-017-01706-x
33.
YangYLQinPWangBLiuYXuGHPengLet al. Broad cross-species infection of cultured cells by bat HKU2-related swine acute diarrhea syndrome coronavirus and identification of its replication in murine dendritic cells in vivo highlight its potential for diverse interspecies transmission. J Virol. (2019) 93:e01448–19. 10.1128/JVI.01448-19
34.
LuoYChenYGengRLiBChenJZhaoKet al. Broad cell tropism of SADS-COV in vitro implies its potential cross-species infection risk. Virol Sin. (2020) 17:1–5. 10.1007/s12250-020-00321-3
35.
WolfeLGGriesemerRA. Feline infectious peritonitis. Pathol Vet. (1966) 3:255–70. 10.1177/030098586600300309
36.
PedersenNC. A review of feline infectious peritonitis virus infection: 1963-2008. J Feline Med Surg. (2009) 11:225–58. 10.1016/j.jfms.2008.09.008
37.
HohdatsuTOkadaSKoyamaH. Characterization of monoclonal antibodies against feline infectious peritonitis virus type II and antigenic relationship between feline, porcine, canine coronaviruses. Arch Virol. (1991) 117:85–95. 10.1007/BF01310494
38.
PedersenNCBlackJWBoyleJFEvermannJFMcKeirnanAJOttRL. Pathogenic differences between various feline coronavirus isolates. Adv Exp Med Biol. (1984) 173:365–80. 10.1007/978-1-4615-9373-7_36
39.
JaimesJAMilletJKStoutAEAndreNMWhittakerGR. A tale of two viruses: the distinct spike glycoproteins of feline coronaviruses. Viruses. (2020) 12:83. 10.3390/v12010083
40.
LookATAshmunRAShapiroLHPeiperSC. Human myeloid plasma membrane glycoprotein CD13 (gp150) is identical to aminopeptidase N. J Clin Invest. (1989) 83:1299–307. 10.1172/JCI114015
41.
HohdatsuTIzumiyaYYokoyamaYKidaKKoyamaH. Differences in virus receptor for type I and type II feline infectious peritonitis virus. Arch Virol. (1998) 143:839–50. 10.1007/s007050050336
42.
DyeCTempertonNSiddellSG. Type I feline coronavirus spike glycoprotein fails to recognize aminopeptidase N as a functional receptor on feline cell lines. J Gen Virol. (2007) 88:1753–60. 10.1099/vir.0.82666-0
43.
TresnanDBLevisRHolmesKV. Feline aminopeptidase N serves as a receptor for feline, canine, porcine, and human coronaviruses in serogroup I. J Virol. (1996) 70:8669–74. 10.1128/jvi.70.12.8669-8674.1996
Summary
Keywords
phylogenetic analysis, epitope, spike glycoprotein, virus genome, swine acute diarrhea syndrome coronavirus (SADS-CoV)
Citation
Guo Y-Y, Wang P-H, Pan Y-Q, Shi R-Z, Li Y-Q, Guo F and Xing L (2021) The Characteristics of Spike Glycoprotein Gene of Swine Acute Diarrhea Syndrome Coronavirus Strain CH/FJWT/2018 Isolated in China. Front. Vet. Sci. 8:687079. doi: 10.3389/fvets.2021.687079
Received
28 March 2021
Accepted
28 June 2021
Published
22 July 2021
Volume
8 - 2021
Edited by
Shao-Lun Zhai, Guangdong Academy of Agricultural Sciences, China
Reviewed by
Lang Gong, South China Agricultural University, China; Jingyun Ma, South China Agricultural University, China
Updates
Copyright
© 2021 Guo, Wang, Pan, Shi, Li, Guo and Xing.
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: Li Xing xingli107@gmail.com
This article was submitted to Veterinary Infectious Diseases, a section of the journal Frontiers in Veterinary Science
†These authors have contributed equally to this work and share first authorship
Disclaimer
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.