key: cord-0948913-0xfvilbf
authors: Ye, Chengjin; Chiem, Kevin; Park, Jun-Gyu; Silvas, Jesus A.; Morales Vasquez, Desarey; Sourimant, Julien; Lin, Michelle J.; Greninger, Alexander L.; Plemper, Richard K.; Torrelles, Jordi B.; Kobie, James J.; Walter, Mark R.; de la Torre, Juan Carlos; Martinez-Sobrido, Luis
title: Analysis of SARS-CoV-2 infection dynamic in vivo using reporter-expressing viruses
date: 2021-10-12
journal: Proc Natl Acad Sci U S A
DOI: 10.1073/pnas.2111593118
sha: 08eeb20d3aa819f10326e7566638326cb059c425
doc_id: 948913
cord_uid: 0xfvilbf

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the current COVID-19 pandemic, is one of the biggest threats to public health. However, the dynamic of SARS-CoV-2 infection remains poorly understood. Replication-competent recombinant viruses expressing reporter genes provide valuable tools to investigate viral infection. Low levels of reporter gene expressed from previous reporter-expressing recombinant (r)SARS-CoV-2 in the locus of the open reading frame (ORF)7a protein have jeopardized their use to monitor the dynamic of SARS-CoV-2 infection in vitro or in vivo. Here, we report an alternative strategy where reporter genes were placed upstream of the highly expressed viral nucleocapsid (N) gene followed by a porcine tescherovirus (PTV-1) 2A proteolytic cleavage site. The higher levels of reporter expression using this strategy resulted in efficient visualization of rSARS-CoV-2 in infected cultured cells and excised lungs or whole organism of infected K18 human angiotensin converting enzyme 2 (hACE2) transgenic mice. Importantly, real-time viral infection was readily tracked using a noninvasive in vivo imaging system and allowed us to rapidly identify antibodies which are able to neutralize SARS-CoV-2 infection in vivo. Notably, these reporter-expressing rSARS-CoV-2, in which a viral gene was not deleted, not only retained wild-type (WT) virus-like pathogenicity in vivo but also exhibited high stability in vitro and in vivo, supporting their use to investigate viral infection, dissemination, pathogenesis, and therapeutic interventions for the treatment of SARS-CoV-2 in vivo.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the current COVID-19 pandemic, is one of the biggest threats to public health. However, the dynamic of SARS-CoV-2 infection remains poorly understood. Replication-competent recombinant viruses expressing reporter genes provide valuable tools to investigate viral infection. Low levels of reporter gene expressed from previous reporter-expressing recombinant (r)SARS-CoV-2 in the locus of the open reading frame (ORF)7a protein have jeopardized their use to monitor the dynamic of SARS-CoV-2 infection in vitro or in vivo. Here, we report an alternative strategy where reporter genes were placed upstream of the highly expressed viral nucleocapsid (N) gene followed by a porcine tescherovirus (PTV-1) 2A proteolytic cleavage site. The higher levels of reporter expression using this strategy resulted in efficient visualization of rSARS-CoV-2 in infected cultured cells and excised lungs or whole organism of infected K18 human angiotensin converting enzyme 2 (hACE2) transgenic mice. Importantly, real-time viral infection was readily tracked using a noninvasive in vivo imaging system and allowed us to rapidly identify antibodies which are able to neutralize SARS-CoV-2 infection in vivo. Notably, these reporter-expressing rSARS-CoV-2, in which a viral gene was not deleted, not only retained wild-type (WT) virus-like pathogenicity in vivo but also exhibited high stability in vitro and in vivo, supporting their use to investigate viral infection, dissemination, pathogenesis, and therapeutic interventions for the treatment of SARS-CoV-2 in vivo.

