key: cord-0835009-km44lv0j authors: Chiem, Kevin; Morales Vasquez, Desarey; Park, Jun-Gyu; Platt, Roy Neal; Anderson, Tim; Walter, Mark R.; Kobie, James J.; Ye, Chengjin; Martinez-Sobrido, Luis title: Generation and Characterization of Recombinant SARS-CoV-2 Expressing Reporter Genes date: 2021-03-10 journal: J Virol DOI: 10.1128/jvi.02209-20 sha: abc79ef57cc4e76c340a1b267e1f7844631c0645 doc_id: 835009 cord_uid: km44lv0j The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen responsible for coronavirus disease 2019 (COVID-19), has devastated public health services and economies worldwide. Despite global efforts to contain the COVID-19 pandemic, SARS-CoV-2 is now found in over 200 countries and has caused a death toll of over 1 million human lives as of November 2020. To date, only one Food and Drug Administration (FDA)-approved therapeutic drug (remdesivir) and a monoclonal antibody (MAb), bamlanivimab, are available for the treatment of SARS-CoV-2. As with other viruses, studying SARS-CoV-2 requires the use of secondary approaches to detect the presence of the virus in infected cells. To overcome this limitation, we have generated replication-competent recombinant SARS-CoV-2 (rSARS-CoV-2) constructs expressing fluorescent (Venus or mCherry) or bioluminescent (Nluc) reporter genes. Vero E6 cells infected with reporter-expressing rSARS-CoV-2 can be easily detected via fluorescence or luciferase expression and display a good correlation between reporter gene expression and viral replication. Moreover, rSARS-CoV-2 expressing reporter genes has plaque sizes and growth kinetics comparable to those of wild-type virus, rSARS-CoV-2/WT. We used these reporter-expressing rSARS-CoV-2 constructs to demonstrate their feasibility to identify neutralizing antibodies (NAbs) or antiviral drugs. Our results demonstrate that reporter-expressing rSARS-CoV-2 represents an excellent option to identify therapeutics for the treatment of SARS-CoV-2, where reporter gene expression can be used as a valid surrogate to track viral infection. Moreover, the ability to manipulate the viral genome opens the feasibility of generating viruses expressing foreign genes for their use as vaccines for the treatment of SARS-CoV-2 infection. IMPORTANCE Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen that causes coronavirus disease 2019 (COVID-19), has significantly impacted the human health and economic status worldwide. There is an urgent need to identify effective prophylactics and therapeutics for the treatment of SARS-CoV-2 infection and associated COVID-19. The use of fluorescent-protein- or luciferase-expressing reporter viruses has significantly advanced viral research. Here, we generated recombinant SARS-CoV-2 (rSARS-CoV-2) constructs expressing fluorescent (Venus and mCherry) or luciferase (Nluc) reporter genes and demonstrated that these viruses represent an excellent option to track viral infections in vitro. Importantly, reporter-expressing rSARS-CoV-2 constructs display growth kinetics and plaque phenotypes similar to those of their wild-type counterpart (rSARS-CoV-2/WT), demonstrating their usefulness for identifying drugs and/or neutralizing antibodies (NAbs) for the therapeutic treatment of SARS-CoV-2. Henceforth, these reporter-expressing rSARS-CoV-2 constructs can be used to interrogate large libraries of compounds and/or monoclonal antibodies (MAb), in high-throughput screening settings, to identify those with therapeutic potential against SARS-CoV-2. activities of NAbs and antivirals, respectively, obtained in our reporter-based microneutralization assays were similar to those observed in classical microneutralization assays using rSARS-CoV-2/WT (23) . These results demonstrate that our reporter-expressing rSARS-CoV-2 represent an excellent tool for studying the biology of the virus and for the identification of therapeutics for the treatment of SARS-CoV-2 infection and also for in vivo studies. Furthermore, because of reporter gene expression, these rSARS-CoV-2 constructs expressing reporter genes represent an ideal option to screen large libraries of biologicals to identify those with antiviral activity. Our results also demonstrate the feasibility of generating rSARS-CoV-2 expressing foreign genes that could be used to generate vaccines for the treatment of SARS-CoV-2 infections and/or associated COVID-19 disease. Generation of rSARS-CoV-2 expressing reporter genes. The pBeloBAC11 plasmid carrying the full-length viral genome of SARS-CoV-2 was previously described (21) . To generate the reporter-expressing rSARS-CoV-2, open reading frame (ORF) 7a was replaced with a Venus, mCherry, or Nluc gene in the pBeloBAC11 plasmid carrying the remaining viral genome to produce pBeloBAC11-SARS-CoV-2-del7a/Venus, -del7a/ mCherry, or -del7a/Nluc plasmids for viral rescue. We then used our previously described BAC-based reverse genetics approach to rescue rSARS-CoV-2-Venus, -mCherry, and -Nluc (Fig. 1A) . We confirmed the rescue of rSARS-CoV-2 expressing Venus, mCherry, or Nluc reporter genes by reverse-transcription-PCR (RT-PCR) using total RNA from mock-, rSARS-CoV-2/WT-, or rSARS-CoV-2 reporter virus-infected cells with primers specific for the viral NP, ORF 7a, or the individual reporter genes (Fig. 1B) . As expected, primers specific for SARS-CoV-2 NP amplified a band of ;1,260 bp from the RNA extracted from rSARS-CoV-2-infected but not mock-infected cells (Fig. 1B) . Bands that were amplified using primers in ORF 7a resulted in the expected ;566 bp in cells infected with rSARS-CoV- At 72 h postinfection, total RNA collected from cells was used to amplify, using RT-PCR, the viral NP, the ORF 7a region, and the different reporter genes (Venus, mCherry, or Nluc). Primers used for this RT-PCR analysis are shown on the left. The molecular sizes (base pairs) of the RT-PCR amplified products are shown on the right. 2/WT and ;920, 911, and 815 bp in the case of cells infected with rSARS-CoV-2-Venus, -mCherry and -Nluc, respectively, based on the different sizes of the reporter genes (Fig. 1B) . Primers specific for the reporter genes resulted in RT-PCR amplification of bands only from cells infected with the respective reporter-expressing rSARS-CoV-2 (Fig. 1B) . These results demonstrate that substitution of the viral ORF 7a with Venus, mCherry, or Nluc genes results in the successful recovery of rSARS-CoV-2 containing these reporter genes. Characterization of rSARS-CoV-2 expressing reporter genes. Next, we characterized the reporter-expressing rSARS-CoV-2 by evaluating the expression levels of Venus, mCherry, or Nluc in cell cultures and compared them to those of cells infected with rSARS-CoV-2/WT (Fig. 2) . rSARS-CoV-2 reporter viruses expressing Venus and mCherry were directly visualized under a fluorescence microscope ( Fig. 2A) . Indirect immunofluorescence microscopy using a MAb against SARS-CoV NP was used to detect rSARS- Vero E6 cells were mock infected or infected (MOI of 0.01) with WT and Venus-or mCherry-expressing rSARS-CoV-2. At 48 h postinfection, cells were fixed and permeabilized, visualized for Venus (left) or mCherry (right) expression, and immunostained with a SARS-CoV NP MAb (1C7). DAPI was used for nuclear staining. Merged images for Venus (left) or mCherry (right), viral NP, and DAPI are presented. Representative images are shown. Magnification, Â20. Bar, 100 mm. (B) Nluc expression. Vero E6 cells were mock infected or infected (MOI of 0.01) with WT and Nluc-expressing rSARS-CoV-2. At 48 h postinfection, Nluc expression in tissue culture supernatants was analyzed using a Synergy LX microplate reader (BioTek). (C) Western blotting. Vero E6 cells were mock infected or infected (MOI 0.01) with WT and Venus (left)-, mCherry (center)-, or Nluc (right)-expressing rSARS-CoV-2. At 48 h postinfection, viral NP and reporter gene protein expression levels were analyzed using specific antibodies. An antibody against beta-actin was used as an internal control. The sizes of molecular markers are shown on the right of each blot. CoV-2/WT infection ( Fig. 2A) . As expected, Venus and mCherry expression was observed only in Vero E6 cells infected with rSARS-CoV-2 expressing Venus or mCherry, respectively, but not in cells infected with rSARS-CoV-2/WT ( Fig. 2A) . Importantly, only cells infected with rSARS-CoV-2-Venus or rSARS-CoV-2-mCherry were detected using green or red filters, respectively (data not shown). As expected, the viral NP was detected in cells infected with rSARS-CoV-2-WT, -Venus, or -mCherry ( Fig. 2A) . Expression of Nluc in rSARS-CoV-2-Nluc-infected cells was evaluated from tissue culture supernatants at 48 h postinfection (Fig. 2B ). High levels of Nluc expression were detected in culture supernatants of cells infected with rSARS-CoV-2-Nluc but not from mock-or rSARS-CoV-2/WT-infected cells (Fig. 2B) . These results demonstrate that Vero E6 cells infected with rSARS-CoV-2-Venus, -mCherry, or -Nluc express the corresponding reporter genes and that viral infections can be detected by fluorescence (rSARS-CoV-2-Venus or -mCherry) or luciferase (rSARS-CoV-2-Nluc) without the need for antibodies, which were required for the detection of rSARS-CoV-2/WT. We next evaluated reporter protein expression levels by Western blotting assay using cell lysates from either mock-, rSARS-CoV-2-WT-, or rSARS-CoV-2-Venus-, -mCherry-, or -Nluc-infected cells using MAbs against the viral NP, the reporter genes, or actin as a loading control (Fig. 2C ). As expected, reporter gene expression was detected in lysates of cells infected with the respective reporter-expressing rSARS-CoV-2 but not from mock-or rSARS-CoV-2-WT-infected cells. Viral NP expression was detected in lysates from all virusinfected cells but not mock-infected cells (Fig. 2C ). Next, we assessed reporter gene expression over a period of 96 h in cells that were mock infected (data not shown) or cells infected with WT or reporter-expressing rSARS-CoV-2 (Fig. 3 ). Venus and mCherry expression levels were determined using fluorescence microscope (Fig. 3A) , while Nluc activity in tissue culture supernatants from infected cells was detected using a luminometer (Fig. 3B ). Venus and mCherry expression were detected as early as 24 h postinfection, and fluorescent protein expression increased over time until 96 h postinfection, when a decrease in fluorescence was observed because of cytopathic effect (CPE) caused by viral infection (bright-field [BF] microscopy) (Fig. 3A ). Similar CPE, but not fluorescence expression, was also observed in cells infected with rSARS-CoV-2/WT (Fig. 3A) . Levels of Nluc expression were also detected as early as 24 h postinfection and increased in a time-dependent matter (Fig. 3B ). To assess whether deletion of ORF 7a and insertion of reporter genes compromised viral fitness in cultured cells, we compared growth kinetics of reporter-expressing rSARS-CoV-2 to those of rSARS-CoV-2/WT (Fig. 3C ). We found that all the reporterexpressing rSARS-CoV-2 exhibited growth kinetics and peak viral titers of infection similar to those of rSARS-CoV-2/WT (Fig. 3C) , suggesting that deletion of ORF 7a and insertion of the reporter genes did not significantly affect viral fitness, at least in cultured cells. These results also support previous findings with SARS-CoV where deletion of ORF 7a and insertion of reporter genes did not impact viral fitness in vitro (24, 25) . These results were further confirmed when we evaluated the plaque phenotypes of the rSARS-CoV-2 constructs expressing fluorescent reporter genes and compared them to those of rSARS-CoV-2/WT (Fig. 3D ). Similar plaque sizes were observed in Vero E6 cells infected with rSARS-CoV-2/WT and rSARS-CoV-2 expressing Venus or mCherry (Fig. 3D) . Notably, Venus-positive or mCherry-positive plaques were detected only in cells infected with rSARS-CoV-2-Venus or -mCherry, respectively, and not in rSARS-CoV-2/WT-infected cells (Fig. 3D) . Importantly, fluorescent plaques overlapped those detected by immunostaining using the SARS-CoV NP MAb 1C7. Similar to the growth kinetics data, we found no significant differences between the plaque sizes of reporterexpressing rSARS-CoV-2 and rSARS-CoV-2/WT (Fig. 3D) . Reporter-based microneutralization assay for the identification of antivirals. To determine the feasibility of using our reporter-expressing rSARS-CoV-2 for the identification of antivirals, we evaluated the ability of remdesivir to inhibit SARS-CoV-2 in reporter-based microneutralization assays (Fig. 4 ). Remdesivir has been described to inhibit SARS-CoV-2 infection and is the only FDA-approved antiviral for the treatment of (23) . This demonstrates the feasibility of using these reporter-expressing rSARS-CoV-2 constructs and the reporterbased assay to easily identify compounds with antiviral activity based on fluorescentprotein or luciferase expression and without the need for MAbs to detect the presence of the virus in infected cells. Reporter-based microneutralization assay for the identification of NAbs. We next evaluated the feasibility of using our reporter-expressing rSARS-CoV in a reporterbased microneutralization assay to identify NAbs against SARS-CoV-2. As proof of concept, we used a human MAb (1212C2) which we recently described to potently bind and neutralize SARS-CoV-2 infection both in vitro and in vivo (27) . The half-maximal neutralizing concentrations (NT 50 s) of 1212C2 against rSARS-CoV-2-Venus (1.94 ng) ( Fig. 5A and B) , -mCherry (5.02 ng) ( Fig. 5C and D) , and -Nluc (3.67 ng) (Fig. 5E ) were similar to those observed with rSARS-CoV-2/WT (4.88 ng) (Fig. 5F ) and to recently reported values (27) . Genetic stability of rSARS-CoV-2 in vitro. The genetic stability of reporter-expressing recombinant viruses is important to demonstrate their viability in in vitro and/or in vivo studies. To evaluate the ability of our rSARS-CoV-2 constructs to maintain fluorescent reporter gene expression, viruses were consecutively passaged in Vero E6 cells, and Venus and mCherry expression was determined by plaque assay using fluorescence microscopy (Fig. 6A) . To that end, we evaluated fluorescence expression of over 40 plaques before immunostaining with anti-SARS-CoV NP MAb 1C7. We found that Venus and mCherry fluorescence expression from our rSARS-CoV-2 reporter viruses was genetically stable with nearly 100% of the plaques analyzed under a fluorescence microscope (Fig. 6A) . We also evaluated the complete genome sequences of the reporter-expressing rSARS-CoV-2 constructs used in our studies (passage 3 [P3]) with those of additional passages (P4 and P5) using next-generation sequencing (Fig. 6B ). In the case of rSARS-CoV-2-Venus (Fig. 6B, top) , few variants were found at low frequencies after two additional passages (P5), indicating no significant changes and/or deletions in the viral genome. However, for rSARS-CoV-2-mCherry (Fig. 6B, middle) , variants containing a mutation at position 21784 in the S gene were found in our viral stock (P3), and the frequency of this mutation increased after additional passages (P4 and P5). Two additional mutations at positions 23525 and 24134 (both in the S gene) were also found at P5. In the case of rSARS-CoV-2-Nluc (Fig. 6B, bottom) , a mutation at position 24755 (S gene) was found in our viral stock (P3). The frequency of this mutation increased up to 100% after 2 additional passages (P5). Other, less abundant mutations at positions 13419 (nsp12, RNA-dependent RNA polymerase), 23525 (S gene), and 26256 (envelope [E] gene) were also found after 2 additional passages (P5) (Fig. 6B, bottom) . Most likely, these mutations arose due to viral adaptation to Vero E6 cells. However, it is also possible that these mutations could be related to the nature of the Journal of Virology reporter gene, since they were different between the three reporter-expressing rSARS-CoV-2 constructs. Reporter-expressing viruses represent a powerful tool for both basic research and translational studies (16, 17, (28) (29) (30) (31) . Several research groups, including ours, have described recombinant viruses expressing reporter genes to easily study the biology of viral infections, to evaluate the efficacy of antivirals or NAbs, and for in vivo studies in validated animal models (16, (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) . Both fluorescent and luciferase proteins have been used to generate reporterexpressing viruses. However, the innate and differing properties of reporter genes dictate which one might be inserted into a recombinant virus. While fluorescent proteins provide an efficient way to track viral infections using microscopy, luciferase proteins are more readily quantifiable and therefore more amenable to HTS studies (16, 17, 46) . For this reason, in this study we generated rSARS-CoV-2 expressing fluorescent proteins (Venus and mCherry) or luciferase (Nluc) (Fig. 1) . These reporter genes were selected either because of their distinctive fluorescent properties (Venus and mCherry) or because of their small size, stability, high bioluminescence activity, and ATP independence (Nluc) (47) . Although reporter-expressing rSARS-CoV-2 viruses similar to those reported here were recently described (18) (19) (20) , this is the first report of a replication-competent rSARS-CoV-2 expressing mCherry. Recombinant viruses expressing a red fluorescent protein represent an advantage over those expressing GFP or mNeonGreen (18) (19) (20) in that many genetically modified cell lines and/or animals express green fluorescent proteins. Another limitation of green fluorescent proteins during in vivo imaging is the absorption of the fluorophores' excitation and emission by hemoglobin and autofluorescence of tissues (48) (49) (50) (51) . Recombinant viruses expressing red fluorescent proteins represent a better option to combine with genetically modified GFP-expressing cell lines and/or animals and, based on their reduced autofluorescence background, to more accurately capture the dynamics of viral infection and replication. Reporter-expressing replicating competent viruses can be used to monitor viral infections, assess viral fitness, and evaluate and/or identify antivirals and/or NAbs, where reporter gene expression can be used as a valid surrogate for virus detection in infected cells. Expression of Venus, mCherry, or Nluc from our rSARS-CoV-2 was confirmed by directly visualizing fluorescence expression under a fluorescence microscope (Venus and mCherry) or luciferase activity (Nluc) using a microplate reader ( Fig. 2 and 3 ). Western blot analyses using specific antibodies against each of the reporter genes further confirmed expression from their respective rSARS-CoV-2 constructs ( Fig. 2 and 3) . Notably, despite deletion of ORF 7a and insertion of a reporter gene, the three reporter-expressing rSARS-CoV-2 constructs displayed growth kinetics and plaque phenotype similar to those of their WT counterpart (Fig. 3) . As expected, viral infection was visualized in real time, without the need for secondary approaches (e.g., MAbs) to detect the presence of the virus in infected cells. Overall, reporter gene expression displayed kinetics that correlated with levels of viral replication, further demonstrating the feasibility of using these reporter genes as a valid surrogate to assess viral infection. Therapeutic treatment of SARS-CoV-2 infection has been limited to the use of remdesivir (3) and the MAb bamlanivimab and, more recently, to two FDA-approved vaccines (4, 5) . Notably, there is a possibility, similar to the situation with other respiratory viruses (e.g., influenza virus), of the emergence of drug-resistant SARS-CoV-2 variants that will impose a significant additional challenge in the ongoing COVID-19 FIG 6 Legend (Continued) from ;40 to 50 viral plaques per passage. Representative images of immunostaining and fluorescence obtained from each P3 to P5 viral plaque are shown. Magnification, Â4. Bar, 750 mm. (B) Sequence analysis. Reporter-expressing rSARS-CoV-2 nonreference allele frequencies from virus stock (P3) and after two consecutive passages in Vero cells (P4 and P5) were determined using next-generation sequencing, using modified rSARS-CoV-2/WT reference genomes. Nonreference alleles that were below 1% of reads are not shown, and those greater than 20% are indicated in red. pandemic (52) . Thus, it is imperative to discover not only new antivirals and other therapeutic approaches but also prophylactics for the treatment of SARS-CoV-2 infections. To that end, rapid and sensitive screening assays to identify compounds with antiviral activity or to assess efficacy of vaccine candidates for the therapeutic and prophylactic treatment of SARS-CoV-2 infections, respectively, are urgently needed. In this study, we demonstrate that reporter-expressing rSARS-CoV-2 constructs represent an excellent option for the rapid identification and characterization of both antivirals (Fig. 4) and NAbs (Fig. 5) for the therapeutic and/or prophylactic treatment of SARS-CoV-2 infections. Importantly, EC 50 s (antivirals) and NT 50 s (NAbs) obtained with our reporterexpressing viruses were comparable to those obtained using rSARS-CoV-2/WT or described by others (20, 23, 27) , demonstrating the feasibility of using our reporterbased microneutralization assays for the rapid identification of antivirals or NAbs ( Fig. 4 and 5, respectively) . Furthermore, our results indicate that Venus-, mCherry-, and Nluc-expressing rSARS-CoV-2 constructs are stable up to 5 passages in vitro in Vero E6 cells, including expression of the reporter gene (Fig. 6) . We have not yet conducted studies to evaluate the feasibility of using these reporter-expressing rSARS-CoV-2 constructs in vivo. It is possible that, as with other respiratory viruses, rSARS-CoV-2 expressing reporter genes could also be used to study the biology of viral infections in validated animal models of viral infection. Our SARS-CoV-2 reverse genetics based on the use of a BAC have allowed us to rescue rSARS-CoV-2/WT (21) and rSARS-CoV-2 stably expressing reporter genes. In the case of our reporter-expressing rSARS-CoV-2 constructs, we removed ORF 7a and replaced it with various reporter genes without a significant impact on viral replication. The feasibility of removing viral genes and inserting reporter genes demonstrates the genetic plasticity of the SARS-CoV-2 genome and opens the possibility of generating recombinant viruses expressing other genes of interest for the development of SARS-CoV-2 vaccines that could be used for the control of the ongoing COVID-19 pandemic. Biosafety. All experiments involving infectious SARS-CoV-2 were performed in a biosafety level 3 (BSL3) laboratory at the Texas Biomedical Research Institute. Protocols involving SARS-CoV-2 were approved by the Texas Biomedical Research Institute's Institutional Biosafety Committee (IBC). Cell lines. African green monkey kidney epithelial cells (Vero E6; CRL-1586) were grown and maintained in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 1% PSG (100 U/ml penicillin, 100 mg/ml streptomycin, and 2 mM L-glutamine) at 37°C with 5% CO 2 . Generation of pBeloBAC11-SARS-CoV-2 carrying reporter genes. The plasmid pBeloBAC11 (New England Biolabs) containing the entire viral genome of SARS-CoV-2 was described previously (21, 53) . Briefly, the entire genome sequence of SARS-CoV-2 USA/WA1/2020 (GenBank accession no. MN985325) was chemically synthesized (Bio Basic) in five fragments and cloned into pUC57 plasmids containing unique restriction sites. Silent mutations were introduced to the spike (S) and matrix (M) genes to remove BstBI and Mlul restriction sites, respectively, that were used for the assembly of the entire SARS-CoV-2 genome into the pBeloBAC11 plasmid. These nucleotide changes were also used as genetic markers to distinguish the natural USA/WA1/2020 and the recombinant SARS-CoV-2 (21) . The five fragments containing the entire SARS-CoV-2 genome were assembled into pBeloBAC11 using standard molecular biology techniques. To remove the 7a gene and introduce the Venus, mCherry, or Nluc reporter genes, the region flanking the 7a viral gene and each individual reporter gene were amplified by extension and overlapping PCR using specific oligonucleotides in a shuttle plasmid. The modified 7a viral genes were inserted into the pBeloBAC11 plasmid containing the remaining SARS-CoV-2 viral genome using BamHI and RsrII restriction sites to generate pBeloBAC11-SARS-CoV-2-del7a/Venus, pBeloBAC11-SARS-CoV-2-del7a/mCherry, and pBeloBAC11-SARS-CoV-2-del7a/Nluc for the rescue of rSARS-CoV-2-Venus, rSARS-CoV-2-mCherry, and rSARS-CoV-2-Nluc, respectively. Plasmids and pBeloBAC11 constructs were validated by Sanger sequencing (ACGT Inc.). Rescue of rSARS-CoV-2 expressing reporter genes. The rSARS-CoV-2/WT and rSARS-CoV-2 expressing reporter genes were rescued as previously described (21, 22) . Briefly, confluent monolayers of Vero E6 cells (1.2 Â 10 6 cells/well, 6-well plate format, triplicates) were transfected using Lipofectamine 2000 (LPF2000; Thermo Fisher) with 4 mg/well of pBeloBAC11-SARS-CoV-2/WT or pBeloBAC11-SARS-CoV-2-del7a/Venus, -del7a/mCherry, or -del7a/Nluc. An empty pBeloBAC11 plasmid was included as an internal control. At 14 h, transfection medium was replaced with postinfection medium (DMEM with 2% FBS), and 24 h later, cells were scaled up into T75 flasks. At 72 h, P0 virus-containing tissue culture supernatants were collected and stored at 280°C. Viral rescues were confirmed by infecting fresh Vero E6 cells (1.2 Â 10 6 cells/well, 6-well plates, triplicates) and assessing fluorescence or Nluc expression. P0 viruses were passaged three times, and viral stocks were generated and titrated for in vitro experiments. Viral titers (in PFU per milliliter) were determined by plaque assay in Vero E6 cells (1.2 Â 10 6 cells/well, 6-well plate format). Sequencing. Viral RNAs from Vero E6 cells (1.2 Â 10 6 cells/well, 6-well plate format) infected at a multiplicity of infection (MOI) of 0.01 were extracted using TRIzol reagent (Thermo Fisher Scientific), according to the manufacturer's specifications. Libraries were generated with a KAPA RNA HyperPrep kit, 100 ng of RNA, and 7 mM adapter. The Illumina HiSeq X was used for sequencing. Raw reads were filtered using Trimmomatic v0.39 (54) . SARS-CoV-2 templates were made for each reporter gene by modifying SARS-CoV-2 USA/WA1/2020 (GenBank accession no. MN985325.1). Modifications included deleting ORF 7a, adding T21895C and T26843A mutations, and inserting the appropriate reporter gene (Venus, mCherry, or Nluc) at position 27937. Reads were mapped to the modified SARS-CoV-2 templates with Bowtie v2.4.1 (55) , and the total genomic coverage was quantified using MosDepth v0.2.6 (56). Allele frequencies were estimated with LoFreq* v2.1.3.1 (57) , and low-frequency variants with less than a 100Â read depth or a 1% minor allele frequency were eliminated. RT-PCR. Total RNA from mock-or virus-infected (MOI of 0.01) Vero E6 cells (1.2 Â 10 6 cells/well, 6well plate format) were extracted using TRIzol reagent (Thermo Fisher Scientific). Superscript II reverse transcriptase (Invitrogen) and Expand high-fidelity PCR (Sigma-Aldrich) were used to synthesize and amplify the cDNAs, respectively, using primers specific for the viral nucleoprotein (NP) or ORF 7a region and Venus, mCherry, or Nluc. Immunofluorescence assays (IFA). Confluent monolayers of Vero E6 cells (1.2 Â 10 6 cells/well, 6well format, triplicates) were mock infected or infected (MOI of 0.01) with rSARS-CoV-2 expressing Venus or mCherry or rSARS-CoV-2/WT. At 48 h postinfection, cells were fixed with 10% neutral buffered formalin at 4°C for 16 h for fixation and viral inactivation and permeabilized with phosphate-buffered saline (PBS) containing 0.5% (vol/vol) Triton X-100 for 5 min at room temperature. Cells were washed with PBS and blocked with 2.5% bovine albumin serum (BSA) in PBS for 1 h before incubation with 1 mg/ml of SARS-CoV anti-NP MAb 1C7 in 1% BSA in PBS for 1 h at 37°C. Cells infected with rSARS-CoV-2-Venus or -mCherry were washed with PBS and stained with either Alexa Fluor 594 goat anti-mouse IgG (Invitrogen; 1:1,000) or fluorescein isothiocyanate (FITC)-conjugated goat anti-mouse IgG (Dako; 1:200), respectively. Cell nuclei were stained with 49,6-diamidino-2-phenylindole (DAPI; Research Organics). Representative images were captured using a fluorescence microscope (EVOS M5000 imaging system) at a magnification of Â20. Protein gel electrophoresis and Western blots. Vero E6 cells (1.2 Â 10 6 cells/well, 6-well plate format, triplicates) were mock-infected or infected (MOI of 0.01) with rSARS-CoV-2/WT or rSARS-CoV-2 expressing Venus, mCherry, or Nluc. At 48 h postinfection, cells were lysed with 1Â passive lysis buffer (Promega), and proteins were separated by denaturing electrophoresis in 12% SDS-polyacrylamide gels and transferred to a nitrocellulose membrane (Bio-Rad) with a Mini-Protean Tetra vertical electrophoresis cell at 100 V for 1 h at 4°C. Membranes were blocked in PBS containing 10% dried skim milk and 0.1% Tween 20 for 1 h and then incubated overnight at 4°C with the following primary MAbs or polyclonal antibodies (PAbs): SARS-CoV NP (mouse MAb 1C7; Thomas Moran, Icahn School of Medicine at Mount Sinai), Venus (rabbit PAb sc-8334; Santa Cruz Biotechnology), mCherry (rabbit PAb; RayBiotech), and Nluc (rabbit PAb; Promega). A MAb against actin (MAb AC-15; Sigma) was included as a loading control. Primary antibodies bound to the membrane were detected using horseradish peroxidase (HRP)-conjugated secondary antibodies against mouse or rabbit (GE Healthcare). Proteins were detected by chemiluminescence using SuperSignal West Femto Maximum Sensitivity substrate (Thermo Scientific) based on the manufacturer's specifications and imaged in a ChemiDoc imaging system (Bio-Rad). Plaque assays and immunostaining. Confluent monolayers of Vero E6 cells (1.2 Â 10 6 cells/well, 6well plate format, triplicates) were infected with WT or reporter-expressing rSARS-CoV-2 for 1 h at 37°C. After viral absorption, infected cells were overlaid with agar and incubated at 37°C for 72 h. Afterwards, cells were submerged in 10% neutral buffered formalin at 4°C for 16 h for fixation and viral inactivation, and then the agar overlays were gently removed. To observe Venus and mCherry fluorescence expression, PBS was added to each well and plates were imaged under a fluorescence microscope (EVOS M5000 imaging system). For immunostaining, plates were permeabilized with 0.5% Triton X-100-PBS for 10 min at room temperature, blocked with 2.5% BSA-PBS for 1 h at room temperature, and then incubated at 37°C for 1 h with the anti-SARS 2 NP MAb 1C7. Plaques were developed for visualization using the Vectastain ABC kit and DAB HRP substrate kit (Vector Laboratories), in accordance with the manufacturer's recommendations. Viral growth kinetics and titrations. Vero E6 cells (1.2 Â 10 6 cells/well, 6-well plate format, triplicates) were infected (MOI of 0.01) with rSARS-CoV-2/WT or rSARS-CoV-2 expressing Venus, mCherry, or Nluc. After viral adsorption for 1 h at 37°C, cells were washed with PBS, provided with fresh postinfection medium, and then placed in a 37°C incubator with a 5% CO 2 atmosphere. At various times postinfection (12, 24, 48, 72 , and 96 h), cells were imaged for Venus or mCherry expression under a fluorescence microscope (EVOS M5000 imaging system). Viral titers in the tissue culture supernatants at each time point were determined by titration and immunostaining, as previously described, using the anti-SARS-CoV NP MAb 1C7. Nluc expression in tissue culture supernatants was quantified using Nano-Glo luciferase substrate (Promega) following the manufacturer's recommendations. Mean values and standard deviations (SD) were determined using GraphPad Prism software (version 8.2). Reporter-based microneutralization assay for the identification of antivirals. Vero E6 cells (96well plate format, 4 Â 10 4 cells/well, quadruplicates) were infected with ;100 to 200 PFU of rSARS-CoV-2/WT or rSARS-CoV-2 expressing Venus, mCherry, or Nluc for 1 h at 37°C. After viral adsorption, cells were washed and incubated in 100 ml of infection medium (DMEM with 2% FBS) containing 3-fold serial dilutions (starting concentration of 50 mM) of remdesivir or 0.1% dimethyl sulfoxide (DMSO) as a vehicle Cells infected with rSARS-CoV-2/WT or rSARS-CoV-2 expressing fluorescent Venus or mCherry were incubated at 37°C for 24 h, while cells infected with rSARS-CoV-2 expressing Nluc were incubated at 37°C for 48 h. For rSARS-CoV-2/WT and rSARS-CoV-2 expressing fluorescent Venus and mCherry, cells were submerged in 10% neutral buffered formalin at 4°C for 16 h for fixation and viral inactivation. Cells were washed with 100 ml/well of PBS three times, permeabilized with 100 ml/well of 0.5% Triton X-100 in PBS at room temperature for 15 min PBS at 37°C for 1 h. Next, cells were staining with the anti-NP MAb 1C7 (1 mg/ml Viral infections were determined using fluorescent images of each well and quantified using Cell Profiler (Broad Institute) cell image analysis software. In the case of cells infected with rSARS-CoV-2 expressing Nluc, tissue culture supernatants were collected at 48 h postinfection, and Nluc expression was measured confluent monolayers of Vero E6 cells (96-plate format, 4 Â 10 4 cells/well, quadruplicates) were infected with ;100 to 200 PFU of rSARS-CoV-2/WT or rSARS-CoV-2 expressing Venus, mCherry, or Nluc for 1 h at 37°C. After viral adsorption, cells were washed and incubated with 100 ml of infection medium (DMEM with 2% FBS) containing 3-fold serial dilutions (starting concentration of 500 ng) of 1212C2 or PBS and 1% Avicel (Sigma-Aldrich) Cells were next incubated with the anti-NP MAb 1C7 (1 mg/ml) in 1% BSA-PBS at 37°C for 1 h. Cells were next washed three times with PBS and incubated with a secondary FITC-conjugated goat anti-mouse IgG (Dako; 1:200) in 1% BSA for 1 h at 37°C. Cell nuclei were stained with DAPI (Research Organics). Viral infections were determined using fluorescent images of each well and quantified using Cell Profiler (Broad Institute). In the case of cells infected with rSARS-CoV-2 expressing Nluc, tissue culture supernatants were collected at 48 h postinfection, and Nluc expression was measured using a lucif Individual wells from three independent experiments conducted in quadruplicate were used to calculate the average and standard deviation (SD) of viral inhibition using Microsoft Excel software. Nonlinear regression curves, the half-maximal effective concentration (EC 50 ) of remdesivir, and the half-maximal neutralizing concentration (NT 50 ) of 1212C2 were determined using GraphPad Prism software China Novel Coronavirus Investigating and Research Team. 2020. A novel coronavirus from patients with pneumonia in China An interactive web-based dashboard to track COVID-19 in real time C4591001 Clinical Trial Group. 2020. Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine The Advisory Committee on Immunization Practices' interim recommendation for use of Moderna COVID-19 vaccine -United States Epidemiology, genetic recombination, and pathogenesis of coronaviruses Coronaviruses and SARS-CoV-2: a brief overview SARS and MERS: recent insights into emerging coronaviruses A novel coronavirus associated with severe acute respiratory syndrome Origin and evolution of pathogenic coronaviruses The pandemic in pictures: how coronavirus is changing the world Members of Steering Committee, Society of Global Health, Chinese Preventive Medicine Association. 2020. Fighting against the common enemy of COVID-19: a practice of building a community with a shared future for mankind COVID-19: a promising cure for the global panic The COVID-19 pandemic: important considerations for contact lens practitioners COVID-19 outbreak: an overview Replication-competent fluorescent-expressing influenza B virus A novel fluorescent and bioluminescent bireporter influenza A virus to evaluate viral infections SARS-CoV-2 reverse genetics reveals a variable infection gradient in the respiratory tract An infectious cDNA clone of SARS-CoV-2 A nanoluciferase SARS-CoV-2 for rapid neutralization testing and screening of anti-infective drugs for COVID-19 Rescue of SARS-CoV-2 from a single bacterial artificial chromosome Rescue of recombinant Zika virus from a bacterial artificial chromosome cDNA clone Rapid in vitro assays for screening neutralizing antibodies and antivirals against SARS-CoV-2 SARS coronavirus accessory proteins 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 Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro Therapeutic activity of an inhaled potent SARS-CoV-2 neutralizing human monoclonal antibody in hamsters Development of a rapid antiviral screening assay based on eGFP reporter virus of Mayaro virus Generation and validation of a highly sensitive bioluminescent HIV-1 reporter vector that simplifies measurement of virus release Development and characterization of a stable eGFP enterovirus 71 for antiviral screening Development and characterization of a stable luciferase dengue virus for high-throughput screening Replication-competent influenza A virus that encodes a split-green fluorescent protein-tagged PB2 polymerase subunit allows live-cell imaging of the virus life cycle Replication-competent influenza A viruses expressing reporter genes Replication-competent influenza A and B viruses expressing a fluorescent dynamic timer protein for in vitro and in vivo studies In vivo imaging of influenza virus infection in immunized mice Influenza A virus encoding secreted Gaussia luciferase as useful tool to analyze viral replication and its inhibition by antiviral compounds and cellular proteins Multi-spectral fluorescent reporter influenza viruses (Color-flu) as powerful tools for in vivo studies In vivo bioluminescent imaging of influenza a virus infection and characterization of novel cross-protective monoclonal antibodies Visualizing real-time influenza virus infection, transmission and protection in ferrets Analysis of in vivo dynamics of influenza virus infection in mice using a GFP reporter virus Hemagglutinin-pseudotyped green fluorescent protein-expressing influenza viruses for the detection of influenza virus neutralizing antibodies Replication-competent influenza A viruses expressing a red fluorescent protein Insertion of a GFP reporter gene in influenza virus Generation of a variety of stable influenza A reporter viruses by genetic engineering of the NS gene segment Multi-modal imaging with a toolbox of influenza A reporter viruses Cellular bioluminescence imaging Engineered luciferase reporter from a deep sea shrimp utilizing a novel imidazopyrazinone substrate Comparative study reveals better far-red fluorescent protein for whole body imaging Bright far-red fluorescent protein for whole-body imaging Whole-body optical imaging of green fluorescent protein-expressing tumors and metastases Improved in vivo whole-animal detection limits of green fluorescent protein-expressing tumor lines by spectral fluorescence imaging Remdesivir targets a structurally analogous region of the Ebola virus and SARS-CoV-2 polymerases Generation of recombinant SARS-CoV-2 using a bacterial artificial chromosome Trimmomatic: a flexible trimmer for Illumina sequence data Fast gapped-read alignment with Bowtie 2 Mosdepth: quick coverage calculation for genomes and exomes LoFreq: a sequence-quality aware, ultrasensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets We thank Thomas Moran at the Icahn School of Medicine at Mount Sinai for providing us with the SARS-CoV cross-reactive NP MAb 1C7. We also thank BEI Resources for providing the SARS-CoV-2 USA-WA1/2020 isolate (NR-52281) and Marina McDew-White and Robbie Diaz for constructing the NGS libraries. Finally, we also thank members at our institutes for their efforts in keeping them fully operational during the COVID-19 pandemic and the BSC and IACUC for reviewing our protocols in a time-efficient manner.