COVID-19 vaccines have been authorized in multiple countries, and more are under rapid development. Careful design of a vaccine prioritization strategy across sociodemographic groups is a crucial public policy challenge given that 1) vaccine supply will be constrained for the first several months of the vaccination campaign, 2) there are stark differences in transmission and severity of impacts from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) across groups, and 3) SARS-CoV-2 differs markedly from previous pandemic viruses. We assess the optimal allocation of a limited vaccine supply in the United States across groups differentiated by age and essential worker status, which constrains opportunities for social distancing. We model transmission dynamics using a compartmental model parameterized to capture current understanding of the epidemiological characteristics of COVID-19, including key sources of group heterogeneity (susceptibility, severity, and contact rates). We investigate three alternative policy objectives (minimizing infections, years of life lost, or deaths) and model a dynamic strategy that evolves with the population epidemiological status. We find that this temporal flexibility contributes substantially to public health goals. Older essential workers are typically targeted first. However, depending on the objective, younger essential workers are prioritized to control spread or seniors to directly control mortality. When the objective is minimizing deaths, relative to an untargeted approach, prioritization averts deaths on a range between 20,000 (when nonpharmaceutical interventions are strong) and 300,000 (when these interventions are weak). We illustrate how optimal prioritization is sensitive to several factors, most notably, vaccine effectiveness and supply, rate of transmission, and the magnitude of initial infections.

Coronaviruses are adept at evading host antiviral pathways induced by viral double-stranded RNA, including interferon (IFN) signaling, oligoadenylate synthetase–ribonuclease L (OAS-RNase L), and protein kinase R (PKR). While dysregulated or inadequate IFN responses have been associated with severe coronavirus infection, the extent to which the recently emerged SARS-CoV-2 activates or antagonizes these pathways is relatively unknown. We found that SARS-CoV-2 infects patient-derived nasal epithelial cells, present at the initial site of infection; induced pluripotent stem cell-derived alveolar type 2 cells (iAT2), the major cell type infected in the lung; and cardiomyocytes (iCM), consistent with cardiovascular consequences of COVID-19 disease. Robust activation of IFN or OAS-RNase L is not observed in these cell types, whereas PKR activation is evident in iAT2 and iCM. In SARS-CoV-2–infected Calu-3 and A549^ACE2 lung-derived cell lines, IFN induction remains relatively weak; however, activation of OAS-RNase L and PKR is observed. This is in contrast to Middle East respiratory syndrome (MERS)-CoV, which effectively inhibits IFN signaling and OAS-RNase L and PKR pathways, but is similar to mutant MERS-CoV lacking innate immune antagonists. Remarkably, OAS-RNase L and PKR are activated in MAVS knockout A549^ACE2 cells, demonstrating that SARS-CoV-2 can induce these host antiviral pathways despite minimal IFN production. Moreover, increased replication and cytopathic effect in RNASEL knockout A549^ACE2 cells implicates OAS-RNase L in restricting SARS-CoV-2. Finally, while SARS-CoV-2 fails to antagonize these host defense pathways, which contrasts with other coronaviruses, the IFN signaling response is generally weak. These host–virus interactions may contribute to the unique pathogenesis of SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) research and antiviral discovery are hampered by the lack of a cell-based virus replication system that can be readily adopted without biosafety level 3 (BSL-3) restrictions. Here, the construction of a noninfectious SARS-CoV-2 reporter replicon and its application in deciphering viral replication mechanisms and evaluating SARS-CoV-2 inhibitors are presented. The replicon genome is replication competent but does not produce progeny virions. Its replication can be inhibited by RdRp mutations or by known SARS-CoV-2 antiviral compounds. Using this system, a high-throughput antiviral assay has also been developed. Significant differences in potencies of several SARS-CoV-2 inhibitors in different cell lines were observed, which highlight the challenges of discovering antivirals capable of inhibiting viral replication in vivo and the importance of testing compounds in multiple cell culture models. The generation of a SARS-CoV-2 replicon provides a powerful platform to expand the global research effort to combat COVID-19.