COVID-19 | SARS-CoV-2 | in vivo imaging | reporter viruses C oronaviruses (CoVs) are enveloped, single-stranded, positivesense RNA viruses that belong to the Coronaviridae family that can cause mild to severe respiratory infections in humans (1) . Two CoVs have been associated with severe respiratory syndrome in the past two decades: severe acute respiratory syndrome CoV (SARS-CoV) in 2002-2003 and Middle East respiratory syndrome CoV (MERS-CoV) in 2012 to the present (2) . Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in the Chinese city of Wuhan in December 2019 and is the causative agent of the COVID-19 pandemic (3, 4) . As of June 2021, SARS-CoV-2 has been reported to be responsible for over 150 million human infection cases and more than 3 million deaths around the World (https://covid19.who.int/).

Like SARS-CoV and MERS-CoV, SARS-CoV-2 mainly replicates in the upper (nasal turbinate) and lower (lungs) respiratory tract, resulting, in some cases, in fatal respiratory illness (5, 6) . However, the intrahost dissemination and pathogenesis of SARS-CoV-2 are not well understood. Several animal models of SARS-CoV-2 infection have been established and have already provided very valuable information to understand the mechanism of tissue and cell tropism, replication, and pathogenesis (7) (8) (9) (10) (11) . However, assessing the presence of SARS-CoV-2 in infected animals, organs, or tissues has required collection and processing of samples upon euthanasia, which complicates studies examining the longitudinal dynamic of a viral infection within an infected host. Recombinant (r)SARS-CoV-2 expressing reporter genes could overcome this problem and allow tracking of viral infection in vivo and in real time by monitoring the expression of the reporter gene. We and others have documented the feasibility of generating reporter-expressing rSARS-CoV-2 using a reverse genetic system (12, 13) . These rSARS-CoV-2 have been genetically engineered to express the reporter gene by substituting the viral open reading frame (ORF) 7a protein with the reporter gene of interest, an experimental approach first employed to generate reporter-expressing rSARS-CoV (14) . Despite these reporterexpressing rSARS-CoV-2 showing plaque phenotype, replication, and growth kinetics comparable to those of wild-type virus (rSARS-CoV-2/WT) in vitro (12, 13, 15) , it is unclear whether the reporter-expressing rSARS-CoV-2 lacking ORF7a recapitulate viral pathogenicity in vivo and whether reporter gene expression levels could be efficiently tracked ex vivo using tissues or organs from infected animals, or in a whole organism in vivo.

In this study, we cloned fluorescent (Venus) and luciferase (Nano luciferase, Nluc) reporter genes upstream of the SARS-CoV-2

To date, due to the insufficient expression level of reporter genes from previous recombinant (r)SARS-CoV-2, in which the viral open reading frame (ORF) 7a gene was replaced by reporter genes, tracking the SARS-CoV-2 infection dynamic has been challenging. Here, we engineered rSARS-CoV-2 expressing fluorescent (Venus) or luciferase (Nano luciferase) reporter genes from the viral nucleocapsid (N) locus, without deleting any viral gene. These novel reporter-expressing rSARS-CoV-2, which give a higher expression level of reporter genes, allow us to monitor SARS-CoV-2 replication dynamic both in vitro and in vivo. These new reporter-expressing rSARS-CoVs-2 represent an excellent option to assess viral replication, tropisms, and pathogenicity as well as the rapid in vivo evaluation of effective countermeasures for the treatment of SARS-CoV-2 infection. nucleocapsid (N) gene separated by the porcine tescherovirus (PTV-1) 2A proteolytic cleavage site to generate new reporter-expressing rSARS-CoV-2 without the deletion of the ORF7a protein. In vitro, rSARS-CoV-2 expressing reporter genes from the viral N locus replicated and made viral plaques similar to those of rSARS-CoV-2/ WT. Reporter-expressing rSARS-CoV-2 generated using this 2A strategy expressed higher levels of reporter gene expression compared to those rSARS-CoV-2 generated by substituting the viral ORF7a protein with the reporter gene of interest. Importantly, rSARS-CoV-2/Venus-2A and rSARS-CoV-2/Nluc-2A showed rSARS-CoV-2/WT−like pathogenicity in vivo. The higher level of Venus expression from rSARS-CoV-2/Venus-2A allowed us to detect viral infection in the lungs of infected K18 human angiotensin converting enzyme 2 (hACE2) transgenic mice using an in vivo imaging system (IVIS). Moreover, Venus expression from rSARS-CoV-2/Venus-2A was stable up to seven passages in vitro in cultured Vero E6 cells and in vivo up to day 6 postinfection. Importantly, levels of Venus expression correlated well with viral titers detected in the lungs, demonstrating the feasibility of using Venus expression as a valid surrogate marker to evaluate SARS-CoV-2 infection. Using rSARS-CoV-2/Nluc-2A, we were able to track the dynamic of viral infection in real time and longitudinally assess SARS-CoV-2 infection in vivo. Finally, we testified to the feasibility of using the rSARS-CoV-2/Nluc-2A to rapidly and accurately identify antibodies that neutralize viral infection in vivo.

Our data demonstrate that these next-generation rSARS-CoV-2 expressing reporter genes we have generated can be used to easily monitor viral infection in cultured cells and in validated animal models of infection. Importantly, our new rSARS-CoV-2/Venus-2A or rSARS-CoV-2/Nluc-2A retain similar virulence to that of rSARS-CoV-2/WT in K18 hACE2 transgenic mice and can be used to investigate viral replication, tropism, and viral dissemination and pathogenesis in vivo and to rapidly identify therapeutics for the treatment of SARS-CoV-2 infection and associated COVID-19 disease.

Generation of rSARS-CoV-2 Expressing Venus. We have recently described the generation and characterization of rSARS-CoV-2 where a reporter gene of interest replaced the viral ORF7a protein (13) . However, these rSARS-CoV-2 showed low levels of reporter gene expression during viral infection. To increase expression levels of reporter gene during SARS-CoV-2 infection and avoid the deletion of ORF7a protein, we implemented a strategy we previously used to generate recombinant influenza viruses and mammarenaviruses (16, 17) . In this approach, the sequences of the reporter Venus fluorescent protein and the porcine teschovirus 1 (PTV-1) 2A self-cleaving peptide (18) were cloned upstream of the SARS-CoV-2 N gene (Fig. 1A) , and the sequences harboring the fusion Venus-2A-N were cloned into the previously described bacterial artificial chromosome (BAC) containing the entire SARS-CoV-2 genome (19, 20) . We rescued rSARS-CoV-2/Venus-2A virus according to our previously described protocol (19, 20) . Vero E6 cells infected with the tissue culture supernatant (TCS) from Vero E6 cells transfected with the BAC containing the viral genome of rSARS-CoV-2/Venus-2A resulted in the expression of Venus in the same cells expressing the viral N protein (Fig. 1B) . Venus was readily detected in whole cell lysate from rSARS-CoV-2/Venus-2A−, but not rSARS-CoV-2/WT−infected Vero E6 cells, while the viral N protein was detected in lysate obtained from both rSARS-CoV-2/Venus-2A− and rSARS-CoV-2/WT−infected Vero E6 cells (Fig. 1C) . We confirmed the genetic identity of rSARS-CoV-2/Venus-2A using RT-PCR to amplify the Venus sequence and the entire sequence between ORF8 and N. The Venus fragment was amplified from cells infected with rSARS-CoV-2/ Venus-2A, while the fragment between ORF8 and N was detected in cells infected with either rSARS-CoV-2/WT or rSARS-CoV-2/ Venus-2A. As predicted, the amplified fragment from rSARS-CoV-2/Venus-2A−infected cells had a higher molecular size than the one obtained from rSARS-CoV-2/WT−infected cells (Fig. 1D ).

In Vitro Characterization of rSARS-CoV-2/Venus-2A. We next assessed the fitness of rSARS-CoV-2/Venus-2A in Vero E6 cells by evaluating its growth kinetics, plaque phenotype, and reporter expression and compared them to our previously described rSARS-CoV-2/Δ7a-Venus and rSARS-CoV-2/WT. The rSARS-CoV-2/Venus-2A showed similar growth kinetics to those of rSARS-CoV-2/WT or rSARS-CoV-2/Δ7a-Venus ( Fig. 2A) . Likewise, rSARS-CoV-2/Venus-2A, rSARS-CoV-2/WT, and rSARS-CoV-2/Δ7a-Venus exhibited similar plaque phenotypes (Fig. 2B) . Notably, plaques formed by rSARS-CoV-2/Venus-2A, but not those of rSARS-CoV-2/Δ7a-Venus, could be readily detected using a fluorescent imaging system, most likely because of the higher levels of Venus expressed from the locus of the viral N than those from the locus of the viral ORF7a (Fig. 2B ). We next examined Venus expression in Vero E6 cells infected (MOI [multiplicity of infection] 0.001) with either rSARS-CoV-2/Venus-2A or rSARS-CoV-2/Δ7a-Venus. Both infections showed peak Venus expression at 48 h postinfection (hpi), but Venus expression levels from rSARS-CoV-2/Venus-2A−infected Vero E6 cells were higher than those infected with rSARS-CoV-2/Δ7a-Venus (Fig. 2C ). To confirm this, we prepared whole cell lysates from mock, rSARS-CoV-2/Venus-2A−, and rSARS-CoV-2/Δ7a-Venus−infected Vero E6 cells and analyzed them by Western blot. Venus expression levels in rSARS-CoV-2/Venus-2A−infected Vero E6 cells were greater than those of rSARS-CoV-2/ Δ7a-Venus−infected cells at all time points, whereas expression levels of N protein were comparable at all time points in Vero E6 cells infected with rSARS-CoV-2/Venus-2A and rSARS-CoV-2/ Δ7a-Venus (Fig. 2D ).

In Vivo Characterization of rSARS-CoV-2/Venus-2A. K18 transgenic mice expressing hACE2 have been shown to be a good animal model of SARS-CoV-2 infection (8, 21) . We therefore examined whether SARS-CoV-2 infection could be tracked ex vivo using Venus expression. To that end, K18 hACE2 transgenic mice were infected intranasally with 10 5 plaque-forming units (pfu) of rSARS-CoV-2/Venus-2A, rSARS-CoV-2/Δ7a-Venus, or rSARS-CoV-2/WT (Fig. 3 ). Mice were killed at 1, 2, 4, and 6 d postinfection (dpi), and their lungs were excised and imaged ex vivo using an IVIS (AMI spectrum). Venus expression was readily detected in all lungs obtained from mice infected with rSARS-CoV-2/Venus-2A but not from those infected with rSARS-CoV-2/Δ7a-Venus or rSARS-CoV-2/WT (Fig. 3A) . Quantitative analyses showed that Venus intensity peaks at 2 dpi and decreases over the course of infection in the lungs of infected mice (Fig. 3B) . Nevertheless, gross lesions on the lung surface of mice infected with rSARS-CoV-2/ Venus-2A were comparable to those observed in rSARS-CoV-2/ WT− or rSARS-CoV-2/Δ7a-Venus−infected mice ( Fig. 3 C and D) .

In addition, infection with rSARS-CoV-2/Venus-2A and rSARS-CoV-2/Δ7a-Venus resulted in comparable viral titers to those observed in K18 hACE2 transgenic mice infected with rASRS-CoV-2/ WT in all organs at all times pi (Fig. 3E) , suggesting that the undetectable Venus expression in the lungs of K18 hACE2 mice infected with rSARS-CoV-2/Δ7a-Venus is unlikely due to lower levels of viral replication in vivo. Notably, we observed a correlation between virus replication and fluorescence intensity in the lungs (Fig. 3F) .

To investigate the pathogenicity of rSARS-CoV-2/Venus-2A, we infected (10 5 pfu) K18 hACE2 transgenic mice intranasally with rSARS-CoV-2/Venus-2A, rSARS-CoV-2/Δ7a-Venus, or rSARS-CoV-2/WT, and monitored changes in body weight and survival rate for 12 d after viral infection. All infected mice showed significant body weight loss starting from 4 dpi (Fig. 3G ).

Mice infected with rSARS-CoV-2/Δ7a-Venus succumbed to viral infection by 8 dpi, and mice infected with rSARS-CoV-2/Venus-2A or rSARS-CoV-2/WT succumbed to infection by 9 dpi (Fig. 3H ).

In Vitro and In Vivo Stability of rSARS-CoV-2/Venus-2A. Because reporter-expressing rSARS-CoV-2 applications necessitate genetic and phenotypic stability, we evaluated the stability of rSARS-CoV-2/Venus-2A in vitro and in vivo. To that end, we passaged the rSARS-CoV-2/Venus-2A seven times in Vero E6 cells, and TCSs from selected passages were subjected to the analysis of plaque assay using a fluorescent imaging system and immunostaining with the 1C7C7 N protein monoclonal antibody (MAb) (Fig. 4 A, Left) . Viruses present in TCS from P1, P3, P5, and P7 retained 100% Venus expression (Fig. 4 A, Right) . For the in vivo stability evaluation of rSARS-CoV-2/Venus-2A, lung homogenates from K18 hACE2 transgenic mice at 1, 2, 4, and 6 dpi infected with rSARS-CoV-2/Venus-2A were subjected to plaque assay using a fluorescent imaging system and immunostaining with the 1C7C7 N protein MAb (Fig. 4 B, Left) . Lung homogenates from 1 and 2 dpi retained 100% Venus expression, and homogenates from 4 and 6 dpi retained 99% and 98%, respectively, Venus expression ( Fig. 4 B, Right) .

Visualization of SARS-CoV-2 Replication Dynamic In Vivo. While fluorescent Venus expression allowed us to conduct ex vivo imaging of lungs from SARS-CoV-2−infected mice, it did not allow us to track viral infection in the entire mouse using IVIS. To circumvent this problem, we engineered an Nluc-expressing rSARS-CoV-2 using our reverse genetic system and the same 2A strategy (Fig. 5A ). This rSARS-CoV-2/Nluc-2A exhibited a similar plaque phenotype and comparable growth kinetics in Vero E6 cells to the rSARS-CoV-2/WT and our previously described rSARS-CoV-2/Δ7a-Nluc ( Fig. 5 B and C) . Notably, Nluc expression levels were increased by more than 30-fold in Vero E6 cells infected with rSARS-CoV-2/Nluc-2A compared with cells infected with rSARS-CoV-2/Δ7a-Nluc at 72 hpi (Fig. 5D ).

Since the rSARS-CoV-2/Nluc-2A expressed significantly higher levels of Nluc than those of our previously described rSARS-CoV-2/Δ7a-Nluc in vitro, we evaluated whether SARS-CoV-2 infection could be tracked in vivo using Nluc expression directed by rSARS-CoV-2/Nluc-2A. To that end, we infected (10 5 pfu) K18 hACE2 transgenic mice intranasally with rSARS-CoV-2/Nluc-2A or rSARS-CoV-2/WT (Fig. 6 ). Mice were anesthetized, retroorbitally injected with Nluc substrate, and then imaged under an IVIS at 1, 2, 4, and 6 dpi. Nluc expression was readily detected in mice infected with rSARS-CoV-2/Nluc-2A but not those infected with rSARS-CoV-2/WT (Fig. 6A) , as previously shown in vitro (Fig. 5D) . Quantitative analyses showed that Nluc intensity continued increasing at later dpi (Fig. 6B ). Gross lesions on the lung surface of mice infected with rSARS-CoV-2/Nluc-2A were comparable to those in the WT rSARS-CoV-2−infected mice (Fig. 6 C and D) . Importantly, viral titers detected in the rSARS-CoV-2/Nluc-2A− infected mice were comparable to those infected with rSARS-CoV-2/WT in all organs tested at different dpi (Fig. 6E) , even though Nluc activity was only detected in the organs from rSARS-CoV-2/Nluc-2A−infected mice (Fig. 6F) . Unexpectedly, the Nluc intensity did not correlate with the viral titers in the lungs at 4 and 6 dpi (Fig. 6G) , indicating that Nluc accumulates in vivo. Meanwhile, we also compared the pathogenicity of rSARS-CoV-2/ Nluc-2A and rSARS-CoV-2/WT in K18 hACE2 transgenic mice. Both rSARS-CoV-2/Nluc-2A and rSARS-CoV-2/WT showed similar virulence (Fig. 6 H and I) .

Vivo. Although several vaccines against COVID-19 have been already authorized for emergency use by US Food and Drug Administration (FDA), identification and characterization of SARS-CoV-2 neutralizing Abs (NAbs) represent a valuable therapeutic option to counteract the putative emergence of variants of concern (VoC). Currently, most NAb screenings are performed in tissue cultured cells rather than in vivo, which is heavily reliant upon viral titration of animal organs, a process that is time and labor intensive. We therefore investigated whether the use of rSARS-CoV-2/Nluc-2A could expedite the screening process and facilitate the investigation of how NAb affects the kinetics of virus infection. To that end, we treated K18 hACE2 transgenic mice prophylactically with 1212C2, a previously described SARS-CoV-2 NAb (22), for 12 h, and then infected them (10 5 pfu) with rSARS-CoV-2/Nluc-2A. Noninvasive longitudinal imaging of the mice revealed that 1212C2-treated mice dramatically restricted rSARS-CoV-2/Nluc-2A multiplication, since no significant Nluc signal was detected at any time point examined (Fig. 7 A  and B) . These results were also supported by the observation of reduced lung surface lesions (Fig. 7 C and D) , Nluc expression (Fig. 7E) , and virus titers (Fig. 7F) in the nasal turbinate, lungs, and brains of 1212C2-treated mice. Moreover, 1212C2 was able to protect mice from clinical symptoms of rSARS-CoV-2/Nluc-2A infection as determined by changes in body weight (Fig. 7G ) and survival rate (Fig. 7H ).

In this study, we report a strategy to generate replicationcompetent reporter-expressing (e.g., Venus or Nluc) rSARS-CoV-2 using the well-documented BAC-based reverse genetic system (20) . Our rSARS-CoV-2/Venus-2A and rSARS-CoV-2/ Nluc-2A both exhibited rSARS-CoV-2/WT−like growth properties in vitro and in vivo without displaying attenuation and allowed us to monitor virus infection ex vivo in the lungs of infected mice (rSARS-CoV-2/Venus-2A) and the dynamic of viral replication in the entire mouse (rSARS-CoV-2/Nluc-2A) using noninvasive longitudinal in vivo imaging. Importantly, we demonstrate the feasibility of using rSARS-CoV-2/Nluc-2A to rapidly identify prophylactics and/or therapeutics in vivo.

We and others have previously demonstrated the feasibility of using this 2A approach to generate recombinant viruses expressing reporter genes fused to a viral protein (16) (17) (18) 23) . To generate these novel reporter rSARS-CoV-2 expressing a higher level of reporter gene, we placed the PTV-1 2A self-cleaving peptide between the reporter gene and the viral N gene (Fig. 1) . This results in the expression of a polyprotein that is posttranslationally cleaved at the 2A site leading to the individual expression of the reporter gene and the viral N protein (24) , and the expression of Venus in rSARS-CoV-2/Venus-2A was extremely increased in vitro and in vivo (Figs. 2 and 3, respectively). The rationale of cloning the reporter gene fused to the viral N gene to increase reporter gene expression was based on the N protein being one of the most abundant structural proteins produced during SARS-CoV-2 and other CoVs infections (25) (26) (27) . Importantly, this 2A approach does not remove any viral gene from the viral genome. Although recent data from our laboratory (28) suggest that ORF7a is not essential for SARS-CoV-2 replication in vitro and in vivo (29) , it is still largely unknown whether the lack of ORF7a affects some unknown aspects of SARS-CoV-2 pathogenesis. In addition, this also increases the instability concern regarding the reporter gene, as shown in a recent study where a reporter gene fused into the C terminus of ORF7a was not stable (30) . Contrarily, our rSARS-CoV-2/Venus-2A exhibited 100% stability after seven passages in cultured cells and retained 98% stability in vivo at 6 dpi (Fig. 4) .

