Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Nature Microbiology volume 8, pages 1820–1833 (2023)
53
Metrics details
The pathogenic and cross-species transmission potential of SARS-CoV-2-related coronaviruses (CoVs) remain poorly characterized. Here we recovered a wild-type pangolin (Pg) CoV GD strain including derivatives encoding reporter genes using reverse genetics. In primary human cells, PgCoV replicated efficiently but with reduced fitness and showed less efficient transmission via airborne route compared with SARS-CoV-2 in hamsters. PgCoV was potently inhibited by US Food and Drug Administration approved drugs, and neutralized by COVID-19 patient sera and SARS-CoV-2 therapeutic antibodies in vitro. A pan-Sarbecovirus antibody and SARS-CoV-2 S2P recombinant protein vaccine protected BALB/c mice from PgCoV infection. In K18-hACE2 mice, PgCoV infection caused severe clinical disease, but mice were protected by a SARS-CoV-2 human antibody. Efficient PgCoV replication in primary human cells and hACE2 mice, coupled with a capacity for airborne spread, highlights an emergence potential. However, low competitive fitness, pre-immune humans and the benefit of COVID-19 countermeasures should impede its ability to spread globally in human populations.
The coronavirus disease 2019 (COVID-19) pandemic is caused by severe acute respiratory disease syndrome coronavirus 2 (SARS-CoV-2), which emerged in southern China, causing >750 million infections and >6.9 million deaths. SARS-CoV-2 most probably emerged from animal reservoirs in Southeast Asia and is efficiently transmitted between humans1,2. Although an intermediate species remains elusive, Malayan pangolins (Manis javanica) represent one of several potential species3,4,5,6. Pangolins harbour diverse Sarbecoviruses that are genetically related to SARS-CoV-2, as well as strains circulating in horseshoe bats (Rhinolophus)3,4,7,8,9,10,11. While none of these strains are considered the direct progenitor of SARS-CoV-2, the pangolin (Pg) CoV strain GD, bat CoVs BANAL-52 and BANAL-103 strains encode receptor binding domains (RBD) that display remarkable homology with the SARS-CoV-2 RBD3,11. These closely related SARS-2-like CoVs lack a polybasic furin cleavage site (RRAR) at the S1/S2 junction of the SARS-CoV-2 spike protein, which is thought to promote SARS-CoV-2 spike RBD-ACE2 receptor binding12, pathogenesis13, transmission14 and emergence15. However, it is unclear whether these related strains can efficiently infect and replicate in primary human cells, undergo airborne transmission between hosts and cause severe disease in mammals.
Micro-variation in the SARS-CoV-2 RBD alters species-specific receptor usage and transmission potential in mammals, and humans transmit virus to other mammalian species such as mink, feline and Cervidae populations16,17,18. The discovery of SARS-2-like bat CoV BANAL-52 and PgCoV GD support previous hypotheses that many zoonotic CoVs are pre-determined to replicate efficiently in multiple mammalian hosts, including humans19,20. The earliest known PgCoV strains, including the GD and GX clusters, were identified in 2 batches of smuggled pangolins. For the GD strain, the same metagenomic data from a single source sample were reported in different publications3,6,21, triggering confusion about origin and authenticity. While human infections have not been reported, pseudotyped viruses encoding the PgCoV GD spike can utilize human ACE2 receptor for entry9. However, PgCoV replication and transmission efficiency, tropism and fitness in primary human respiratory epithelial cells remain uncertain21, as does its ability to replicate and transmit in small animal models. Lastly, the performance of current COVID-19 antiviral drugs, vaccines and therapeutic antibodies to protect from live virus PgCoV zoonotic events in vivo is important for global health preparedness.
The PgCoV GD genome shares 90.1% identity with SARS-CoV-2 and encodes 90.7% and 96.7% of homologous residues in the full-length spike protein and RBD, respectively (Supplementary Fig. 1a,b). To evaluate whether the PgCoV genome is viable in cell culture, a synthetic complementary (c)DNA clone was generated from the full-length PgCoV GD strain sequence22,23,24,25, including derivatives that replaced ORF7a with either the green fluorescent protein (GFP) or the nano-luciferase (nLuc) reporter genes. After assembly of the multicomponent cDNA clone, full-length viral RNA genomes were transcribed in vitro using a T7-RNA polymerase and then electroporated into Vero-E6 cells. Cytopathic effects (CPEs) were extensive by day 3–4 post infection and recombinant viruses formed visible plaques (Fig. 1a,b). To evaluate PgCoV GD transcription, northern blot analysis demonstrated full length and eight subgenomic mRNA (sgRNA), phenocopying SARS-CoV-2 transcription profiles (Fig. 1c). As expected, the molecular weights of sgRNA2 to sgRNA7 were visibly increased in PgCoV derivatives encoding the larger GFP and nLuc reporter transgenes as compared with the wild-type (WT) virus.
a, Schematic design of PgCoV GD1 infectious cDNA clone. b, Plaque formation of the three recombinant PgCoVs. Scale bar, 8 mm. c, Northern blot analysis of genomic and subgenomic mRNAs isolated from SARS-CoV-2 and PgCoV infected cells at 24 h. d, Western blot analysis of semi-purified PgCoV and SARS-CoV-2 virions cultured in Vero or Vero-furin cell lines identified as full-length (FL), S1/S2 cleaved spike protein (S) and nucleocapsid protein (N). Samples were loaded on the basis of an equal amount of the N protein; this western blot was repeated twice with the same result. e, Efficient entry of PgCoV-nLuc and SARS-CoV-2-nLuc recombinant viruses into Vero-81 and Vero-furin cells at an MOI of 2. After 1 h infection, viruses were removed and cells were treated with neutralization antibodies to minimize secondary rounds of infection. The RLU representing the nLuc expression level was measured at 12 h post infection (n = 8 replicates per group, data are mean ± s.d.) and analysed using unpaired t-test; ****P < 0.0001. f–h, Multistep growth curves of PgCoV-WT (f), SARS-CoV-2-WT (g) and SARS-CoV-2-∆PRRA (h) in Vero and Vero-furin cells.
Source data
Group 1b Sarbecoviruses may encode a furin cleavage site at the S1/S2 boundary and/or a less efficient site in S2 (refs. 26,27,28). To evaluate the expression and processing of the PgCoV spike glycoprotein which lacks a polybasic cleavage site at the S1/S2 boundary, western blot analysis of PgCoV-WT and SARS-CoV-2 WT (strain WA1) virions were compared in infected wild-type Vero-E6 or in a previously described furin-overexpressing Vero-E6 cell line (Fig. 1d)29. SARS-CoV-2 grown in wild-type Vero cells displayed both full-length and S1/S2 processed spike glycoproteins, whereas the S1/S2 cleavage products were prominent in virions acquired from Vero-furin cells. In contrast, the PgCoV glycoprotein spike was insensitive to S1/S2 cleavage, even in furin-overexpressing cell lines, although the presence of a ~140–160 kDa spike product in PgCoV infected cells may reflect limited cleavage at a second S2 furin cleavage site, reported in SARS-CoV and PgCoV strains27. Next, we compared PgCoV-nLuc and SARS-CoV-2-nLuc entry efficiency in Vero and Vero-furin cells, the latter expressing high levels of furin (Fig. 1e)29. Both viruses display enhanced entry in the Vero-furin cells compared with Vero-E6 cells, as evidenced by ~5- and ~10-fold significantly increased levels of nLuc expression in PgCoV and SARS-CoV-2 infected cultures, respectively, through 48 h.
To compare virus growth, we performed a multistep growth curve analysis in both Vero cell types (Fig. 1f,g). To investigate the role of the furin cleavage site in SARS-CoV-2 entry, we generated a recombinant SARS-CoV-2 mutant ∆PRRA lacking the S1/D2 furin cleavage site (Supplementary Fig. 1c) and compared the growth kinetics of PgCoV GD, SARS-CoV-2 and SARS-CoV-2 ∆PRRA in culture (Fig. 1h). In Vero-furin cells, all three viruses replicated to similar titres of ~7 × 107 plaque forming unit (p.f.u.) ml−1 at 24 h. At early times, SARS-CoV-2 replicated efficiently in both cell lines, although virus titre decay rates were accelerated in Vero-furin cell lines. Importantly, PgCoV GD growth characteristics most closely approximated the phenotypes seen with the furin-site knockout virus SARS-2-∆PRRA in both cell lines. Together, these data support earlier studies that showed that the S1/S2 furin cleavage facilitates SARS-CoV-2 entry and replication efficiency in vitro13, but that other factors may also contribute.
The PgCoV GD S RBD is highly homologous with SARS-CoV-2 (Supplementary Fig. 1d). To examine the receptor specificities of PgCoV GD, we used non-permissive cells stably transfected with ACE2 receptors derived from pangolin, human and many other mammalian species (Fig. 2). Using PgCoV GD GFP fluorescent viruses, we demonstrate robust virus replication in all the cell lines except non-permissive cells, although some differences in fluorescence were noted (Fig. 2a). Using western blot analyses with anti-ACE2 antibody, we further demonstrate high levels of mammalian ACE2 orthologue receptor expression from each species, but endogenous ACE2 expression was not detected in non-permissive cells (Fig. 2b).
a, GFP signals were observed in DBT9 cells expressing ACE2 orthologues or empty vector (E.V.) infected with icPgCoV-GFP at 24 h post infection. Scale bar, 100 μM. b, Western blot analysis of ACE2 orthologue or E.V. expression from 15 mammalian species in DBT9 cells. a and b were repeated twice with the same results. c, In vitro potency of RDV, EIDD and nirmatrelvir against icPgCoV in A549-hACE2 cells based on RLU of icPgCoV replication in the presence of decreasing concentrations of RDV, EIDD and nirmatrelvir in A549-hACE2 cells. Quantitation of icPgCoV-nLuc replication as measured by nLuc in technical triplicates. Shown are the IC50 values of RDV, EIDD and nirmatrelvir, n = 3 replicates per group, data are mean ± s.d.
Source data
Given the broad host range of PgCoV GD, we evaluated virus sensitivity to US Food and Drug Administration (FDA)-approved drugs such as remdesivir (RDV) and molnupiravir (EIDD), which have broad antiviral activity against many coronaviruses30,31. Using A549-hACE2 cells32, the PgCoV GD nLUC reporter virus was highly sensitive to current FDA-approved drugs such as RDV, EIDD and the oral protease inhibitor nirmatrelvir, demonstrating half-maximal inhibitory concentration (IC50) values of 0.099, 0.166 and 0.0128 μM, respectively (Fig. 2c). Together, these data not only demonstrate the broad host range potential of PgCoV GD but also revealed its high sensitivity to FDA-approved small-molecule inhibitors that target conserved nsp12 polymerase and nsp5 protease activities.
To evaluate the ability of PgCoV to replicate in human respiratory tissues, primary human cells from large airway epithelia (LAE) and nasal epithelia (HNE) from different donors (n = 2 per cell type) were infected with either PgCoV-WT or SARS-CoV-2 WA1 strain at an MOI of 0.5 (Fig. 3). Although donor variation influenced replication efficiency, PgCoV exhibited slightly increased or equivalent growth, as compared with SARS-CoV-2. In HNE, PgCoV replicated to titres of 6.6 × 106 and 1.62 × 107 p.f.u. ml−1 at 5 d post infection, whereas SARS-CoV-2 peaked at ~2-fold lower titres of 3.16 × 106 and 6.6 × 106 p.f.u. ml−1, respectively (Fig. 3a). In one LAE culture, PgCoV and SARS-CoV-2 titres plateaued between ~104 and 105 p.f.u. ml−1, while another patient code did not efficiently support the growth of either virus (Fig. 3b). Similar variation in Sarbecovirus growth across different patient codes has been reported previously22. Immunofluorescence staining of PgCoV antigens demonstrated that ciliated nasal and bronchial epithelial cells were the primary targets of viral infection (Fig. 3d), recapitulating SARS-CoV-2 tropism for ciliated cells in these tissues22.
a–c, Comparison of PgCoV-WT and SARS-CoV-2 (WA1) growth in primary HNE cells (a), LAE cells (b) and AT2 organoids (c). The growth curve values were measured in triplicate for each timepoint, data are mean ± s.d. d, Immunofluorescence co-staining of PgCoV-WT antigens in representative HNE and LAE cultures at 72 h post infection. Scale bar, 50 μm. e, Live image of PgCoV-GFP virus-infected AT2 organoids at 48 h (left; scale bar, 40 μm) and immunofluorescence co-staining (right; scale bar, 25 μm) of PgCoV-WT (green) with SFTPC (red) and AGER (white) on infected AT2 organoids fixed at 72 h. f, Percentage of NTD coding sequences in RNA samples collected from PgCoV and SARS-CoV-2 competition experiments, n = 3 samples per timepoint, data are mean ± s.d.
Source data
SARS-CoV-2 primarily infects alveolar type II (AT2) pneumocytes, which maintain alveolar homoeostasis and gas exchange, and whose loss is correlated with severe disease outcomes during COVID-19 infection22. To evaluate the ability of PgCoV replication in distal lung tissues, human AT2 organoid cultures from two donors33 were infected with PgCoV-WT or SARS-CoV-2 at an MOI of 0.5 (Fig. 3c). In the AT2 organoids from donor 5, SARS-CoV-2 replicated to 7.41 × 105 p.f.u. ml−1, a log higher than the peak PgCoV titre of 4.57 × 104 p.f.u. ml−1. In organoid patient code 6, both viruses displayed similar growth, peaking at titres of ~104 p.f.u. ml−1. The replication of PgCoV in AT2 organoids was also confirmed by immunofluorescence analysis (Fig. 3e), which demonstrated expression of viral antigens in close proximity with surfactant protein C (SFTPC)33. Together, these data indicate that PgCoV can replicate efficiently in human respiratory tissues from the proximal to the distal lung.
To compare replication fitness, competition assays were performed in LAE cultures simultaneously infected at an MOI of 1 each for PgCoV-WT and either SARS-CoV-2 strain WA1 or B.1.1.7 (UK variant). After 4 d, viral RNAs from progeny virions collected at 48 h intervals were quantified by deep sequencing of the variable spike N-terminal domain (NTD) (Fig. 3f). After 48 h, both WA1 and B.1.1.7 dominated the culture, with more than 98% of the RNA amplicons derived from these human strains. Together, these data suggest that the competitive fitness of PgCoV is lower than these early SARS-CoV-2 strains.
PgCoV efficiently uses the hamster ACE2 receptor for entry (Fig. 2), so hamsters were intranasally infected with 103 p.f.u. of PgCoV or SARS-CoV-2. By 2–3 d, PgCoV replicated to titres of ~7.0 × 104 p.f.u. g−1 and ~6.0 × 105 p.f.u. g−1 in lung and nasal turbinate tissues, respectively, similar to titres achieved by SARS-CoV-2 (Fig. 4). At 4 d post infection, all hamsters developed mild broncho-interstitial pneumonia, with pathology scores ranging from 3 of 16 to 6 of 16, depending on the animal (Fig. 4d and Supplementary Table 1). Immunohistochemistry signals were prominent in bronchiolar and alveolar epithelial cells (pneumocytes) and presumed to be macrophages. To evaluate the airborne transmission potential, 18 or 6 hamsters were infected with PgCoV GD or SARS-CoV-2, as described in earlier studies34,35,36. After 1 d, each infected hamster was housed adjacent to an uninfected recipient hamster (Fig. 4a–c), separated by a wire-mesh divider (5 cm) to acquire airborne transmission. Viral titres in nasal wash demonstrated that SARS-CoV-2 and PgCoV GD replicated to similar titres on day 2, although PgCoV appeared to clear faster at later times. While SARS-CoV-2-infected donors rapidly transmitted virus to 100% of the uninfected recipients after 3–4 d (6/6), only three hamsters in the PgCoV donor group successfully transmitted the infectious virus to an uninfected recipient (3/18) by day 6, demonstrating a significantly lower airborne transmission rate (~16.6%), as compared with SARS-CoV-2 (P < 0.001, Fisher exact test) (Fig. 4a–c).
a–c, Viral titres in nasal washes were determined from infected and exposed hamster pairs in SARS-CoV-2 (a) and 3 replicate experiments involving PgCoV groups (b and c (groups 1 and 2)). Uninfected groups of exposed animals were housed in separate cages with a 5-cm distance from the infected hamster. d.p.i., days post infection. d, PgCoV disease pathology in the lung (left) and immunohistology staining for PgCoV N protein (right). Boxes on left identify tissue regions examined under higher magnification on right.
The PgCoV GD spike protein has ~10% amino acid variation when compared with the prototype SARS-CoV-2 WA1 strain that encodes significant variation in spike NTD (Supplementary Fig. 1a). The PgCoV GD RBD also contains 6 amino acid changes (R346T, A372T, I402V, K417R, Q498H and H519N) as compared with the SARS-CoV-2 RBD, including residues that bind ACE2 (Supplementary Fig. 1b,d). To evaluate the antigenic relationships, we compared IC50 neutralization titres using SARS-CoV-2 and PgCoV GD nLuc live reporter viruses and a panel of commercially available and broad-spectrum SARS-CoV-2 antibodies (Fig. 5a and Supplementary Table 2). The 13 RBD-binding human Abs exhibit similar neutralization profiles against both CoVs, with less than a 10-fold difference in IC50 values. Under conditions where RBD class 1 and class 2 antibodies were potent neutralizers, RBD class III and S2 Ab either failed to neutralize or showed reduced potency against PgCoV GD. For example, class III antibodies, Ab COV2-2499 and COV2-2130 either exhibited a 19.4-fold reduction or failed to neutralize PgCoV, most probably associated with the R346T and Q498H mutations (Fig. 5b). Paradoxically, some RBD-specific class 1 and class 2 Abs (B38, LY-COV016, LY-COV555, ADG-2 and ADG-3) showed slightly more potent neutralization titres against PgCoV. While most commercial antibodies potently neutralized PgCoV, REGN10933 (Class 1) and 10987 (Class 3) were ~5-fold less potent. In contrast, all three NTD-binding Abs neutralized SARS-CoV-2 but not PgCoV GD, which encodes variation in the NTD Ab-binding epitopes (Fig. 5c). In addition, several COVID-19 patient convalescent sera showed similar half half-maximal inhibitory dilution (ID50) values (Fig. 5d,e), demonstrating that most, but not all COVID-19 patients, elicit a robust neutralizing response against PgCoV GD and SARS-CoV-2.
Neutralization assays were performed using SARS-CoV-2 nLUC and PgCoV-nLUC recombinant viruses and panels of neutralizing monoclonal antibodies and convalescent patient sera. a, Summary of IC50 values from a panel of SARS-CoV-2 neutralizing antibodies. b, RBD antibody binding sites (pink) impacted by PgCoV natural variation (red). c, NTD antibody binding sites (pink) impacted by PgCoV natural variation (red). Binding sites in b and c were annotated on the basis of a SARS-CoV-2 spike structure (PBD ID: 6zp7). d, COVID-19 patient sera 8-point neutralization curves against SARS-CoV-2 (WA1) and PgCoV; each sample was run in triplicate in the assay, data are mean ± s.d. e, COVID-19 ID50 neutralizing titres against SARS-CoV-2 and PgCoV. Values were analysed using paired t-test; NS, P > 0.05.
To evaluate PgCoV growth in standard laboratory mice, 10-month-old BALB/c mice were infected with 1 × 104 p.f.u. PgCoV intranasally. In agreement with murine ACE2 receptor usage (Fig. 2), PgCoV replicated efficiently in the lungs of mice, but demonstrated little mortality, weight loss and gross lung pathology lesions as compared with SARS-CoV MA10 (Fig. 6a–d). As the pan-Sarbecovirus neutralizing antibody ADG-2 exhibited the highest potency against PgCoV in vitro (Fig. 5 and Supplementary Table 2), mice were treated with PBS, ADG-2 or a dengue isotype IgG control and challenged with either PgCoV or mouse-adapted SARS-CoV-2 MA10. In controls, both viruses exhibited high viral titres in the lungs on days 2 and 4 post infection, whereas ADG-2-treated mice were protected (Fig. 6a–c). Consistent with our previous finding37, ADG-2 significantly reduced PgCoV GD and SARS-CoV-2 MA10 weight loss, lung viral titre and gross lung discoloration (Fig. 6a–d). These data indicate that the pan-Sarbecovirus neutralizing Ab ADG-2 protects mice from PgCoV replication in the lung.
N = 10 per group of BALB/c mice were treated with either ADG-2 or an isotype control antibody and then infected with PgCoV or SARS-CoV-2 MA10. PgCoV and MA10 infection (104 p.f.u. per mouse) in the presence or absence of ADG-2 neutralizing antibody in BALB/c mice, n = 5 mice per group were collected at days 2 and 4 for titre and histological analysis. a, Weight loss. b, Lung viral titres at days 2 and 4 post infection. c, Lung discoloration scores. Another cohort of BALB/c mice (n = 10 mice per group) were vaccinated with SARS-CoV-2 S2P protein along with alum adjuvant and then challenged with PgCoV GD or mouse-adapted SARS-CoV-2 MA10 (104 p.f.u. per mouse per virus), n = 5 mice per group were collected at days 2 and 4 for titre and histological analysis. d, Weight loss in PgCoV and SARS-CoV-2 MA10 challenged mice. e, Virus titres in the lung. f, Lung congestion scores of the challenged mice. We next evaluated PgCoV GD pathogenesis and DH1047 prophylaxis studies in K18-hACE2 mice, n = 15 mice per group. g, Weight loss following PgCoV GD challenge. h, PgCoV GD titres in the lung, nasal turbinates and brain (n = 5 per timepoint). i, Lung discoloration scores following PgCoV challenge. j, Percent survival following PgCoV GD challenge. k, Lung respiratory function (PenH) as measured by whole-body plethysmography. l–n, In another set of hACE2-K18 mice (n = 7 for the DH1047 group and n = 5 for the CH65 group), we demonstrated that the prophylactic administration of DH1047 protected against PgCoV weight loss (l), and reduces virus titre (m) and lung congestion (n). Data are mean ± s.d.
To evaluate whether a SARS-CoV-2 vaccine elicits cross-protection against PgCoV, we also immunized and boosted 10-week-old BALB/c mice with soluble SARS-CoV-2 Spike-2P (S2P) protein with alum adjuvant38,39 and demonstrated high neutralizing antibody titres in serum against both SARS-CoV-2 and PgCoV GD (Supplementary Fig. 2). After intranasal challenge with either PgCoV GD or SARS-CoV-2 MA10 at 3 weeks post boost (7 weeks post prime) (Fig. 6d–f), the S2P vaccine protected mice from SARS-CoV-2 MA10 and PgCoV GD. Little if any virus replication was noted in the lung, indicating that the SARS-CoV-2 S2P recombinant protein vaccine elicited robust cross-protection against both viruses.
Although wild-type mice provide insight into the pathogenic potential and mouse ACE2 receptor usage in vivo, the lack of clinical disease is inadequate to access human disease potential. To test a model more relevant to human populations, we evaluated PgCoV GD growth in the K18 human ACE2 mouse model, which has been repeatedly shown to progress to severe clinical disease and death following SARS-CoV-2 infection40. Following PgCoV GD infection, K18-hACE2-expressing mice lost weight (~25%), developed lung respiratory dysfunction, demonstrated lung discoloration and had high-titre virus replication on days 2 through 7–8 post infection (Fig. 6). Mortality rates approached 60% (Fig. 6j). To evaluate countermeasure performance, animals were then treated with hmAB DH1047, which has been shown to broadly neutralize many group 1a and 1b Sarbecoviruses and variants of concern41. Not only did prophylactic DH1047 administration protect from serious weight loss and clinical disease following PgCoV GD challenge, but it also significantly reduced lung discoloration scores and virus titres by about 3 logs on day 4 post infection (Fig. 6l–n). The PgCoV severe disease model in hACE2-K18 mice offers novel opportunities for countermeasure testing.
Coronaviruses have emerged multiple times from zoonotic reservoirs in the twenty-first century and SARS-CoV-2 mostly probably originated from strains that circulated in bats or other mammals11,42. Despite undersampling, SARS-CoV-2-like animal strains are indigenous to Southeast Asia and encode up to 96.8% identify with SARS-CoV-2 (ref. 11). One metagenomic study reported a nearly completed PgCoV genome6 on the basis of earlier work on smuggled pangolins confiscated in Guangdong province in March 201943. Another group reanalysed the same pangolin samples and published a different full-length PgCoV genome GD3,21. Meanwhile, the genomes of a distant cluster of PgCoV strains GX were identified in another batch of smuggled pangolins4,44. The GD and GX PgCoV strains have 90.4% and 85.5% identity to the SARS-CoV-2 genome, respectively. While neither of these PgCoV strains is considered the immediate ancestor of SARS-CoV-2, PgCoV strain genetic diversity is substantially higher than previously recognized yet dramatically undersampled, increasing the probability that pangolins may have participated in the SARS-CoV-2 pandemic45. Unfortunately limited virus sampling, coupled with the illegal trade of these endangered species, has hampered a full characterization of the SARS-like strain diversity that circulates in pangolins6,44,46.
The PgCoV GD RBD preserves 13/14 hACE2 interaction sites, suggesting that recombination events between an unidentified bat Sarbecovirus and PgCoV GD may have contributed to SARS-CoV-2 emergence5. However, this hypothesis is less compelling because BANAL-52 and related group 1b strains encode near-identical RBD-ACE2 interacting residues with SARS-CoV-2 (ref. 11). Using pseudotypes, PgCoV spikes are reported to have a high affinity for hACE2 entry receptors47,48,49,50, although the in vitro and in vivo infection characteristics of the PgCoV GD strain remain unclear. In this study, we recovered a recombinant wild-type PgCoV GD strain and characterized its ACE2 receptor breadth, host range, tropism, pathogenesis in two host species, transmissibility and antigenicity. Supporting structural studies47, the PgCoV GD strain and other closely related clade 1b SARS-CoV-2 zoonotic strains11 may efficiently engage a large number of mammalian ACE2 receptors for docking and entry, replicate efficiently in multiple mammalian host species, and cause clinical disease and mortality in k18-hACE2-expressing mice, but not in standard laboratory mice. Importantly, PgCoV GD replicates efficiently in primary human epithelial cells from the proximal to distal human lung and undergoes airborne transmission between hamster pairs, albeit inefficiently as compared with SARS-CoV-2 Wuhan.
Current models of coronavirus emergence acknowledge the existence of zoonotic strains that encode spike glycoproteins with broad orthologue receptor specificities (Supplementary Fig. 3)19,20,51. For SARS-CoV and MERS-CoV, intensive farming practices involving civets and camels were critical for virus amplification and evolution, leading to repeat introductions into human populations52. As with many zoonosis, the reservoir species and index case for the SARS-CoV-2 pandemic remains elusive, although pangolins and racoon dogs remain suspects. Pangolins, which harbour distinct SARS-CoV-2-related strains2,3,4,5,6, are an endangered solitary species but are valuable as illegal trade commodities. On the basis of the broad ACE2 receptor usage of PgCoV GD, documented airborne transmission potential and efficient growth in primary nasal airway epithelial cells, we suggest that individual pangolins, or perhaps some other rare wildlife species, was productively infected and served as a nearly untraceable pass-through species that transmitted virus to humans. Such rare events might be enhanced if early human cases were immunosuppressed, potentially generating complex mutational variants during persistent infections53,54,55 (Supplementary Fig. 3).
Although speculative, a variety of data support this hypothesis. Recombinant PgCoV GD and SARS-CoV-2 share a similar host range44,56 and use pangolin, human and many other mammalian ACE2 receptors for docking and entry47,50. PgCoV, but not SARS-CoV WA1, can efficiently utilize the murine ACE2 receptor for entry and can infect mice18. However, several SARS-CoV-2 variants of concern (VOC) quickly evolved the N501Y mutation that enhanced growth in mice17,57. In contrast, the PgCoV Q498H RBD mutation probably influenced rodent susceptibility, as a similar Q498Y mutation in the SARS-CoV-2 WA1 strain conferred mACE2 receptor usage and growth in the laboratory mouse18,44. Humans have transmitted SARS-CoV-2 to other mammals, such as minks, tigers and American white-tailed deer58,59,60. As mink to human and deer to human transmission is documented, new animal reservoirs are emerging that support the evolution and transmission of new SARS-CoV-2 variants61,62,63,64. Sarbecovirus infections in pangolins and other susceptible mammals (for example, racoon dogs, civets) are understudied. One recent report noted that among the 334 pangolins confiscated in Malaysia from 2009 to 2019, none contained CoV RNA, while another study in confiscated pangolins showed 4/163 infected with virus45. Clearly, pangolins are not a heavily colonized reservoir species but rather, could serve as rare solitary hosts for zoonotic SARS-CoV-2-like viruses. The extensive overlap in potential host range suggests that the SARS-CoV-2 and PgCoV GD strains probably emerged from a common source, most probably generalists in bats. Given the close homology between PgCoV GD RBD and the different Banal strains (Supplementary Fig. 1d), our data support earlier hypotheses that Sarbecoviruses in pangolins present an emerging threat to human health (Supplementary Fig. 3).
Protease cleavage sites in the S glycoprotein play critical roles in coronavirus entry, pathogenesis and transmission. In contrast to SARS-CoV-2, PgCoV GD and related strains lack the S1/S2 RRAR polybasic furin cleavage site (682–686). Despite this deficiency, PgCoV GD infectivity and early growth phenotypes were more efficient that SARS-CoV-2 in cell lines overexpressing furin. Of interest, another furin cleavage site (K814) in SARS-CoV-2, RaTG13 and PgCoV GD regulated S2 into S2’ cleavage and infectivity27,28. Removal of the SARS-CoV-2 S1/S2 polybasic cleavage site (682–686) also enhanced virus titre stability and growth in Vero cells in vitro. Interpretation of loss-of-function studies is complicated by off-target effects on spike protein/structure–function and unknown epistatic interactions65. For example, the SARS-CoV-2 S1/S2 furin cleavage site resides within a disordered region flanked by serine residues that can be phosphorylated by proline-directed and basophilic kinases to inhibit proteolytic cleavage66. In addition, an upstream O-linked glycosylation site also regulates cleavage efficacy67,68. Loss-of-function studies deleting the RRAR furin cleavage site or upstream QTQTN residues that contain the O-linked glycosylation site resulted in virus mutants with reduced in vitro replication in airway epithelial cells and in vivo pathogenicity and/or transmission14,69. However, these loss-of-function mutations also attenuated cleavage at overlapping S1/S2 protease sites such as TMPRSS2, complicating data interpretation13. In the absence of the S1/S2 polybasic furin cleavage site, PgCoV GD replicated as efficiently as SARS-CoV-2 in primary human nasal airway epithelial and large bronchial airway epithelial cells. In contrast, PgCoV GD replication efficiency in AT2 organoids ranged from equally efficient to less efficient than SARS-CoV-2 in a donor-dependent manner. The 2003 SARS-CoV strain, which lacks an S1/S2 furin cleavage site, replicated efficiently in primary human epithelial cells, was efficiently transmitted between civets, raccoon dogs and humans, and caused ~10% mortality rates in humans70. As humans living and working near bat hibernacula experience occasional bat SARS-like CoV infections71, it appears that Sarbecoviruses use multiple protease-dependent pathways to ensure efficient transmission.
Antigenic variation in spike (~3%) has reduced vaccine and therapeutic antibody performance against SARS-CoV-272,73. Although the PgCoV GD spike is ~10% different from SARS-CoV-2 spike, the virus encode many conserved neutralizing epitopes in the RBD but not in the NTD domain. Importantly, COVID-19 convalescent patient sera and many commercial therapeutic antibodies effectively neutralized PgCoV in vitro. Unlike another PgCoV strains44, PgCoV GD caused lethal outcomes in K18-hACE2 mice, but not in wild-type mice. Using these models, the pan-Sarbecovirus mAB ADG-2 (ref. 37), hmAB DH1047 and SARS-CoV-2 S2P recombinant protein vaccines protected animals from PgCoV GD replication and disease. Thus, the current COVID-19 vaccines and immunotherapeutics should protect humans from PgCoV GD infection. Although the mechanism is probably complex, competition studies in primary human cells demonstrated reduced PgCoV GD fitness as compared with SARS-CoV-2.
Virus discovery is critical for One Health, policy, global health preparedness and the development of broadly effective countermeasures. Field collections, coupled with virus recovery in cultured cells and animals oftentimes under Biosafety Level 2 (BSL2) conditions, are historically used to identify viruses in the natural world74,75. Today, high-throughput sequencing technologies and the use of molecular approaches and synthetic genomics to recover recombinant viruses from wild-type genome length sequences have altered this historic paradigm76. The PgCoV GD strain represents a replica of a wild-type PgCoV strain circulating in nature and was isolated by researchers wearing protective clothing in a highly controlled BSL3 environment. No genetic approaches were implemented to alter the intrinsic biology and pathogenesis of the virus. Lacking ORF7, PgCoV GD derivatives encoding reporter genes such as GFP and nLUC were generated to provide high-throughput diagnostic reagents to understand virus tropism, neutralization and countermeasure performance. The paradigm allows for the recovery of a single, targeted virus from a virologically complex field sample in a safe, controlled BSL3 facility, thereby reducing the risk of the inadvertent recovery of some unidentified contaminant41,57,77. Although the PgCoV GD strain provides an optimal model to evaluate the impact of a furin cleavage site at the S1/S2 border on replication, pathogenesis and transmission, we urge constraint as such studies should require transparent, independent and rigorous review57.
This research complies with all relevant ethical regulations, and study protocols were approved by institutional boards and committees. Synthetic reconstruction of the authentic, wild-type Pangolin Coronavirus GD virus strain was based on its published sequence and undertaken with the approval of the Institutional Biosafety Committee of the University of North Carolina at Chapel Hill (Application no. 80022). All constructs were handled under proper biosafety level conditions as specified by the institutional approval, with reconstructed virus handled exclusively at biosafety level 3, with personnel wearing full-body personal protective equipment and HEPA-filtered respiratory protection. Current US Gain-of-Function (GOF) regulations note that wild-type pathogens that are circulating in or have been recovered from nature are not enhanced potential pandemic pathogens (PPPs), regardless of their pandemic potential.
Simian kidney cell line Vero-E6 (ATCC, CRL1586) was maintained in Eagle’s minimum essential medium (Gibco) supplemented with 10% fetal calf serum (FBS, Hyclone). A clonal furin-overexpressing Vero cell line was generated by using the Sleeping Beauty Transposon System as previously described by our group29. Briefly, the plasmid encoding the furin gene was transfected into the Vero cells, and the cells were selected by appropriate antibiotics. A549-ACE2 cells were maintained in Dulbecco’s modified Eagle medium (DMEM, Gibco) with 10% FBS. Using the Sleeping Beauty Transposon System34, we expressed a panel of ACE2 orthologues from the bat Rhinolophus affinis (GenBank, MT394215.1), the bat R. alcyone (GenBank, ALJ94035), camel (GenBank, XM_010993415.2), cat (GenBank, AB211997), cow (GenBank, XM_024987850.1), ferret (GenBank, NM_001310190), horse (GenBank, XM_001490191.5), rabbit (XM_002719845.3), hamster (GenBank, XM_005074209.3), dog (GenBank, NM_001165260.1), pangolin (GenBank, XM_017650263.2), mouse (GenBank, NP_001123985.1), mink (GenBank, XP_044091952) and civet (GenBank, AY881174) in non-permissive DBT cells. All ACE2 constructs included a FLAG tag at the C terminus, except the mouse and R. alcyone receptors. After selection with puromycin or gentamycin (mink, R. alcyone, mouse), ACE2 expression was confirmed by western blot analysis using a polyclonal ACE2 antiserum.
The primary human pulmonary cell cultures were purchased from the UNC Marsico Lung Institute cell bank. The generation of primary human pulmonary cell cultures was described previously22. Primary HNE cells, human bronchial epithelial LAE cells and bronchiolar (small airway epithelial (SAE)) cells were isolated from freshly excised normal human tissues obtained from transplant donors under UNC Institutional Review Board (IRB)-approved protocol (no. 03-1396) and cultured in air–liquid interface (ALI) media, as previously described22,78. The age and gender of the donors included males and females 21–70 years of age.
The molecular clone-derived SARS-CoV-2 WT and nLuc viruses were generated as previously described18,22. Briefly, genomic cDNA sequences of these viruses were ligated in vitro and subjected to T7 transcription. The transcribed RNA was electroporated into Vero cells. SARS-CoV-2-∆PRRA was generated by deleting the four residues ‘PRAA’ coding sequence from the SARS-CoV-2 WT infectious cDNA clone. All work was conducted in a high-containment BSL3 laboratory and personnel wore powered air purifying respirator (PAPR), Tyvek suits, aprons and booties, and were double gloved.
The antiviral activity of RDV, EIDD and nirmatrelvir against SARS-CoV-2-related pangolin-CoV GD (MP789) was measured in a A549-hACE2 cell line-based assay. The human cell line, A549-hACE2, was maintained in DMEM (Gibco), 20% fetal bovine serum (Hyclone) and 1X antibiotic-antimycotic (A.A., Gibco). At 24 h after plating 2 × 104 cells per well, fresh medium was added. In triplicate, cells were exposed to serial dilutions of compound and control compounds (RDV, EIDD, nirmatrelvir) in ‘infection medium’ (that is, modified growth medium similar to that above but with 5% FBS) and immediately infected for 1 h with SARS-CoV-2-related pangolin-CoV nLUC added at 800 p.f.u. per well. Virus was then removed, and cells were rinsed once and infection medium containing dilutions of drug or vehicle was added. At 24 h post infection, virus replication was measured by nLUC assay (Promega) and then read on a Promega Glomax plate reader (Promega). The IC50 value was defined in Graphpad Prism 9.0 (Graphpad) as the concentration at which there was a 50% decrease in viral replication, using non-infected wells (100% inhibition) and infected and vehicle-treated wells (0% inhibition) as controls. The experiment was repeated once.
Seven cDNA fragments covering the entire genome of the PgCoV GD strain (GISAID accession no. EPI_ISL_410721) were chemically synthesized by Bio Basic. Fragments were cloned into vector plasmid pUC57, and junctions were divided by non-palindromic sites BsaI (GGTCTCN^NNNN), BsmBI (CGTCTCN^NNNN) or SapI (GCTCTTCN^NNN) with unique 3- or 4-nucleotide cohesive ends. The cohesive ends in each fragment are indicated in Fig. 1a. To assist the replication of the full-length viral genome, we introduced a T7-promoter sequence into the 5’-end of fragment A and a 25 nt poly-A tail into the 3’-end of the fragment G. Each fragment was verified by Sanger sequencing. To enhance the efficiency of recovering the PgCoV in the cell culture, an sgRNA-N construct encoding a 5’ leader sequence, an N gene, a 3’ UTR and a 25 nt poly-A tail was assembled downstream of a T7 promoter. Two reporter viruses containing a GFP or an nLuc gene were generated by replacing part of the ORF7a gene with the reporter genes.
Seven genomic cDNA fragments were digested with appropriate restriction enzymes, resolved in 1% agarose gels, excised and purified using a QIAquick Gel Extraction kit (Qiagen). A full-length genomic cDNA was obtained by ligating seven fragments in an equal molar ratio using T4 DNA ligase (NEB). The full-length PgCoV genomic RNA or sgRNA-N was synthesized using the T7 mMESSAGE T7 transcription kit (Thermo Fisher) at 32 °C for 5 h. The RNA transcripts were then mixed and electroporated into Vero-E6 cells at 450 V and 50 µF using a Gene Pulser II electroporator (Bio-Rad). The cells were cultured as usual in the medium for 3–4 d.
Monolayers of Vero-E6 and Vero-furin cells were cultured in black-walled 96-well plates (Corning 3904) overnight. The cells were infected with SARS-CoV-2-nLuc or PgCoV-nLuc viruses at an MOI of 1. After incubation for 1 h, the inoculum was removed, the cells were washed two times with PBS and maintained in DMEM containing 5% FBS and the mixture of SARS-CoV-2 nAbs REGN10933 and REGN10987 at a concentration of 1,000 times the IC50 against SARS-CoV-2. After incubation at 37 °C for 8 h, viral entry was quantified by measuring nLuc activity using a Nano-Glo luciferase assay system (Promega) according to manufacturer specifications.
Monolayers of ACE2 orthologue-expressing A549 cell lines were cultured in black-walled 96-well plates (Corning 3904) overnight. The cells were infected with a PgCoV-GFP virus at an MOI of 0.5. For fluorescence imaging, the cells were seeded in clear-walled 96-well plates (Corning, 3598) and images captured in a Nikon fluorescence scope. Experiments were repeated twice.
Extracellular PgCoV and SARS-CoV-2 virions were collected from Vero or Vero-furin cells. Samples were lysed with modified RIPA buffer and inactivated at 98 °C. Protein samples were electrophoresed in 4–20% continuous SDS–PAGE gel (Bio-Rad) and transferred onto a PVDF membrane (Bio-Rad). Spike proteins of both viruses were probed using a mAb targeting a conserved region in S2 (Abcam, ab272504), and the N protein was probed using a mouse antiserum produced in our laboratory. The western blot images were captured and quantified using the Thermo Fisher iBright imaging system and software.
Vero-E6 cells were infected with SARS-CoV-2 isolate, icPgCoV-WT, icPgCoV-GFP, icPgCoV-nLuc or mock at an MOI of 1. At 24 h post infection, we extracted the total cellular RNA using TRIzol reagent (Thermo Fisher). Poly-A-containing messenger RNA was isolated from the total RNA using an Oligotex mRNA mini kit (Qiagen). Messenger RNA (0.6–0.7 μg) was separated on an agarose gel and transferred to a BrightStar-Plus membrane using a NorthernMax-Gly kit (Invitrogen). Blots were hybridized with a biotin-labelled oligomer (5’-BiodT/GGCTCTGTTGGGAATGTTTTGTATGCG/BiodT-3’), then detected with a chemiluminescent nucleic acid detection module (Thermo Fisher) using the iBright western blot imaging system (Thermo Fisher).
Vero-E6 cells were plated at 20,000 cells per well in black-walled 96-well plates (Corning 3904). Human serum samples were tested at a starting dilution of 1:40 and mAb samples were tested at a starting concentration of 30 to 0.1 μg ml−1 and were serially diluted 3-fold for up to eight dilution spots. Diluted antibodies and sera were then mixed with 200 p.f.u. per well of PgCoV-nLuc or SARS-nLuc virus and the mixtures were incubated at 37 °C with 5% CO2 for 1 h. Following incubation, growth media were removed and virus–antibody mixtures were added to the cells in duplicate. Virus-only controls were included in each plate and all samples were run in duplicate. Following infection, plates were incubated at 37 °C with 5% CO2 for 48 h. After the 48 h incubation, cells were lysed and luciferase activity was measured via a Nano-Glo luciferase assay system (Promega) according to manufacturer specifications. Neutralization titres were defined as the sample dilution at which a 50% reduction in the relative light units (RLU) was observed relative to the average of the virus control wells.
Hamster studies were performed in accordance with University of Wisconsin-Madison Institutional Animal Care and Use Committee (IACUC) protocol no. V006426. Syrian hamsters (females, 4–6 weeks old) were purchased from Envigo and allowed to acclimate for a minimum of 3 d at BSL3 agriculture containment at the Influenza Research Institute (University of Wisconsin). Hamsters were infected with 103 p.f.u. of PgCoV GD or a D-form SARS-CoV-2 strain NCGM02 intranasally under isoflurane anaesthesia. At the indicated timepoints, a subset of hamsters was euthanized by deep anaesthesia using isoflurane inhalation and cervical dislocation, and tissue samples were collected for virus titre. No animals were excluded from the study and animals were randomly assigned to groups. Sections were stained and scored by pathologists blinded to the experimental groups.
To evaluate indirect virus transmission between hamsters, groups of hamsters (n = 6 per group) were infected with 103 p.f.u. of PgCoV or SARS-CoV-2 viruses intranasally under isoflurane anaesthesia. Infected animals were placed in specially designed cages inside an isolator unit34,35. After 24 h, naïve hamsters were placed on the other side of the cage, with 5-cm separation by a double-layered divider to allow free airflow. The isolator unit provided one-directional airflow; therefore, the infected hamsters were placed in the front of the isolator unit. Metal shrouds were placed over the cages so that only the front and back of the cage was open. Nasal washes were collected at 3-d intervals for the infected hamsters and 2-d intervals for the exposed animals starting on day 2 after infection or exposure.
Mouse studies complied with all relevant ethical regulations and were performed after UNC IACUC protocol no. 20-074 review and approval. For the mAB ADG-2 prophylactic study, 10-week-old female BALB/c mice (Envigo, 047) were inoculated with 200 µg of indicated antibodies intraperitoneally 12 h before infection. For the SARS-2 S2P vaccine study, 10-week-old female BALB/c mice (Envigo, 047) were primed with 500 μg of S2P protein in 10 µg aluminum hydroxide suspension intramuscularly and boosted 4 weeks later. ADG-2-injected or S2P-vaccinated mice were anaesthetized with a mixture of ketamine/xylazine and infected with 104 p.f.u. of either PgCoV or mouse-adapted SARS-CoV-2 (MA10). Alternatively, hACE2-K18 mice were prophylactically treated with mAB DH1047 at 200 μg intraperitoneally or with control isotype and infected with PgCoV GD. Mice were monitored daily for clinical signs of disease (weight loss) and mortality. At indicated timepoints, mice were euthanized via isoflurane overdose and samples for titre and histopathological analyses were collected. Titre samples were collected in PBS with glass beads and stored at −80 °C until further processing; viral lung, nasal or brain titres were determined via plaque assay as described above. Histological samples were stored in 10% formalin at 4 °C for 7 d and formalin was replaced with fresh 10% formalin before removal out of the BSL3. No animals were excluded from the study and animals were randomly assigned to groups.
Primary airway cell cultures and AT2 organoids were fixed twice for 30 min in 4% formaldehyde in PBS, and washed and stored in PBS. The viral signals were stained by a mouse antiserum against SARS-CoV-2 N protein (1:500 dilution) using species-specific secondary antibodies as previously described22. The HNE and LAE cultures were also imaged for α-tubulin (Millipore MAB1864; 3 μg ml−1), MUC5AC (Thermo Scientific 45M1; 4 μg ml−1), MUC5B (polyclonal rabbit against a MUC5B peptide (MAN5BII), 1:1,000) and CCSP (Sigma 07-623; 1:2,000) as indicated. Filamentous actin was localized with phalloidin (Invitrogen, A22287) and DNA with Hoechst 33342 (Invitrogen). An Olympus FV3000RS confocal microscope in Galvo scan mode was used to acquire 5-channel Z stacks by the 2-phase sequential scan. Representative stacks were acquired with a ×60 oil objective (XYZ = 212 μm × 212 μm × ~25 μm) and shown as Z-projections or single-slice, XZ cross-sections to distinguish individual cell features to characterize the infected cell types. A ×20 objective was used to acquire 2D, single-channel, apical snapshots of nine fields (636 μm × 636 μm; combined area of 3.64 mm2), selected in evenly spaced grids across each sham-infected donor culture, and ImageJ was used to measure the relative apical culture surface covered by multiciliated cells.
No statistical methods were used to pre-determine sample sizes, but our sample sizes are similar to those reported in previous publications from our group (18, 19, 20, 22). Technical and biological replicates are described in the figure legends. Data distribution was assumed to be normal, but this was not formally tested. For analysis of PgCoV and SARS-CoV-2 nLuc virus entry (Fig. 1e), we applied an unpaired t-test. For COVID-19 patient sera cross-neutralization against PgCoV and SARS-CoV-2 (Fig. 5e), we applied paired t-test. For analysis of mouse studies, the comparison of weight-change curves was performed using a repeated measurements two-way ANOVA with Tukey’s post hoc test using Prism v.9.0 (GraphPad). Results of transmission studies were compared using Fischer exact test, while paired t-test and unpaired t-test were used in Figs. 5e and 1e, respectively.
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
The experiment data that support the findings of this study are available from the corresponding author upon request. In addition, we will provide Biodefense and Emerging Infections Research Resources Repository (BEI) with the recombinant wild-type PgCoV virus, including indicator viruses encoding nLUC and GFP. Users would be requested to demonstrate approval for receiving the virus from their institutional biosafety group and work with the virus under Biosafety Level 3 (BSL3) conditions. Requests for the PgCoV molecular clone should be directed to the corresponding author and will require a description of the planned experimental studies listed under a Material Transfer Agreement. As experiments of concern could be performed with this construct, we will review the request for the molecular clone with NIH for feedback and approval before providing the reagent. The materials transfer timeline shall be determined by the appropriate time for reviewing the documentations by the University of North Carolina and NIH. Source data are provided with this paper.
Worobey, M. et al. The Huanan Seafood Wholesale Market in Wuhan was the early epicenter of the COVID-19 pandemic. Science 377, 951–959 (2022).
Article CAS PubMed Google Scholar
Holmes, E. C. et al. The origins of SARS-CoV-2: a critical review. Cell 184, 4848–4856 (2021).
Article CAS PubMed PubMed Central Google Scholar
Xiao, K. et al. Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins. Nature 583, 286–289 (2020).
Article CAS PubMed Google Scholar
Lam, T. T.-Y. et al. Identifying SARS-CoV-2-related coronaviruses in Malayan pangolins. Nature 583, 282–285 (2020).
Article CAS PubMed Google Scholar
Zhang, T., Wu, Q. & Zhang, Z. Probable pangolin origin of SARS-CoV-2 associated with the COVID-19 outbreak. Curr. Biol. 30, 1346–1351.e2 (2020).
Article CAS PubMed PubMed Central Google Scholar
Liu, P. et al. Are pangolins the intermediate host of the 2019 novel coronavirus (SARS-CoV-2)? PLoS Pathog. 16, e1008421 (2020).
Article CAS PubMed PubMed Central Google Scholar
Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273 (2020).
Article CAS PubMed PubMed Central Google Scholar
Hu, D. et al. Genomic characterization and infectivity of a novel SARS-like coronavirus in Chinese bats. Emerg. Microbes Infect. 7, 154 (2018).
Article PubMed PubMed Central Google Scholar
Zhou, H. et al. A novel bat coronavirus closely related to SARS-CoV-2 contains natural insertions at the S1/S2 cleavage site of the spike protein. Curr. Biol. 30, 2196–2203.e3 (2020).
Article CAS PubMed PubMed Central Google Scholar
Delaune, D. et al. A novel SARS-CoV-2 related coronavirus in bats from Cambodia. Nat. Commun. 12, 6563 (2021).
Article CAS PubMed PubMed Central Google Scholar
Temmam, S. et al. Bat coronaviruses related to SARS-CoV-2 and infectious for human cells. Nature 604, 330–336 (2022).
Article CAS PubMed Google Scholar
Wrobel, A. G. et al. SARS-CoV-2 and bat RaTG13 spike glycoprotein structures inform on virus evolution and furin-cleavage effects. Nat. Struct. Mol. Biol. 27, 763–767 (2020).
Article CAS PubMed PubMed Central Google Scholar
Johnson, B. A. et al. Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis. Nature 591, 293–299 (2021).
Article CAS PubMed PubMed Central Google Scholar
Zhu, Y. et al. A genome-wide CRISPR screen identifies host factors that regulate SARS-CoV-2 entry. Nat. Commun. 12, 961 (2021).
Article CAS PubMed PubMed Central Google Scholar
Andersen, K. G., Rambaut, A., Lipkin, W. I., Holmes, E. C. & Garry, R. F. The proximal origin of SARS-CoV-2. Nat. Med. 26, 450–452 (2020).
Article CAS PubMed PubMed Central Google Scholar
Liu, K. et al. Binding and molecular basis of the bat coronavirus RaTG13 virus to ACE2 in humans and other species. Cell 184, 3438–3451.e10 (2021).
Article CAS PubMed PubMed Central Google Scholar
Gu, H. et al. Adaptation of SARS-CoV-2 in BALB/c mice for testing vaccine efficacy. Science 369, 1603–1607 (2020).
Article CAS PubMed PubMed Central Google Scholar
Dinnon, K. H. et al. A mouse-adapted model of SARS-CoV-2 to test COVID-19 countermeasures. Nature 586, 560–566 (2020).
Article PubMed PubMed Central Google Scholar
Menachery, V. D. et al. A SARS-like cluster of circulating bat coronaviruses shows potential for human emergence. Nat. Med. 21, 1508–1513 (2015).
Article CAS PubMed PubMed Central Google Scholar
Menachery, V. D. et al. SARS-like WIV1-CoV poised for human emergence. Proc. Natl Acad. Sci. USA 113, 3048–3053 (2016).
Article CAS PubMed PubMed Central Google Scholar
Chan, Y. A. & Zhan, S. H. Single source of pangolin CoVs with a near identical Spike RBD to SARS-CoV-2. Preprint at bioRxiv https://doi.org/10.1101/2020.07.07.184374 (2020).
Hou, Y. J. et al. SARS-CoV-2 reverse genetics reveals a variable infection gradient in the respiratory tract. Cell 182, 429–446.e14 (2020).
Article CAS PubMed PubMed Central Google Scholar
Yount, B. et al. Reverse genetics with a full-length infectious cDNA of severe acute respiratory syndrome coronavirus. Proc. Natl Acad. Sci. USA 100, 12995–13000 (2003).
Article CAS PubMed PubMed Central Google Scholar
Yount, B., Denison, M. R., Weiss, S. R. & Baric, R. Systematic assembly of a full-length infectious cDNA of mouse hepatitis virus strain A59. J. Virol. 76, 11065–11078 (2002).
Article CAS PubMed PubMed Central Google Scholar
Scobey, T. et al. Reverse genetics with a full-length infectious cDNA of the Middle East respiratory syndrome coronavirus. Proc. Natl Acad. Sci. USA 110, 16157–16162 (2013).
Article CAS PubMed PubMed Central Google Scholar
Peacock, T. P. et al. The furin cleavage site in the SARS-CoV-2 spike protein is required for transmission in ferrets. Nat. Microbiol. 6, 899–909 (2021).
Article CAS PubMed Google Scholar
Zhang, Y. et al. A second functional furin site in the SARS-CoV-2 spike protein. Emerg. Microbes Infect. 11, 182–194 (2022).
Article CAS PubMed PubMed Central Google Scholar
Coutard, B. et al. The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade. Antivir. Res. 176, 104742 (2020).
Article CAS PubMed Google Scholar
Tse, L. V. et al. Generation of mature DENVs via genetic modification and directed evolution. mBio 13, e00386-22 (2022).
Article PubMed PubMed Central Google Scholar
Sheahan, T. P. et al. An orally bioavailable broad-spectrum antiviral inhibits SARS-CoV-2 in human airway epithelial cell cultures and multiple coronaviruses in mice. Sci. Transl. Med. 12, eabb5883 (2020).
Article CAS PubMed Google Scholar
Sheahan, T. P. et al. Broad-spectrum antiviral GS-5734 inhibits both epidemic and zoonotic coronaviruses. Sci. Transl. Med. 9, eaal3653 (2017).
Article PubMed PubMed Central Google Scholar
Pruijssers, A. J. et al. Remdesivir inhibits SARS-CoV-2 in human lung cells and chimeric SARS-CoV expressing the SARS-CoV-2 RNA polymerase in mice. Cell Rep. 32, 107940 (2020).
Article CAS PubMed PubMed Central Google Scholar
Katsura, H. et al. Human lung stem cell-based alveolospheres provide insights into SARS-CoV-2-mediated interferon responses and pneumocyte dysfunction. Cell Stem Cell 27, 890–904.e8 (2020).
Article CAS PubMed PubMed Central Google Scholar
Hou, Y. J. et al. SARS-CoV-2 D614G variant exhibits efficient replication ex vivo and transmission in vivo. Science 370, 1464–1468 (2020).
Article CAS PubMed PubMed Central Google Scholar
Imai, M. et al. Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature 486, 420–428 (2012).
Article CAS PubMed PubMed Central Google Scholar
Imai, M. et al. Syrian hamsters as a small animal model for SARS-CoV-2 infection and countermeasure development. Proc. Natl Acad. Sci. USA 117, 16587–16595 (2020).
Article CAS PubMed PubMed Central Google Scholar
Rappazzo, C. G. et al. Broad and potent activity against SARS-like viruses by an engineered human monoclonal antibody. Science 371, 823–829 (2021).
Article CAS PubMed PubMed Central Google Scholar
DiPiazza, A. T. et al. COVID-19 vaccine mRNA-1273 elicits a protective immune profile in mice that is not associated with vaccine-enhanced disease upon SARS-CoV-2 challenge. Immunity 54, 1869–1882.e6 (2021).
Article CAS PubMed PubMed Central Google Scholar
Corbett, K. S. et al. SARS-CoV-2 mRNA vaccine design enabled by prototype pathogen preparedness. Nature 586, 567–571 (2020).
Article CAS PubMed PubMed Central Google Scholar
Zheng, J. et al. COVID-19 treatments and pathogenesis including anosmia in K18-hACE2 mice. Nature 589, 603–607 (2021).
Article CAS PubMed Google Scholar
Martinez, D. R. et al. A broadly cross-reactive antibody neutralizes and protects against sarbecovirus challenge in mice. Sci. Transl. Med. 14, eabj7125 (2022).
Article CAS PubMed PubMed Central Google Scholar
Wacharapluesadee, S. et al. Evidence for SARS-CoV-2 related coronaviruses circulating in bats and pangolins in Southeast Asia. Nat. Commun. 12, 972 (2021).
Article CAS PubMed PubMed Central Google Scholar
Liu, P., Chen, W. & Chen, J.-P. Viral metagenomics revealed sendai virus and coronavirus infection of Malayan pangolins (Manis javanica). Viruses 11, 979 (2019).
Article CAS PubMed PubMed Central Google Scholar
Liu, M.-Q. et al. A SARS-CoV-2-related virus from Malayan pangolin causes lung infection without severe disease in human ACE2-transgenic mice. J. Virol. 97, e01719–e01722 (2023).
Article PubMed PubMed Central Google Scholar
Peng, M.-S. et al. The high diversity of SARS-CoV-2-related coronaviruses in pangolins alerts potential ecological risks. Zool. Res. 42, 834–844 (2021).
Article PubMed PubMed Central Google Scholar
Li, X. et al. Evolutionary history, potential intermediate animal host, and cross-species analyses of SARS-CoV-2. J. Med. Virol. 92, 602–611 (2020).
Article CAS PubMed PubMed Central Google Scholar
Niu, S. et al. Molecular basis of cross‐species ACE2 interactions with SARS‐CoV‐2‐like viruses of pangolin origin. EMBO J. 40, e107786 (2021).
Article CAS PubMed PubMed Central Google Scholar
Zhao, X. et al. Broad and differential animal angiotensin-converting enzyme 2 receptor usage by SARS-CoV-2. J. Virol. 94, e00940–20 (2020).
Article PubMed PubMed Central Google Scholar
Li, Y. et al. SARS-CoV-2 and three related coronaviruses utilize multiple ACE2 orthologs and are potently blocked by an improved ACE2-Ig. J. Virol. 94, e01283-20 (2020).
Article CAS PubMed PubMed Central Google Scholar
Wrobel, A. G. et al. Structure and binding properties of pangolin-CoV spike glycoprotein inform the evolution of SARS-CoV-2. Nat. Commun. 12, 837 (2021).
Article CAS PubMed PubMed Central Google Scholar
Ge, X.-Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535–538 (2013).
Article CAS PubMed PubMed Central Google Scholar
Cui, J., Li, F. & Shi, Z.-L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Microbiol. 17, 181–192 (2019).
Article CAS PubMed Google Scholar
Caccuri, F. et al. Competition for dominance within replicating quasispecies during prolonged SARS-CoV-2 infection in an immunocompromised host. Virus Evol. 8, veac042 (2022).
Article PubMed PubMed Central Google Scholar
Cele, S. et al. SARS-CoV-2 prolonged infection during advanced HIV disease evolves extensive immune escape. Cell Host Microbe 30, 154–162.e5 (2022).
Article CAS PubMed PubMed Central Google Scholar
Xia, S., Wang, L., Zhu, Y., Lu, L. & Jiang, S. Origin, virological features, immune evasion and intervention of SARS-CoV-2 Omicron sublineages. Sig. Transduct. Target. Ther. 7, 241 (2022).
Article CAS Google Scholar
Letko, M., Marzi, A. & Munster, V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Nat. Microbiol. 5, 562–569 (2020).
Article CAS PubMed PubMed Central Google Scholar
Martinez, D. R. et al. Chimeric spike mRNA vaccines protect against Sarbecovirus challenge in mice. Science 373, 991–998 (2021).
Article CAS PubMed PubMed Central Google Scholar
Hammer, A. S. et al. SARS-CoV-2 Transmission between mink (Neovison vison) and humans, Denmark. Emerg. Infect. Dis. https://doi.org/10.3201/eid2702.203794 (2021).
McAloose, D. et al. From people to Panthera: natural SARS-CoV-2 infection in tigers and lions at the Bronx Zoo. mBio 11, e02220 (2020).
Article CAS PubMed PubMed Central Google Scholar
Chandler, J. C. et al. SARS-CoV-2 exposure in wild white-tailed deer (Odocoileus virginianus). Proc. Natl Acad. Sci. USA 118, e2114828118 (2021).
Article CAS PubMed PubMed Central Google Scholar
Oude Munnink, B. B. et al. Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans. Science 371, 172–177 (2021).
Article CAS PubMed Google Scholar
Pickering, B. et al. Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission. Nat. Microbiol. 7, 2011–2024 (2022).
Article CAS PubMed PubMed Central Google Scholar
Hale, V. L. et al. SARS-CoV-2 infection in free-ranging white-tailed deer. Nature 602, 481–486 (2022).
Article CAS PubMed Google Scholar
Willgert, K. et al. Transmission history of SARS-CoV-2 in humans and white-tailed deer. Sci. Rep. 12, 12094 (2022).
Article CAS PubMed PubMed Central Google Scholar
Housden, B. E. et al. Loss-of-function genetic tools for animal models: cross-species and cross-platform differences. Nat. Rev. Genet. 18, 24–40 (2017).
Article CAS PubMed Google Scholar
Örd, M., Faustova, I. & Loog, M. The sequence at Spike S1/S2 site enables cleavage by furin and phospho-regulation in SARS-CoV-2 but not in SARS-CoV1 or MERS-CoV. Sci. Rep. 10, 16944 (2020).
Article PubMed PubMed Central Google Scholar
Sanda, M., Morrison, L. & Goldman, R. N- and O-glycosylation of the SARS-CoV-2 spike protein. Anal. Chem. 93, 2003–2009 (2021).
Article CAS PubMed PubMed Central Google Scholar
Vu, M. N. et al. QTQTN motif upstream of the furin-cleavage site plays key role in SARS-CoV-2 infection and pathogenesis. Proc. Natl Acad. Sci. USA 119, e2205690119 (2022).
Article CAS PubMed PubMed Central Google Scholar
Hoffmann, M., Kleine-Weber, H. & Pöhlmann, S. A multibasic cleavage site in the spike protein of SARS-CoV-2 is essential for infection of human lung cells. Mol. Cell 78, 779–784.e5 (2020).
Article CAS PubMed PubMed Central Google Scholar
Wang, M. et al. SARS-CoV infection in a restaurant from palm civet. Emerg. Infect. Dis. https://doi.org/10.3201/eid1112.041293 (2005).
Li, H. et al. Human–animal interactions and bat coronavirus spillover potential among rural residents in Southern China. Biosaf. Health 1, 84–90 (2019).
Article PubMed PubMed Central Google Scholar
Pajon, R. et al. SARS-CoV-2 Omicron variant neutralization after mRNA-1273 booster vaccination. N. Engl. J. Med. 386, 1088–1091 (2022).
Article PubMed Google Scholar
Collie, S., Champion, J., Moultrie, H., Bekker, L.-G. & Gray, G. Effectiveness of BNT162b2 vaccine against Omicron variant in South Africa. N. Engl. J. Med. 386, 494–496 (2022).
Article PubMed Google Scholar
Hsiung, G. D. The impact of cell culture sensitivity on rapid viral diagnosis: a historical perspective. Yale J. Biol. Med. 62, 77–88 (1989).
Google Scholar
Murphy, F. A. in Advances in Virus Research Vol. 95 (eds Kielian, M.,Maramorosch, K. & Mettenleiter, T. C.) 197–220 (Elsevier, 2016).
Mokili, J. L., Rohwer, F. & Dutilh, B. E. Metagenomics and future perspectives in virus discovery. Curr. Opin. Virol. 2, 63–77 (2012).
Article CAS PubMed PubMed Central Google Scholar
Walls, A. C. et al. Elicitation of broadly protective sarbecovirus immunity by receptor-binding domain nanoparticle vaccines. Cell 184, 5432–5447.e16 (2021).
Article CAS PubMed PubMed Central Google Scholar
Fulcher, M. L. & Randell, S. H. in Epithelial Cell Culture Protocols 2nd edn (eds Randell, S. H. & Fulcher, M. L.) 109–121 (Humana Press, 2013).
Download references
This work was supported by grants from the US National Institutes of Health, including CA260543 (R.S.B.), AI151797 (R.S.B.), AI116484 (R.S.B.), AI158571 (R.S.B.), AI165977 (Y.K.), R01HL146557 (P.R.T.) R01 HL153375 (P.R.T.), AI171292 (R.S.B.) and R01 AI157155 (R.S.B. and J.E.C.), and the Center for Research on Influenza Pathogenesis and Transmission (CRIPT) 75N93021C00014 (Y.K.). D.R.M. was funded by a Burroughs Wellcome Fund Postdoctoral Enrichment Program Award and a Hanna H. Gray Fellowship from the Howard Hughes Medical Institute. This work was supported in part by the Japan Program for Infectious Diseases Research and Infrastructure from the Japan Agency for Medical Research and Development (JP21wm0125002). B. Haynes (Duke University) kindly provided DH1047 for use in this study.
Yixuan J. Hou
Present address: Moderna Inc., Cambridge, MA, USA
Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Yixuan J. Hou, Sarah R. Leist, Rita M. Meganck, David R. Martinez, Alexandra Schäfer, Nicholas J. Catanzaro, Ande West, Catlin E. Edwards, Boyd Yount, John Powers, Lily Adams, Edgar F. Kong, Melissa Mattocks & Ralph S. Baric
Influenza Research Institute, Department of Pathobiological Sciences, School of Veterinary Medicine, University of Wisconsin, Madison, WI, USA
Shiho Chiba, Peter Halfmann & Yoshihiro Kawaoka
Department of Cell Biology, Regeneration Next Initiative, Duke University Medical Center, Durham, NC, USA
Vishwaraj Sontake, Aleksandra Tata & Purushothama R. Tata
Marsico Lung Institute, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Rhianna E. Lee, Samuel C. Gallant & Scott H. Randell
Vanderbilt Vaccine Center, Vanderbilt University Medical Center, Nashville, TN, USA
Seth J. Zost & James E. Crowe Jr
Division of Virology, Department of Microbiology and Immunology, Institute of Medical Science, University of Tokyo, Tokyo, Japan
Yoshihiro Kawaoka
Department of Microbiology and Immunology, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA
Ralph S. Baric
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.J.H., Y.K. and R.S.B. conceived the study; Y.J.H., S.C., B.Y. and S.R.L. performed most of the experiments; R.M.M., V.S., C.E.E., R.E.L., S.C.G., A.T., P.H., A.S., D.R.M., N.J.C., L.A., A.W., S.J.Z., E.F.K., M.M. and J.P. assisted in experiments; S.H.R., P.R.T. and J.E.C. provided essential materials and reagents; Y.K., P.H. and R.S.B. supervised the research; J.E.C. provided key antibody reagents; Y.J.H., N.J.C. and R.S.B. wrote the manuscript and generated the illustrations; Y.K. and R.S.B. provided support, reviewed and edited the manuscript. All authors read and approved the final manuscript.
Correspondence to Ralph S. Baric.
J.E.C. has served as consultant for Luna Biologics, Merck and GlaxoSmithKline, is a member of the Scientific Advisory Board of Meissa Vaccines and is Founder of IDBiologics. R.S.B. is a member of the advisory board of VaxArt and Invivyd and has collaborations with Takeda, Pfizer, Moderna, Ridgeback Biosciences, Gilead and Eli Lily. The Crowe laboratory at Vanderbilt University Medical Center has received sponsored research agreements from Takeda Pharmaceuticals, IDBiologics and AstraZeneca. Vanderbilt University has applied for patents pertaining to some of the antibodies described here. All other authors declare no competing interests.
Nature Microbiology thanks the anonymous reviewers for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Figs. 1–3, and Tables 1 and 2.
Unprocessed western blot images for Fig. 1d.
Unprocessed northern blot images for Fig. 1c.
Unprocessed western blot images for Fig. 2b.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Hou, Y.J., Chiba, S., Leist, S.R. et al. Host range, transmissibility and antigenicity of a pangolin coronavirus. Nat Microbiol 8, 1820–1833 (2023). https://doi.org/10.1038/s41564-023-01476-x
Download citation
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41564-023-01476-x
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.