Increasing evidence indicates that superspreading plays a dominant role in COVID-19 transmission. Recent estimates suggest that the dispersion parameter k for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is on the order of 0.1, which corresponds to about 10% of cases being the source of 80% of infections. To investigate how overdispersion might affect the outcome of various mitigation strategies, we developed an agent-based model with a social network that allows transmission through contact in three sectors: “close” (a small, unchanging group of mutual contacts as might be found in a household), “regular” (a larger, unchanging group as might be found in a workplace or school), and “random” (drawn from the entire model population and not repeated regularly). We assigned individual infectivity from a gamma distribution with dispersion parameter k. We found that when k was low (i.e., greater heterogeneity, more superspreading events), reducing random sector contacts had a far greater impact on the epidemic trajectory than did reducing regular contacts; when k was high (i.e., less heterogeneity, no superspreading events), that difference disappeared. These results suggest that overdispersion of COVID-19 transmission gives the virus an Achilles’ heel: Reducing contacts between people who do not regularly meet would substantially reduce the pandemic, while reducing repeated contacts in defined social groups would be less effective.

RNA-dependent RNA polymerases (RdRps) of the Nidovirales (Coronaviridae, Arteriviridae, and 12 other families) are linked to an amino-terminal (N-terminal) domain, called NiRAN, in a nonstructural protein (nsp) that is released from polyprotein 1ab by the viral main protease (M^pro). Previously, self-GMPylation/UMPylation activities were reported for an arterivirus NiRAN-RdRp nsp and suggested to generate a transient state primed for transferring nucleoside monophosphate (NMP) to (currently unknown) viral and/or cellular biopolymers. Here, we show that the coronavirus (human coronavirus [HCoV]-229E and severe acute respiratory syndrome coronavirus 2) nsp12 (NiRAN-RdRp) has Mn²⁺-dependent NMPylation activity that catalyzes the transfer of a single NMP to the cognate nsp9 by forming a phosphoramidate bond with the primary amine at the nsp9 N terminus (N3825) following M^pro-mediated proteolytic release of nsp9 from N-terminally flanking nsps. Uridine triphosphate was the preferred nucleotide in this reaction, but also adenosine triphosphate, guanosine triphosphate, and cytidine triphosphate were suitable cosubstrates. Mutational studies using recombinant coronavirus nsp9 and nsp12 proteins and genetically engineered HCoV-229E mutants identified residues essential for NiRAN-mediated nsp9 NMPylation and virus replication in cell culture. The data corroborate predictions on NiRAN active-site residues and establish an essential role for the nsp9 N3826 residue in both nsp9 NMPylation in vitro and virus replication. This residue is part of a conserved N-terminal NNE tripeptide sequence and shown to be the only invariant residue in nsp9 and its homologs in viruses of the family Coronaviridae. The study provides a solid basis for functional studies of other nidovirus NMPylation activities and suggests a possible target for antiviral drug development.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a beta-CoV that recently emerged as a human pathogen and is the causative agent of the COVID-19 pandemic. A molecular framework of how the virus manipulates host cellular machinery to facilitate infection remains unclear. Here, we focus on SARS-CoV-2 NSP1, which is proposed to be a virulence factor that inhibits protein synthesis by directly binding the human ribosome. We demonstrate biochemically that NSP1 inhibits translation of model human and SARS-CoV-2 messenger RNAs (mRNAs). NSP1 specifically binds to the small (40S) ribosomal subunit, which is required for translation inhibition. Using single-molecule fluorescence assays to monitor NSP1–40S subunit binding in real time, we determine that eukaryotic translation initiation factors (eIFs) allosterically modulate the interaction of NSP1 with ribosomal preinitiation complexes in the absence of mRNA. We further elucidate that NSP1 competes with RNA segments downstream of the start codon to bind the 40S subunit and that the protein is unable to associate rapidly with 80S ribosomes assembled on an mRNA. Collectively, our findings support a model where NSP1 proteins from viruses in at least two subgenera of beta-CoVs associate with the open head conformation of the 40S subunit to inhibit an early step of translation, by preventing accommodation of mRNA within the entry channel.