The increased expression level of reporter gene also facilitates the use of bioluminescence imaging of the entire infected mouse. The rSARS-CoV-2/Nluc-2A expressed ∼30-fold higher levels of Nluc than rSARS-CoV-2/Δ7a-Nluc in cultured cells (Fig. 5) , which allowed us to track the viral infection as early as 1 dpi (Fig. 6A) . Another advantage is that the rSARS-CoV-2/Nluc-2A could provide real-time and longitudinal information in a noninvasive manner during infection rather than providing a static "snapshot" using the traditional necropsy and titration of tissues or organs from an infected animal. Notably, Nluc activity correlated well with the titers detected in the lungs at 1 and 2 dpi (Fig. 6 E and F) . In contrast, viral titers in lungs were decreased at 4 and 6 dpi, whereas we still detected high levels of Nluc activity. This may reflect Nluc stability (31) , which can lead to the accumulation and gradual increase of Nluc signal during the course of viral infection. This explanation was also supported by the in vitro infection data in which the viral titers were declining and yet the Nluc activity was gradually increasing (Fig. 5D ).

NAb represents a promising prophylactic and/or therapeutic treatment against SARS-CoV-2, particularly for individuals infected with newly identified VoC. However, evaluation of SARS-CoV-2 NAb in vivo relies on the necropsy and viral titration of tissues and/or organs from infected animals. To overcome this limitation, we established a rapid method based on a noninvasive measurement of Nluc expression. By using this IVIS, NAbs could be easily identified as early as 1 dpi and in a relatively high-throughput method. This was further supported by Nluc activity, viral titration, body weight changes, and survival rate (Fig. 7) , indicating our strategy provides a rapid in vivo screening method to identify NAbs against SARS-CoV-2.

