Porcine Deltacoronavirus-Related Viruses in House Sparrows

1. Introduction

Porcine deltacoronavirus (PDCoV) belongs to the Deltacoronavirus genus within the Orthocoronavirinae subfamily of Coronaviridae family (https://ictv.global/taxonomy. Last accessed on July 17, 2025). The viral genome is approximately 25.4 kb and encodes 14-16 nonstructural proteins (nsps), four structural proteins [spike (S), envelope (E), membrane (M), and nucleocapsid (N) proteins] and accessory proteins, like other coronaviruses (CoVs) (16, 30). It was first reported in 2012 that PDCoV was detected from the fecal samples collected from clinically normal pigs in Asia, but its etiological role was not identified until 2014 when it caused gastroenteritis outbreaks in pigs in the United States (U.S.) (12, 27, 30). To date, it has been detected in the U.S., Canada, and many Asian countries (16). Coronaviruses (CoVs) are positive-sense RNA viruses and infect various mammalian and avian species, causing enteric, respiratory, or systemic diseases of variable severity (6). The Coronaviridae family contains four genera: Alphacoronavirus, Betacoronavirus, Gammacoronavirus and Deltacoronavirus. Recently, a new CoV was detected in the European badger and was proposed to form a new genus Epsilon (33). These viruses evolve through an accumulation of point mutations and recombination. The later may allow a virus to create new forms that can infect different species (28). AlphaCoVs and betaCoVs infect mammalian species, and gammaCoVs and deltaCoVs mainly infect birds. GammaCoVs were also detected from some whales (31). All known deltaCoVs, except PDCoV, exclusively infect birds. Besides pigs, PDCoV has been detected in Asian leopard cats and Chinese ferret badgers (5). It also experimentally infected calves, turkey and chicken poults, ducks, geese, ferrets and mice (2, 13, 22, 23, 34). Most importantly, it was also detected from three Haitian children with self-limited illness (fever, respiratory symptoms, and/or abdominal pain), making it a zoonotic pathogen (18). Molecular clock analysis suggested that PDCoV emerged in the 1990s (32) and global PDCoVs share high genetic identity (> 96%) at the genomic level (16). Sequence analyses suggest that PDCoV probably evolved from multiple recombination events among avian deltaCoVs, such as sparrow deltaCoV (SpDCoV) HKU17, munia CoV HKU13, bulbul CoV HKU11, and quail CoV HKU30 strains (16, 29).

House sparrows (Passer domesticus) are opportunistic feeders and are closely associated with swine farms due to the availability of food and nesting sites. Whether PDCoV can be transmitted between sparrows and pigs is unknown. However, such information is important to understand the ecology and evolution of PDCoVs. In this study, we captured and sampled house sparrows (Passer domesticus) near Ohio swine farms and identified PDCoV-related CoVs from the fecal samples. Our results suggest that house sparrows contribute to the complexity of PDCoV evolution.

2. Materials and Methods

2.1. Sample Collection, Screening, and Preparation

House sparrows were trapped, captured, and sampled following the approved IACUC protocol. After sample collection, the birds were released. During the winter of 2020 – 2021, we collected the fecal/cloacal swab and oropharyngeal swab samples from 198 clinically healthy house sparrows captured in mist nets or baited cage traps near eight Ohio swine farms, which tested positive for PDCoVs by the Ohio Department of Agriculture (ODA) (Table 1). Two additional house sparrow samples were collected in the Clintonville neighborhood of Columbus, Ohio. Samples were diluted tenfold with virus transport medium (VTM) [HBSS (Invitrogen, Carlsbad, CA, USA) supplemented with 2% of fetal bovine serum (FBS, HyClone, Logan, UT, USA) and 1% of Anti-Anti (Invitrogen)], vortexed, and then centrifuged at 2000 × g at 4 °C for 30 min. Subsequently, aliquots of the supernatants were made and stored at −80 °C until further analysis.

2.2. Virus Isolation in Embryonated Chicken Eggs (ECEs) and Four Cell Lines

We used 4- to 10-day-old specific pathogen-free (SPF) ECEs from The Ohio State University flock for virus isolation. The sparrow fecal suspension was filtered through 0.22 μm-pore size filters (Millipore, Chicago, IL, USA). The SF80 fecal suspension was also pretreated with a 1:100 dilution of the hyperimmune serum raised against turkey CoV (TCoV) (a gift from Dr. Mo Saif, The Ohio State University) overnight at 4°C. For each ECE, 100 - 200 μL of the sparrow fecal sample was inoculated into the allantoic sac or chorioallantoic membrane (CAM), and ≥ 4 ECEs per sample per inoculation route, as previously reported (1, 3) or illustrated by Jove video (https://app.jove.com/v/20386/egg-windowing-method-opening-fertilized-chicken-egg-shell-to-expose). Phosphate-buffered saline without Mg²⁺ and Ca²⁺ [PBS(−)] (Gibco, Carlsbad, CA, USA) was used as a negative control and PDCoV FD22 strain (GenBank accession no. MZ291567.1) at 1:100 dilution was the positive control (8). Then the ECEs were incubated at 37 °C and 55-60% relative humidity for 5 - 7 days post-inoculation (dpi) or when the ECE was found dead. After putting the ECEs at 4 °C for 2 – 6 h, the allantoic fluids were harvested and centrifuged at 2000 × g for 30 min at 4 °C to remove the cells and debris.

African green monkey kidney cell line Vero E6 (ATCC No. CRL-1586) was cultured in Dulbecco's modified Eagle's medium (DMEM) in the presence of 10% of FBS (Hyclone, Logan, UT, USA) and PSG (100 U/ml penicillin, 100 μg/ml streptomycin, and 2 mM of L-glutamine). The porcine kidney epithelial cell line LLC-PK1 (ATCC No. CL-101), swine testicle (ST) fibroblast-like cell line (ATCC No. CRL1746), and chicken embryo fibroblast cell line DF-1 (ATCC No. CRL-12203) cells were cultured in minimum essential medium (MEM) plus 5% FBS and the PSG, advanced MEM plus 5% FBS and PSG, and DMEM plus 10% FBS and the PSG, respectively (1, 8). For virus isolation, the cells were maintained in the above medium minus FBS but supplemented with 5 μg/mL of trypsin. All, except FBS, the culture reagents were from Thermo Fisher Scientific (Waltham, MA, USA).

2.3. Concentration of CoVs from the Allantoic Fluids

The allantoic fluids (~ 38 or 13 mL), which tested as weak positives for PDCoV RNA by RT-qPCR (see below), were combined and clarified by centrifugation at 2,000 x g for 30 min, followed by filtration of the supernatants through 0.45 µm-pore size filters. Then, these samples were ultracentrifuged at 106,750 x g for 2 h at 4 °C with or without 20% (w/v) sucrose cushion using an ultracentrifuge (Beckman Coulter, Miami, FL). The pellet was resuspended in 200 µL of MEM for RNA extraction.

2.4. Viral RNA Extraction

Viral RNA was extracted from 50 μL of the field sparrow fecal or oropharyngeal sample suspension, the allantoic fluids and cell culture suspension using the 5× MagMAX-96 Viral Kit (Invitrogen, Carlsbad, CA, USA) and the MagMax™ Express machine (Applied Biosystems, Bedford, MA, USA), following the manufacturer’s instructions. Finally, 50 μL of RNA was obtained in an elution buffer. Trizol (Invitrogen, Carlsbad, CA, USA) was also used to extract RNA from the concentrated allantoic fluid samples by ultracentrifugation as mentioned in section 2.3.

2.5. Conventional Reverse Transcription (RT)-PCR and TaqMan Real-Time RT-PCR (RT-qPCR)

Pan-deltaCoV forward primer (5’-TGGGATTAYCCYAAGTGTGA-3’) and reverse primer (5’-CADGCDACACCRTCATCWG-3’, located in nucleotide position 13217-13236 and 13659-13641 of deltaCoV HKU30 strain (GenBank accession no. LC364345), and pan-gammaCoV forward primer (5’-TATGAWGGYGGYTGTATMCC-3’) and reverse primer (5’-GCYCTATCACACTTWGGATARTC-3’, located in nucleotide position 13806-13825 and 14230-14208 of gammaCoV TCoV-ATCC strain (GenBank accession no. EU022526), were designed based on the conserved RNA-dependent RNA polymerase (RdRp) region of avian deltaCoVs and gammaCoVs after performing sequence alignment (see below). The Qiagen OneStep RT-PCR Kit (Qiagen, Valencia, CA, USA) was used. The reaction mixture (25 μL) consisted of 5 μL of sample RNA or water (negative control) or positive control RNA (PDCoV RNA for deltaCoVs and TCoV RNA for gammaCoVs), 10 units of RNAsin (Promega), 5 μL of 5× PCR buffer, 1 μL of deoxynucleotide triphosphates (dNTPs) (10 mM each), 1 μL of enzyme mix, 1 μL each of forward and reverse primers (100 µM each), and RNase-free water. The RT was performed at 50 ℃ for 30 min followed by inactivation at 95 ℃ for 15 min and then PCR reaction with 40 cycles (94 °C for 30 s, 50 °C for 30 s, and 72 °C for 30 s), followed by a final extension step at 72 °C for 10 min. The RT-PCR products (443 bp for deltaCoVs and 425 bp for gammaCoVs) were analyzed using 1.5% agarose gel electrophoresis to confirm the proper sizes of amplicons.

The PDCoV-specific RT-qPCR was performed as previously reported using Qiagen OneStep RT-PCR kit (12), targeting the conserved M gene of PDCoVs with forward primer 5′-ATCGACCACATGGCTCCAA-3′, reverse primer 5′-CAGCTCTTGCCCATGTAGCTT-3′, and probe 6-carboxyfluorescein-CACACCAGTCGTTAAGCATGGCAAGCT-IABkFQ. All primers and probe were synthesized by Integrated DNA Technologies (https://www.idtdna.com/page/about).

2.6. Sequencing

To sequence the panDelta and panGamma RT-PCR products of representative samples, which showed strong bands by agarose electrophoresis, the RT-PCR products of the correct size were purified using the QIAquick Gel extraction kit (Qiagen) for direct Sanger sequencing at The Ohio State University Medical Center. To obtain the nearly complete genomes of SpDCoV SF3, SF80 and SF184, viral RNA was reverse transcribed to cDNA using the SuperScript™ IV cDNA synthesis kit (Invitrogen) employing a strategy combining oligo(dT) and random hexamers’ priming. Subsequently, viral sequences were amplified using conventional PCR with the primer-walking method. The cDNA was amplified with the primers targeting the DeltaCoVs’ conserved regions and virus-specific primers designed based on the obtained viral sequences. PCR was performed using the high-fidelity PrimeSTAR GXL DNA polymerase (TaKaRa, Toshima-ku, Tokyo, Japan). To obtain additional sequences of SF80, nested PCR was also performed using PDCoV-specific primers. The purified PCR products were sent for Sanger sequencing. Finally, the overlapping DNA sequences were assembled for the near complete genomes.

2.7. Sequence Analysis

The Sanger sequence data was visualized using Chromas (https://technelysium.com.au/wp/chromas/). Reads generated from both forward and reverse primers were assembled and primer sequences were removed using the Lasergene software package (v10) (DNASTAR Inc., Madison, WI). The Basic Local Alignment Search Tool (BLAST) was used to identify homologous matches at NCBI (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Sequence alignment was performed using ClustalW, and phylogenetic trees were constructed using the Maximum likelihood method with 1,000 bootstrap replicates using MEGA 12.0 (17).

2.8. Prediction of S Protein 3D Structure

For comparative structural analysis, the amino acid sequences of the S proteins of SpDCoV/SF3, SpDCoV/SF184, and PDCoV/FD22 (GenBank accession numbers XUM38528, XXQ18308, and QWX20050, respectively) were retrieved. Structure prediction was performed using AlphaFold2 (ColabFold v1.5.5, https://colab.research.google.com/github/sokrypton/ColabFold/blob/main/AlphaFold2.ipynb) (11). Each prediction employed the PDCoV S protein structure (PDB ID: 9DEZ) as a template using an unpaired-paired alignment strategy with six recycles to optimize model quality. For each sequence, five top-ranked models were generated. The top three models were optimized using the Amber force field to reduce steric clashes and improve overall geometry. Finally, the best model was chosen for monomer illustration and for overlay of SpDCoV and PDCoV S proteins. The AlphaFold outputs included per-residue confidence (pLDDT) scores and predicted TM (pTM) scores, which evaluate the accuracy of SF3, SF184 and FD22 S structure against that of PDCoV (PDB ID: 9DEZ). Using the FD22 S protein as a reference, SF3 and SF184 S proteins were individually superimposed onto FD22, respectively, to obtain the overlays using the MatchMaker in ChimeraX (version 1.10.1) (25). Two types of overall Root Mean Square Deviation (RMSD) values (all-atom RMSD and pruned-atom RMSD) were calculated to assess the similarity between the predicted S protein structures of SF3 and SF184 and that of FD22 structure, respectively. The pruned-atom RMSD was calculated based on a subset of atoms used for structural alignment, excluding less similar or flexible regions. All structures were visualized with ChimeraX. Similarly, porcine aminopeptidase N (pAPN, GenBank No. QIS62296) was modeled using AlphaFold2 with the crystal structure of pAPN (PDB: 4FKE) as template, and the modeled pAPN together with the S proteins of SF3 and SF184 were subjected to multibody docking using HADDOCK v2.4 (https://rascar.science.uu.nl/haddock2.4/) (7). We selected the corresponding active binding residues of S proteins based on the PDCoV S protein, which were the residues 318F, 320E, 322R, 357R, 396W, 397N, and 398Y (10). The SpDCoV SF3 and SF184 contain some of these residues. Then, the PRODIGY (protein binding energy prediction) webserver (https://rascar.science.uu.nl/prodigy/) was used to calculate the binding affinity, expressed in dissociation constant (Kd), at 38.6 °C, which is the physiological internal body temperature of pigs, between the S and the APN receptor.

2.9. Nucleotide Sequence Accession Numbers

The GenBank accession no. for the SpDCoVs identified in this study are listed: PV691851 (SF3) and PV935252 (SF184). The GenBank accession no. for the viruses included in the phylogenetic analyses are listed in the phylogenetic trees.

3. Results

3.1. Detection of deltaCoVs and gammaCoVs from House Sparrows

We detected 36 deltaCoV-positive (18.0%, 36/200) and 11 gammaCoV-positive (5.5%, 11/200) samples, among which four samples (SF17, SF54, SF60, and SF80) were both positive for deltaCoV and gammaCoV, from the 200 house sparrows’ fecal/cloacal swab samples by the one-step RT-PCR assays (Table 1). The first 100 house sparrow oropharyngeal swab samples were negative. Therefore, we did not test the rest of 100 oropharyngeal swab samples. Two of the 11 gammaCoV-positive and 18 of the 36 deltaCoV-positive samples’ RT-PCR products were sequenced and BLAST confirmed their identity. Both the gammaCoVs SF60G (329-nt) and SF80G (298-nt) shared 99% nucleotide identity with the gammaCoV TCoV-ATCC strain and were not further characterized in this study. Among the deltaCoVs, 17 SpDCoVs were highly similar, sharing 98-100% nucleotide identity, and genetically closely related to quail deltaCoV HKU30 strain (93% - 95% identity) in the partial RdRp region (349 – 362-nt). Interestingly, SpDCoV strain SF80 shared 100% nucleotide identity with many PDCoV strains, including the FD22 strain in our laboratory, in this partial RdRp region (362 nt). To rule out potential contamination with the positive control PDCoV, we tested again the original sample, confirming that sample SF80 contained both a TCoV-like gammaCoV and a PDCoV-like deltaCoV. The PDCoV-like virus consistently produced only a faint band on agarose gel following RT-PCR, indicating an extremely low level of deltaCoV RNA in the sparrow droppings and suggesting inefficient replication in the sparrow.

3.2. Genomic Analysis of SpDCoV Strains SF3 and SF184

The near-complete genomes of two SpDCoVs (SF3 and SF184) were obtained. We could not obtain additional sequences beyond the partial RdRp region for SF80 although ultracentrifugation of the PDCoV RNA-positive allantoic fluids was attempted to concentrate potential virus particles and RT followed by nested PCR was performed to amplify additional genomic regions. Phylogenetic analyses showed that SF3 and SF184 grouped with the US SpDCoV strains (ISU690-4, ISU690-7, ISU73347 and ISU42824), which were detected from house sparrows at swine barns/farms in Illinois or Minnesota in 2017 (4) (Figure 1). These SpDCoVs were genetically closest to quail deltaCoV HKU30. Comparison of genomic organization also showed that SpDCoVs possess four accessory proteins—NS6, NS7a, NS7b, and NS7c—similar to other avian deltacoronaviruses, whereas PDCoV contains NS6 along with co-terminal NS7 and NS7a. (Figure 2). In ORF1ab, E, M, and N trees, these SpDCoVs clustered together and were closely related to PDCoV and HKU17 strains (Supplementary Figures).

3.3. The S Proteins of SpDCoVs Have Specific INDELs Compared with PDCoVs

In the S protein, SpDCoV HKU17 only shares ~45% identity with PDCoVs. The phylogenetic tree of the S proteins (Figure 3) showed that SF184-like SpDCoVs (SF184, ISU73347, and ISU42824) are genetically closest to PDCoVs. Among known strains, SF184, ISU42834 and ISU73347 SpDCoVs share the highest sequence identity (81.5% – 82.6%) with the earliest detected PDCoV strain AH2004 (GenBank: KP757890), followed by SF3-like SpDCoVs (75.0–75.5%), quail HKU30 (72.5%), munia HKU13 (71.6%), and bulbul HKU11 (70.3%) (Supplementary Tables). Compared with point mutations, insertion-deletions (INDELs) are less frequent. Sequence alignment of the S1 subunit of PDCoVs and genetically related avian deltaCoVs revealed PDCoV-specific INDELs in the S1 region (Figure 5), including the insertions of 41NS42 and 522R, which may reflect viral adaptation from avian to mammalian hosts. Insertions of 43SN44 and 412R and a 2-nt-Del of 427RH428 in SF184-like SpDCoVs suggest that these viruses may be the potential intermediate forms in the interspecies transmission of PDCoV from birds to pigs.

PDCoV uses pAPN as the main host receptor but may also utilize other receptors, as pAPN-knockout ST cells still support limited PDCoV replication (20). The C-terminal domain (CTD) of the S1 subunit of PDCoV S protein, spanning amino acid residues 300–419, contains the receptor-binding domain (RBD) that interacts with the conserved domain II and domain IV of pAPN (10, 26). Ji et al. identified five critical residues 318F, 320E, 322R, 396W and 398Y of PDCoV, corresponding to residues 320, 322, 327, 405 and 407 in the aligned sequences, for the receptor binding, with 396W and 398Y the most essential (10) (Figure 4). Using AlphaFold2, high-quality 3D structures of the pAPN was generated with pLDDT scores of 93.8 and pTM scores of 0.899. The predicted SF3 S and pAPN binding Kd was 12 nM, which is higher than that (2.6 nM) between SF184 and pAPN, suggesting that SF184 has a higher binding affinity to pAPN than SF3.

3.4. Prediction of S Protein 3D Structures Revealed That SpDCoV SF184 Strain Is Closer to PDCoV than SF3

Using AlphaFold2, high-quality 3D structures of the S proteins of SpDCoVs SF3 (pLDDT scores of 82.0; pTM scores of 0.663) and Sf184 (pLDDT scores of 82.5; pTM scores of 0.672) strains and PDCoV FD22 strain (pLDDT scores of 84.4; pTM scores of 0.686) were predicted (Figure 5). Structural overlays revealed that the S protein of SF3 aligned with that of FD22 with the pruned-atom RMSD of 0.440 Å (952 atom pairs), and all-atom RMSD of 15.879 Å (1154 pairs); whereas SF184 had a lower pruned-atom RMSD of 0.361 Å (957 pairs) and all-atom RMSD of 7.172 Å (1157 pairs), indicating that the core folds of SF184 are closer to PDCoV FD22 than SF3.

3.5. SpDCoV Strains Were Not Isolated in ECEs or Cell Lines

We attempted to isolate SpDCoVs since the genetically related PDCoV propagated well in ECEs and several cell lines (LLC-PK1, ST, and DF-1) (1, 8, 9, 19). However, no virus was isolated in these cell lines after three passages, as tested by the PDCoV-specific RT-qPCR.

We inoculated the allantoic sac of 5-day-old ECEs to isolate the SpDCoVs from nine samples (SF3, SF9, SF80, SF88, SF94, SF104, SF126, SF135, and SF184). During the first two passages, no ECEs were found dead and the allantoic fluids harvested at 5 dpi were negative by the panDelta RT-PCR. At passage three (P3), only one ECE inoculated with SF80 died at 2 dpi. The allantoic fluid was weakly positive for deltaCoV by panDelta RT-PCR and strongly positive for gammaCoV by the panGamma RT-PCR, suggesting that the gammaCoV SF80G, rather than the deltaCoV SF80 in the original sparrow feces, replicated well in the ECEs.

For the following virus isolation using ECEs, the SF80 fecal suspension was pre-treated with hyperimmune serum against TCoV to neutralize the gammaCoV before inoculation. The ECEs were inoculated with SF80 via two routes, the allantoic sac and CAM (Table 2). At 5-7 dpi, the allantoic fluids tested weakly positive by RT-qPCR, with high Ct values (36-38, corresponding to 1-10 RNA copies/reaction calculated from the standard curve). In all ECE experiments, the negative and positive controls performed as expected: all PDCoV-inoculated ECEs died at 2-3 dpi, and their allantoic fluids tested strongly positive by RT-qPCR (Ct values 22-24, corresponding to 5 log₁₀ RNA copies/reaction). Continuous passaging of SF80 in ECEs resulted in no significant decrease in Ct values, indicating no increase in viral RNA titer. Even after concentrating PDCoV RNA–positive allantoic fluids over 30-fold or 100-fold by ultracentrifugation, no PDCoV RNA was detected. These findings suggest that SF80 likely did not complete a full replication cycle and that no infectious viral particles were produced in the ECEs.

4. Discussion

Coronaviruses are a group of viruses that often cross species barriers and cause severe diseases and epidemics/pandemics in humans and livestock, such as severe acute respiratory syndrome (SARS)-CoV-1 and SARS-CoV-2, and porcine epidemic diarrhea virus (PEDV) (14, 21). The first PDCoV outbreak was reported in Ohio, U.S.A (27). Because house sparrows often visit farms and SpDCoVs closely related to PDCoVs were detected on farms (4), we surveyed for PDCoVs from fecal/cloacal swab and oropharyngeal swab samples collected from house sparrows near swine farms across various locations in Ohio during the winter of 2020 - 2021. DeltaCoVs were detected more frequently than gammaCoVs in fecal samples, whereas none were detected in the oropharyngeal swabs. These findings suggest that both gammCoVs and deltaCoVs primarily shed in the feces or droppings of house sparrows. All sampled sparrows appeared clinically normal. These results suggest that sparrows act as a natural reservoir for gammaCoVs and deltaCoVs, like other wild birds (29). Although SF80 shared 100% nucleotide identity with PDCoV in the partial RdRp region (362-nt), it cannot be amplified by nested RT-PCR using PDCoV-specific primers. In general, 3-7-day-old ECEs were very sensitive to PDCoV and the virus replicated efficiently in them (1, 9). Therefore, the inability of SF80 to replicate in ECEs indicate that this SpDCoV is probably quite different from PDCoV in the rest of genomic regions or the virus load in the original sample was extremely low.

We further sequenced the near-complete genomes of two representative SpDCoV strains (SF3 and SF184) that are genetically closest to the four ISU strains detected on swine barns/farms in Illinois and Minnesota (4). These results indicate that such SpDCoVs are widespread in house sparrows in the Midwest of the U.S. The forward primer and probe sequences used in our RT-qPCR assay match 100% to the M gene of SF3 and SF184, and the reverse primer had 3-nt and 2-nt mismatches to SF3 and SF184 sequences, respectively. This assay detected SF3 and SF184. So, this assay cannot differentiate PDCoV and SpDCoV when used to test environmental samples.

The identity of the CoV donor of the PDCoV S gene remains unknown. Recombinant PDCoVs carrying the S protein of SpDCoV HKU17 or the RBD of ISU73347 lose efficient replication in pigs. (24). In the S1 subunit, PDCoVs have evolved specific INDELs. These findings indicate that PDCoVs did not arise directly from any currently known SpDCoVs. Continued surveillance of diverse wild animal species may eventually identify the ancestral donor of the PDCoV S gene.

Although the disease severity of PDCoV is less than that of PEDV in pigs (15), the uniqueness of PDCoV is that it can infect both avian and multiple mammalian species (16). Phylogenetic analyses suggest that PDCoV crossed the species barrier from birds to mammals, including pigs and small carnivores such as Asian leopard cats (ALCs)/Chinese ferret badgers (CFBs) (5). Current evidence indicates that transmission likely occurred from sparrows/intermediate mammalian hosts to pigs, and subsequently from pigs to ALCs/CFBs, rather than direct spillover from sparrows/intermediate hosts to all three species independently. The PDCoV S protein contains a 125-amino-acid extension at the N-terminal domain (NTD) compared with the S proteins of the viruses detected in ALCs and CFBs from a wet market in Guangxi in 2006 (Figure 4). Aside from this NTD region, the remainder of the S protein (residues 129-1160 of PDCoV HKU15 and residues 4-1035 of ALC/CFB PDCoV S proteins) is 99% (1028/1032) identical, suggesting minimal host adaptation to ALCs/CFBs. Because pig organs were used to feed ALCs/CFBs at wet markets and molecular clock analyses estimated that PDCoV emerged around the 1990s (16). It is more plausible that PDCoV was transmitted from pigs to ALCs/CFBs. In these hosts, PDCoV likely replicated inefficiently and failed to establish sustained infection, as supported by the low viral RNA levels detected and the inability to recover complete viral genomes from the samples. Moreover, the S proteins of ALC and CFB PDCoVs share 100% identity, suggesting a single transmission event between these two species—likely during trade or co-housing in wet markets. Notably, PDCoV has only been detected in those market-traded ALCs/CFBs and not in wild populations (30), and there have been no new reports of PDCoV in these species since then.

Limitations of the current study: To detect and identify deltaCoVs and gammaCoVs, we designed new primers targeting a 400-450-nt region of the RdRp gene, based on sequences from dozens of gammaCoVs and deltaCoVs, respectively. This amplicon size allowed us to obtain sufficient sequence data for preliminary virus identification; however, they may be less sensitive compared with RT-PCR assays that use primers targeting shorter regions. As alphaCoVs and betaCoVs have not been reported in avian species, they were not targeted in this study. Recently, Drs. Saif and Vlasova laboratories at The Ohio State University, supported by WHO CoViNet, developed updated pan-CoV primers for the detection of all four genera of CoVs from different types of samples. These assays have been validated by several CoViNet reference laboratories and can be utilized for future surveillance studies.

5. Conclusions

Our results suggest that PDCoV-related SpDCoVs are widespread in the U.S. Midwest, and sparrows may play a significant role in the evolution and transmission of deltacoronaviruses. PDCoV infection in pigs can range from asymptomatic cases, commonly seen in adults, to severe diarrhea, particularly in neonatal pigs. Given the large global pig population and the frequent interactions between pigs and house sparrows—which can carry and transmit viruses to other avian species—the potential for PDCoV transmission between avian and mammalian hosts raises significant public health concerns. Continuous monitoring of this emerging zoonotic PDCoV in diverse wild animals, livestock, and humans is essential to mitigate the risk of future coronavirus pandemics.