It was recently shown that the major genetic risk factor associated with becoming severely ill with COVID-19 when infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is inherited from Neandertals. New, larger genetic association studies now allow additional genetic risk factors to be discovered. Using data from the Genetics of Mortality in Critical Care (GenOMICC) consortium, we show that a haplotype at a region on chromosome 12 associated with requiring intensive care when infected with the virus is inherited from Neandertals. This region encodes proteins that activate enzymes that are important during infections with RNA viruses. In contrast to the previously described Neandertal haplotype that increases the risk for severe COVID-19, this Neandertal haplotype is protective against severe disease. It also differs from the risk haplotype in that it has a more moderate effect and occurs at substantial frequencies in all regions of the world outside Africa. Among ancient human genomes in western Eurasia, the frequency of the protective Neandertal haplotype may have increased between 20,000 and 10,000 y ago and again during the past 1,000 y.

The contributions of asymptomatic infections to herd immunity and community transmission are key to the resurgence and control of COVID-19, but are difficult to estimate using current models that ignore changes in testing capacity. Using a model that incorporates daily testing information fit to the case and serology data from New York City, we show that the proportion of symptomatic cases is low, ranging from 13 to 18%, and that the reproductive number may be larger than often assumed. Asymptomatic infections contribute substantially to herd immunity, and to community transmission together with presymptomatic ones. If asymptomatic infections transmit at similar rates as symptomatic ones, the overall reproductive number across all classes is larger than often assumed, with estimates ranging from 3.2 to 4.4. If they transmit poorly, then symptomatic cases have a larger reproductive number ranging from 3.9 to 8.1. Even in this regime, presymptomatic and asymptomatic cases together comprise at least 50% of the force of infection at the outbreak peak. We find no regimes in which all infection subpopulations have reproductive numbers lower than three. These findings elucidate the uncertainty that current case and serology data cannot resolve, despite consideration of different model structures. They also emphasize how temporal data on testing can reduce and better define this uncertainty, as we move forward through longer surveillance and second epidemic waves. Complementary information is required to determine the transmissibility of asymptomatic cases, which we discuss. Regardless, current assumptions about the basic reproductive number of severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) should be reconsidered.

The majority of currently circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viruses have mutant spike glycoproteins that contain the D614G substitution. Several studies have suggested that spikes with this substitution are associated with higher virus infectivity. We use cryo-electron microscopy to compare G614 and D614 spikes and show that the G614 mutant spike adopts a range of more open conformations that may facilitate binding to the SARS-CoV-2 receptor, ACE2, and the subsequent structural rearrangements required for viral membrane fusion.

The characterization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral kinetics in hospitalized patients and its association with mortality is unknown. We analyzed death and nasopharyngeal viral kinetics in 655 hospitalized patients from the prospective French COVID cohort. The model predicted a median peak viral load that coincided with symptom onset. Patients with age ≥65 y had a smaller loss rate of infected cells, leading to a delayed median time to viral clearance occurring 16 d after symptom onset as compared to 13 d in younger patients (P < 10⁻⁴). In multivariate analysis, the risk factors associated with mortality were age ≥65 y, male gender, and presence of chronic pulmonary disease (hazard ratio [HR] > 2.0). Using a joint model, viral dynamics after hospital admission was an independent predictor of mortality (HR = 1.31, P < 10⁻³). Finally, we used our model to simulate the effects of effective pharmacological interventions on time to viral clearance and mortality. A treatment able to reduce viral production by 90% upon hospital admission would shorten the time to viral clearance by 2.0 and 2.9 d in patients of age <65 y and ≥65 y, respectively. Assuming that the association between viral dynamics and mortality would remain similar to that observed in our population, this could translate into a reduction of mortality from 19 to 14% in patients of age ≥65 y with risk factors. Our results show that viral dynamics is associated with mortality in hospitalized patients. Strategies aiming to reduce viral load could have an effect on mortality rate in this population.