In the present work, we have documented a strategy to generate replication-competent reporter rSARS-CoV-2 expressing higher levels of reporter gene than those previously described, by substituting the viral ORF7a protein with the report gene. This strategy does not eliminate any viral gene, and these reporterexpressing rSARS-CoV-2 were genetically stable and replicated as efficiently as rSARS-CoV-2/WT both in vitro and in vivo, with comparable pathogenicity in K18 hACE2 transgenic mice. Notably, the robust levels of reporter gene expression of these reporter-expressing rSARS-CoV-2 represent an excellent option to study viral pathogenesis, tissue tropism, and replication kinetics of SARS-CoV-2, including recently identified VoC.

Biosafety. All the in vitro and in vivo experiments with infectious rSARS-CoV-2 were conducted in appropriate biosafety level (BSL) 3 and animal BSL3 (ABSL3) laboratories, respectively, at Texas Biomedical Research Institute (Texas Biomed). Experiments were approved by the Texas Biomed Institutional Biosafety and Animal Care and Use (IACUC) committees.

Cells and Viruses. African green monkey kidney epithelial cells (Vero E6, CRL-1586) were obtained from the American Type Culture Collection and maintained in Dulbecco's modified Eagle medium (DMEM) supplemented with 5% (vol/vol) fetal bovine serum (FBS, VWR) and 1% penicillin−streptomycin−glutamine (PSG) solution (Corning).