COVID-19 transmits by droplets generated from surfaces of airway mucus during processes of respiration within hosts infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. We studied respiratory droplet generation and exhalation in human and nonhuman primate subjects with and without COVID-19 infection to explore whether SARS-CoV-2 infection, and other changes in physiological state, translate into observable evolution of numbers and sizes of exhaled respiratory droplets in healthy and diseased subjects. In our observational cohort study of the exhaled breath particles of 194 healthy human subjects, and in our experimental infection study of eight nonhuman primates infected, by aerosol, with SARS-CoV-2, we found that exhaled aerosol particles vary between subjects by three orders of magnitude, with exhaled respiratory droplet number increasing with degree of COVID-19 infection and elevated BMI-years. We observed that 18% of human subjects (35) accounted for 80% of the exhaled bioaerosol of the group (194), reflecting a superspreader distribution of bioaerosol analogous to a classical 20:80 superspreader of infection distribution. These findings suggest that quantitative assessment and control of exhaled aerosol may be critical to slowing the airborne spread of COVID-19 in the absence of an effective and widely disseminated vaccine.

Several lines of existing evidence support the possibility of airborne transmission of coronavirus disease 2019 (COVID-19). However, quantitative information on the relative importance of transmission pathways of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains limited. To evaluate the relative importance of multiple transmission routes for SARS-CoV-2, we developed a modeling framework and leveraged detailed information available from the Diamond Princess cruise ship outbreak that occurred in early 2020. We modeled 21,600 scenarios to generate a matrix of solutions across a full range of assumptions for eight unknown or uncertain epidemic and mechanistic transmission factors. A total of 132 model iterations met acceptability criteria (R² > 0.95 for modeled vs. reported cumulative daily cases and R² > 0 for daily cases). Analyzing only these successful model iterations quantifies the likely contributions of each defined mode of transmission. Mean estimates of the contributions of short-range, long-range, and fomite transmission modes to infected cases across the entire simulation period were 35%, 35%, and 30%, respectively. Mean estimates of the contributions of larger respiratory droplets and smaller respiratory aerosols were 41% and 59%, respectively. Our results demonstrate that aerosol inhalation was likely the dominant contributor to COVID-19 transmission among the passengers, even considering a conservative assumption of high ventilation rates and no air recirculation conditions for the cruise ship. Moreover, close-range and long-range transmission likely contributed similarly to disease progression aboard the ship, with fomite transmission playing a smaller role. The passenger quarantine also affected the importance of each mode, demonstrating the impacts of the interventions.

The RNA polymerase inhibitor favipiravir is currently in clinical trials as a treatment for infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), despite limited information about the molecular basis for its activity. Here we report the structure of favipiravir ribonucleoside triphosphate (favipiravir-RTP) in complex with the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) bound to a template:primer RNA duplex, determined by electron cryomicroscopy (cryoEM) to a resolution of 2.5 Å. The structure shows clear evidence for the inhibitor at the catalytic site of the enzyme, and resolves the conformation of key side chains and ions surrounding the binding pocket. Polymerase activity assays indicate that the inhibitor is weakly incorporated into the RNA primer strand, and suppresses RNA replication in the presence of natural nucleotides. The structure reveals an unusual, nonproductive binding mode of favipiravir-RTP at the catalytic site of SARS-CoV-2 RdRp, which explains its low rate of incorporation into the RNA primer strand. Together, these findings inform current and future efforts to develop polymerase inhibitors for SARS coronaviruses.

The outbreak of COVID-19 caused by SARS-CoV-2 has resulted in more than 50 million confirmed cases and over 1 million deaths worldwide as of November 2020. Currently, there are no effective antivirals approved by the Food and Drug Administration to contain this pandemic except the antiviral agent remdesivir. In addition, the trimeric spike protein on the viral surface is highly glycosylated and almost 200,000 variants with mutations at more than 1,000 positions in its 1,273 amino acid sequence were reported, posing a major challenge in the development of antibodies and vaccines. It is therefore urgently needed to have alternative and timely treatments for the disease. In this study, we used a cell-based infection assay to screen more than 3,000 agents used in humans and animals, including 2,855 small molecules and 190 traditional herbal medicines, and identified 15 active small molecules in concentrations ranging from 0.1 nM to 50 μM. Two enzymatic assays, along with molecular modeling, were then developed to confirm those targeting the virus 3CL protease and the RNA-dependent RNA polymerase. Several water extracts of herbal medicines were active in the cell-based assay and could be further developed as plant-derived anti–SARS-CoV-2 agents. Some of the active compounds identified in the screen were further tested in vivo, and it was found that mefloquine, nelfinavir, and extracts of Ganoderma lucidum (RF3), Perilla frutescens, and Mentha haplocalyx were effective in a challenge study using hamsters as disease model.

Immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection during the current pandemic remains a field of immense interest and active research worldwide. Although the severity of acute infection may depend on the intensity of innate and adaptive immunity, leading to higher morbidity and mortality, the longevity of IgG antibodies, including neutralizing activity to SARS-CoV-2, is viewed as a key correlate of immune protection. Amid reports and concern that there is a rapid decay of IgG antibody levels within 1 mo to 2 mo after acute infection, we set out to study the pattern and duration of IgG antibody response to various SARS-CoV-2 antigens in asymptomatic and symptomatic patients in a community setting. Herein, we show the correlation of IgG anti-spike protein S1 subunit, receptor binding domain, nucleocapsid, and virus neutralizing antibody titers with each other and with clinical features such as length and severity of COVID-19 illness. More importantly, using orthogonal measurements, we found the IgG titers to persist for more than 4 mo post symptom onset, implying that long-lasting immunity to COVID-19 from infection or vaccination might be observed, as seen with other coronaviruses such as SARS and Middle East respiratory syndrome.

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has emerged as a major threat to global health. Although varied SARS-CoV-2–related coronaviruses have been isolated from bats and SARS-CoV-2 may infect bat, the structural basis for SARS-CoV-2 to utilize the human receptor counterpart bat angiotensin-converting enzyme 2 (bACE2) for virus infection remains less understood. Here, we report that the SARS-CoV-2 spike protein receptor binding domain (RBD) could bind to bACE2 from Rhinolophus macrotis (bACE2-Rm) with substantially lower affinity compared with that to the human ACE2 (hACE2), and its infectivity to host cells expressing bACE2-Rm was confirmed with pseudotyped SARS-CoV-2 virus and SARS-CoV-2 wild virus. The structure of the SARS-CoV-2 RBD with the bACE2-Rm complex was determined, revealing a binding mode similar to that of hACE2. The analysis of binding details between SARS-CoV-2 RBD and bACE2-Rm revealed that the interacting network involving Y41 and E42 of bACE2-Rm showed substantial differences with that to hACE2. Bats have extensive species diversity and the residues for RBD binding in bACE2 receptor varied substantially among different bat species. Notably, the Y41H mutant, which exists in many bats, attenuates the binding capacity of bACE2-Rm, indicating the central roles of Y41 in the interaction network. These findings would benefit our understanding of the potential infection of SARS-CoV-2 in varied species of bats.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, employs two key host proteins to gain entry and replicate within cells, angiotensin-converting enzyme 2 (ACE2) and the cell surface transmembrane protease serine 2 (TMPRSS2). TMPRSS2 was first characterized as an androgen-regulated gene in the prostate. Supporting a role for sex hormones, males relative to females are disproportionately affected by COVID-19 in terms of mortality and morbidity. Several studies, including one employing a large epidemiological cohort, suggested that blocking androgen signaling is protective against COVID-19. Here, we demonstrate that androgens regulate the expression of ACE2, TMPRSS2, and androgen receptor (AR) in subsets of lung epithelial cells. AR levels are markedly elevated in males relative to females greater than 70 y of age. In males greater than 70 y old, smoking was associated with elevated levels of AR and ACE2 in lung epithelial cells. Transcriptional repression of the AR enhanceosome with AR or bromodomain and extraterminal domain (BET) antagonists inhibited SARS-CoV-2 infection in vitro. Taken together, these studies support further investigation of transcriptional inhibition of critical host factors in the treatment or prevention of COVID-19.