The rSARS-CoV-2 were generated based on the backbone of the USA-WA1/2020 strain using a previously described BAC-based reverse genetics system (19, 20) .

Rescue of rSARS-CoV-2. Virus rescue experiments were performed as previously described (32) . Briefly, confluent monolayers of Vero E6 cells (10 6 cells per well, six-well plates, triplicates) were transfected with 4.0 μg per well of SARS-CoV-2 BAC using Lipofectamine 2000. After 24 h, transfection media was exchanged for postinfection media (DMEM supplemented with 2% FBS and 1% PSG), and cells were split and seeded into T75 flasks 72 h posttransfection. After incubation for another 72 h, TCSs were collected, labeled as P0, and stored at −80°C. After being titrated, the P0 virus was used to infect fresh Vero E6 cells at MOI 0.0001 for 72 h to generate P1 stocks. P1 viral stocks were aliquoted and stored at −80°C until being used.

Western Blot. Whole cell lysates, sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), and Western blotting were performed as previously described (16) . Briefly, cells were lysed in passive lysis buffer (Promega) at 4°C for 30 min, followed by centrifugation at 12,000 × g at 4°C for another 30 min. Equivalent amounts of cell lysates were subjected to 12% SDS-PAGE and transferred to nitrocellulose membranes. After blocking with 5% bovine serum albumin in PBS containing 0.1% Tween 20 at room temperature for 1 h, the membranes were incubated with the indicated primary antibodies at 4°C overnight, followed by horseradish peroxidase-conjugated secondary antibody incubation at 37°C for 1 h. β-Actin was used as the loading control. Membranes were developed with ECL detection reagent (Thermo Fisher Scientific) in the ChemiDoc MP Imaging System (Bio-Rad).