As the COVID-19 pandemic is spreading around the world, increasing evidence highlights the role of cardiometabolic risk factors in determining the susceptibility to the disease. The fragmented data collected during the initial emergency limited the possibility of investigating the effect of highly correlated covariates and of modeling the interplay between risk factors and medication. The present study is based on comprehensive monitoring of 576 COVID-19 patients. Different statistical approaches were applied to gain a comprehensive insight in terms of both the identification of risk factors and the analysis of dependency structure among clinical and demographic characteristics. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus enters host cells by binding to the angiotensin-converting enzyme 2 (ACE2), but whether or not renin−angiotensin−aldosterone system inhibitors (RAASi) would be beneficial to COVID-19 cases remains controversial. The survival tree approach was applied to define a multilayer risk stratification and better profile patient survival with respect to drug regimens, showing a significant protective effect of RAASi with a reduced risk of in-hospital death. Bayesian networks were estimated, to uncover complex interrelationships and confounding effects. The results confirmed the role of RAASi in reducing the risk of death in COVID-19 patients. De novo treatment with RAASi in patients hospitalized with COVID-19 should be prospectively investigated in a randomized controlled trial to ascertain the extent of risk reduction for in-hospital death in COVID-19.

Information about macromolecular structure of protein complexes and related cellular and molecular mechanisms can assist the search for vaccines and drug development processes. To obtain such structural information, we present DeepTracer, a fully automated deep learning-based method for fast de novo multichain protein complex structure determination from high-resolution cryoelectron microscopy (cryo-EM) maps. We applied DeepTracer on a previously published set of 476 raw experimental cryo-EM maps and compared the results with a current state of the art method. The residue coverage increased by over 30% using DeepTracer, and the rmsd value improved from 1.29 Å to 1.18 Å. Additionally, we applied DeepTracer on a set of 62 coronavirus-related cryo-EM maps, among them 10 with no deposited structure available in EMDataResource. We observed an average residue match of 84% with the deposited structures and an average rmsd of 0.93 Å. Additional tests with related methods further exemplify DeepTracer’s competitive accuracy and efficiency of structure modeling. DeepTracer allows for exceptionally fast computations, making it possible to trace around 60,000 residues in 350 chains within only 2 h. The web service is globally accessible at https://deeptracer.uw.edu.

With nearly every country combating the 2019 novel coronavirus (COVID-19), there is a need to understand how local environmental conditions may modify transmission. To date, quantifying seasonality of the disease has been limited by scarce data and the difficulty of isolating climatological variables from other drivers of transmission in observational studies. We combine a spatially resolved dataset of confirmed COVID-19 cases, composed of 3,235 regions across 173 countries, with local environmental conditions and a statistical approach developed to quantify causal effects of environmental conditions in observational data settings. We find that ultraviolet (UV) radiation has a statistically significant effect on daily COVID-19 growth rates: a SD increase in UV lowers the daily growth rate of COVID-19 cases by ∼

We use the synthetic control method to analyze the effect of face masks on the spread of COVID-19 in Germany. Our identification approach exploits regional variation in the point in time when wearing of face masks became mandatory in public transport and shops. Depending on the region we consider, we find that face masks reduced the number of newly registered severe acute respiratory syndrome coronavirus 2 infections between 15% and 75% over a period of 20 days after their mandatory introduction. Assessing the credibility of the various estimates, we conclude that face masks reduce the daily growth rate of reported infections by around 47%.

The COVID-19 pandemic is caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and has spread worldwide, with millions of cases and more than 1 million deaths to date. The gravity of the situation mandates accelerated efforts to identify safe and effective vaccines. Here, we generated measles virus (MeV)-based vaccine candidates expressing the SARS-CoV-2 spike glycoprotein (S). Insertion of the full-length S protein gene in two different MeV genomic positions resulted in modulated S protein expression. The variant with lower S protein expression levels was genetically stable and induced high levels of effective Th1-biased antibody and T cell responses in mice after two immunizations. In addition to neutralizing IgG antibody responses in a protective range, multifunctional CD8⁺ and CD4⁺ T cell responses with S protein-specific killing activity were detected. Upon challenge using a mouse-adapted SARS-CoV-2, virus loads in vaccinated mice were significantly lower, while vaccinated Syrian hamsters revealed protection in a harsh challenge setup using an early-passage human patient isolate. These results are highly encouraging and support further development of MeV-based COVID-19 vaccines.