Immunofluorescence Assay. Vero E6 cells (10 6 cells per well, six-well plate format, triplicates) were mock inoculated or inoculated with the TCS (1 mL per well) of Vero E6 cells transfected with the BACs collected at 72 h posttransfection. At 48 h postinoculation, cells were fixed with 10% formaldehyde solution at 4°C overnight and permeabilized using 0.5% (vol/vol) Triton X-100 in PBS for 15 min at room temperature. Then, cells were incubated overnight with 1 μg/mL of a SARS-CoV cross-reactive N MAb (1C7C7) at 4°C, washed stringently with PBS, and stained with a fluorescein isothiocyanate−labeled goat anti-mouse IgG (1:200). Finally, cells were visualized and imaged under an EVOS fluorescent microscope (Thermo Fisher Scientific).

RT-PCR. Total RNA from virus-infected (MOI 0.01) Vero E6 cells (10 6 cells per well, six-well plate format) was extracted with TRIzol Reagent (Thermo Fisher Scientific) according to the manufacturer's instructions. RT-PCR amplification of the viral genome spanning nucleotides 27,895 to 29,534 (according to the SARS-CoV-2 USA-WA1/2020 viral genome) was performed using Super Script II Reverse transcriptase (Thermo Fisher Scientific) and Expanded High Fidelity PCR System (Sigma Aldrich). The amplified DNA products were separated on 0.7% agarose gel.

Plaque Assay and Immunostaining. Confluent monolayers of Vero E6 cells (10 6 cells per well, six-well plate format, triplicates) were infected with serial viral dilutions for 1 h at 37°C. After viral adsorption, cells were overlaid with pi media containing 1% low melting agar and incubated at 37°C. At 72 hpi, cells were fixed overnight with 10% formaldehyde solution. For visualization of Venus, plates were photographed under a ChemiDoc MP Imaging System. For immunostaining, cells were permeabilized with 0.5% (vol/vol) Triton X-100 in PBS for 15 min at room temperature and immunostained using the SARS-CoV cross-reactive N protein 1C7C7 MAb (1 μg/mL) and the Vectastain ABC kit (Vector Laboratories), following the manufacturers' instruction. After immunostaining, plates were scanned and photographed using a ChemiDoc MP Imaging System.

Virus Growth Kinetics. Confluent monolayers of Vero E6 cells (six-well format, 10 6 cells per well, triplicates) were mock infected or infected (MOI 0.01) with rSARS-CoV-2/WT, rSARS-CoV-2/Δ7a-Venus, rSARS-CoV-2/Venus-2A, or rSARS-CoV-2/Nluc-2A. After 1 h of virus adsorption at 37°C, cells were washed with chilled PBS and overlaid with 3 mL of pi medium and incubated at 37°C. At the indicated times (12, 24, 48 , and 72 hpi), viral titers in the TCSs were determined by plaque assay (33) . Presence of Nluc in the TCSs from mock and rSARS-CoV-2/WT−, rSARS-CoV-2/Δ7a-Nluc−, or rSARS-CoV-2/Nluc-2A− infected cells was quantified using Nano-Glo Luciferase Assay System (Promega) following the manufacturers' specification.

Animal Experiments. All animal protocols were approved by Texas Biomed IACUC (1718MU). Five-week-old female K18 hACE2 transgenic mice were purchased from The Jackson Laboratory and maintained in the animal facility at Texas Biomed under specific pathogen-free conditions. For virus infection, mice were anesthetized following gaseous sedation in an isoflurane chamber and inoculated intranasally with a dose of 10 5 pfu per mouse.

For ex vivo imaging of lungs, mice were humanely killed at 1, 2, 4, and 6 dpi to collect lungs. Fluorescent images of lungs were photographed using an IVIS (AMI HTX), and the brightfield images of lungs were taken using an iPhone 6s (Apple). Nasal turbinate, lungs, and brains from mock or infected animals were homogenized in 1 mL of PBS for 20 s at 7,000 rpm using a Precellys tissue homogenizer (Bertin Instruments). Tissue homogenates were centrifuged at 12,000 × g (4°C) for 5 min, and supernatants were collected and titrated by plaque assay and immunostaining as previously described.

For in vivo bioluminescence imaging, mice were anesthetized with isoflurane, injected retroorbitally with 100 μL of Nano-Glo luciferase substrate (Promega), and immediately imaged. The bioluminescence data acquisition and analysis were performed using the Aura program (AMI spectrum). Flux measurements were acquired from the region of interest. The scale used is included in each figure. Immediately after imaging, nasal turbinate, lungs, and brains were collected and homogenized in 1 mL of PBS. Supernatants were collected, and presence of virus was determined by plaque assay and immunostaining, as described above. Nluc activity in the TCSs was determined under a multiplate reader (BioTek Instruments, Inc.) as above.

For the body weight and survival studies, 5-week-old female K18 hACE2 transgenic mice were infected intranasally with 10 5 pfu per animal following gaseous sedation in an isoflurane chamber. After infection, mice were monitored daily for morbidity (body weight) and mortality (survival rate) for 12 d.

Mice showing a loss of more than 25% of their initial body weight were defined as reaching the experimental end point and humanely killed, and the survival curves were plotted according to the method of Kaplan−Meier (34) .

Statistical Analysis. All data are presented as mean ± SD for each group and analyzed by SPSS13.0 (IBM). A P value of less than 0.05 (P < 0.05) was considered statistically significant.

Data Availability. All study data are included in the article. All primer sequences used for RT-PCR are available upon request. The plasmids described in this study are available at the following website: https://www.txbiomed. org/services-2/reverse-genetics-plasmids.

Characteristics of SARS-CoV-2 and COVID-19

SARS and MERS: Recent insights into emerging coronaviruses

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modelling study

Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding

Pathological findings of COVID-19 associated with acute respiratory distress syndrome

Clinical features of familial clustering in patients infected with 2019 novel coronavirus in Wuhan, China

Syrian hamsters as a small animal model for SARS-CoV-2 infection and countermeasure development

Lethality of SARS-CoV-2 infection in K18 human angiotensinconverting enzyme 2 transgenic mice

A mouse model of SARS-CoV-2 infection and pathogenesis

Responses to acute infection with SARS-CoV-2 in the lungs of rhesus macaques, baboons and marmosets

Infection and rapid transmission of SARS-CoV-2 in ferrets

An infectious cDNA clone of SARS-CoV-2

Generation and characterization of recombinant SARS-CoV-2 expressing reporter genes

Severe acute respiratory syndrome coronavirus infection of human ciliated airway epithelia: Role of ciliated cells in viral spread in the conducting airways of the lungs

A nanoluciferase SARS-CoV-2 for rapid neutralization testing and screening of anti-infective drugs for COVID-19

Development of reverse genetics for the prototype New World mammarenavirus Tacaribe virus

A novel fluorescent and bioluminescent bireporter influenza A virus to evaluate viral infections

Recombinant Lassa virus expressing green fluorescent protein as a tool for high-throughput drug screens and neutralizing antibody assays

Generation of recombinant SARS-CoV-2 using a bacterial artificial chromosome

Rescue of SARS-CoV-2 from a single bacterial artificial chromosome

COVID-19 treatments and pathogenesis including anosmia in K18-hACE2 mice

Therapeutic activity of an inhaled potent SARS-CoV-2 neutralizing human monoclonal antibody in hamsters

Analysis of in vivo dynamics of influenza virus infection in mice using a GFP reporter virus

2A to the fore -Research, technology and applications

Quantification of individual subgenomic mRNA species during replication of the coronavirus transmissible gastroenteritis virus

SARS-CoV-2 reverse genetics reveals a variable infection gradient in the respiratory tract

Reverse genetics with a full-length infectious cDNA of the Middle East respiratory syndrome coronavirus

Contribution of SARS-CoV-2 accessory proteins to viral pathogenicity in K18 hACE2 transgenic mice

Identification of multiple large deletions in ORF7a resulting in inframe gene fusions in clinical SARS-CoV-2 isolates

A plasmid DNA-launched SARS-CoV-2 reverse genetics system and coronavirus toolkit for COVID-19 research

Engineered luciferase reporter from a deep sea shrimp utilizing a novel imidazopyrazinone substrate

Rescue of recombinant Zika virus from a bacterial artificial chromosome cDNA clone

Influenza A virus attenuation by codon deoptimization of the NS gene for vaccine development

Logistic-regression, survival analysis, and the Kaplan-Meier curve

ACKNOWLEDGMENTS. We thank Dr. Thomas Moran at The Icahn School of Medicine at Mount Sinai for providing the SARS-CoV cross-reactive 1C7C7 N protein MAb.