Int J Biol Sci 2021; 17(6):1521-1529. doi:10.7150/ijbs.59151
Developing effective siRNAs to reduce the expression of key viral genes of COVID-19
Faculty of Health Sciences, University of Macau, Macao SAR, China
Wu R, Luo KQ. Developing effective siRNAs to reduce the expression of key viral genes of COVID-19. Int J Biol Sci 2021; 17(6):1521-1529. doi:10.7150/ijbs.59151. Available from https://www.ijbs.com/v17p1521.htm
The COVID-19 pandemic has been raging worldwide for more than a year. Many efforts have been made to create vaccines and develop new antiviral drugs to cope with the disease. Here, we propose the application of short interfering RNAs (siRNAs) to degrade the viral genome, thus reducing viral infection. By introducing the concept of the probability of binding efficiency (PBE) and combining the secondary structures of RNA molecules, we designed 11 siRNAs that target the consensus regions of three key viral genes: the spike (S), nucleocapsid (N) and membrane (M) genes of SARS-CoV-2. The silencing efficiencies of the siRNAs were determined in human lung and endothelial cells overexpressing these viral genes. The results suggested that most of the siRNAs could significantly reduce the expression of the viral genes with inhibition rates above 50% in 24 hours. This work not only provides a strategy for designing potentially effective siRNAs against target genes but also validates several potent siRNAs that can be used in the clinical development of preventative medication for COVID-19 in the future.
Keywords: COVID-19, SARS-CoV-2, RNA secondary structure, siRNA, gene silencing
COVID-19 has infected over 122 million people and caused over 2.7 million human deaths worldwide till March 2021 . COVID-19 is caused by a positive-sense single-stranded RNA virus named SARS-CoV-2 which has a genome size of approximately 30 kilobases and belongs to the family Coronaviridae . SARS-CoV-2 uses a spike protein localized on the coat of the viral particle to bind with angiotensin converting enzyme 2 (ACE2) on the membrane of human epithelial cells in the lungs, arteries, heart, kidney, and intestines . Upon entry into the cells, the virus will produce RNA-dependent RNA polymerase (RdRP) to synthesize the antisense or negative strand of its genomic RNA for viral replication. In addition, viral RNA can be transcribed into subgenomic RNAs, which will be translated into structural proteins, including spike (S), nucleocapsid (N), and membrane (M) proteins to assemble into more viral particles .
One unique feature of RNA viruses is that they use RNA polymerase to synthesize their genetic material, and this enzyme is suggested to have a high error rate during RNA replication . This means that when viral RNAs replicate in host cells, new mutated forms of the virus can be generated and subsequently spread to other people. It is estimated that two new mutations can accumulate in the viral genome of SARS-CoV-2 on a monthly basis . Some of these mutations may make the new viral strains more contagious and less responsive to the vaccines that were developed based on the original version of the SARS-CoV-2 virus. For example, several more virulent variants of SARS-CoV-2 have been found in the UK, Brazil, and South Africa, and are now spreading to many other countries .
As the pace of developing new antiviral drugs seems behind the progression of the COVID-19 pandemic, many people have pinned their hope on vaccination [8,9]. However, according to reports, a considerable population is unwilling to be vaccinated, both in developed and developing countries, which will slow the establishment of community immunity . In addition, certain susceptible populations have been suggested as inapplicable recipients of the vaccine, including women who are pregnant or within the lactation period and patients with severe basic illnesses or malignant tumors. Thus, developing prompt and flexible antiviral therapies will contribute to the elimination of the COVID-19 pandemic.
An effective way to eliminate RNA-based viruses is to degrade their genetic material by using RNA interference (RNAi), especially small interfering RNAs (siRNAs) or short hairpin RNAs (shRNAs) [11-13]. RNAi-mediated gene silencing occurs when double-stranded siRNA molecules enter cells and are recognized by an RNA-induced silencing complex (RISC) that uses one strand of the siRNA to bind to the complementary strand of the target RNA via base pairing. Afterwards, the RISC uses its RNase activity to cleave the target RNA . RNAi technology has been widely used to silence gene expression in human cells in the last twenty years [15-18]. The past several years have witnessed the rapid growth of RNAi-based clinical trials for viral infections such as Epstein-Barr virus (EBV) , hepatitis B virus (HBV) , human immunodeficiency virus (HIV) , SARS-CoV  and the raging SARS-CoV-2 .
One of the challenges of using RNAi to reduce gene expression is that the length of an siRNA is very short, containing only 19-23 nucleotides, while the length of the target RNA strand is much longer, ranging from several hundred to thousands of nucleotides long. Thus, it is difficult to choose target sites on the RNA to design siRNAs. Due to the complexity of intracellular RNAi processing, as well as the dynamic conformation of target RNA molecules in different physiological statuses, the existing siRNA design tools can only partially predict plausible targets and still requires extensive validation to achieve better efficacy in a wet lab . In this study, we developed a method to design effective siRNAs against three viral genes of SARS-CoV-2 and tested their efficacies in reducing viral gene expression in human cells.
Materials and methods
Human umbilical vascular endothelial cells (HUVECs) and human lung carcinoma A549 cells were purchased from the American Type Culture Collection (ATCC, Rockville, MD, USA). All cells were cultured in DMEM (#12100046, Thermo Fisher Scientific, USA) supplemented with 10% fetal bovine serum (#10270106, Gibco, USA) and 1% penicillin-streptomycin antibiotics (#15140122, Thermo Fisher Scientific, USA). All cells tested negative for mycoplasma contamination.
Transfection of plasmid and siRNA
An overexpression plasmid for a mutant variant of the S (D614G) gene (#GS-200519_A012, pCDNA3.1(+)-S-D614G) was purchased from Genewiz (Suzhou, China); plasmids of the N (#HH-gene-158, pCDNA3.1(+)-2019-nCoV-N) and M (#HH-gene-135, pCDNA3.1(+)-2019-nCoV-M) genes were purchased from HedgehogBio (Shanghai, China). siRNAs were purchased from GeneralBio (Chuzhou, China). For transfection of plasmids and siRNAs, 1×105 cells were seeded into each well of a 12-well plate and grown overnight. Afterwards, 30 ng of plasmid DNA and 10 pmol of siRNA (at a final concentration of 10 nM) were added to 800 µL Opti-MEM™ (#22600050, Gibco, USA). Then, 1.6 µL Lipofectamine™ 2000 (#11668019, Invitrogen, Waltham, MA, USA) was added to the plasmid:siRNA solution and the mixture was incubated at room temperature for 30 min. Then, the mixture was added to the cells and incubated for 6 hours at 37ºC. After 18 hours of incubation, cells were harvested for RNA isolation and qPCR measurement.
Real-time quantitative polymerase chain reaction (qPCR)
RNA was extracted with TRIzol (#15596026, Thermo Fisher Scientific, USA) and isopropanol precipitation methods. cDNA was synthesized using an iScript™ cDNA Synthesis Kit (#1778890, Bio-Rad, USA), and qPCR was performed using iTaq™ Universal SYBR® Green (#1725122, Bio-Rad, USA) and analyzed with an Applied Biosystems™ 7500 Fast Real-Time PCR Detection System (Thermo Fisher Scientific, USA). The changes in mRNA levels were calculated using the ΔΔCt method with the mRNA level of GAPDH as an internal control. The used primers are listed in Table 1.
|Target gene||Direction||Primer (5' - 3')|
|Spike (S) gene||Forward||GCTGGTGCTGCAGCTTATTA|
|Nucleocapsid (N) gene||Forward||GGGGAACTTCTCCTGCTAGAAT|
|Membrane (M) gene||Forward||TGTGACATCAAGGACCTGCC|
Experiments were repeated at least four times, and the data were analyzed using one-way ANOVA as processed by GraphPad Prism 8 software (San Diego, CA, USA). All the analyzed results are presented as the mean ± SD (standard deviation). Statistical significance is shown as *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001; ns, not significant.
Selection of nine variants from the major prevalent viral strains worldwide
Since the beginning of the COVID-19 pandemic in January 2020, many variants of the SARS-CoV-2 virus have been detected in over 200 countries . Based on the times and locations of the viral strains that appeared in the past 12 months, we prepared a phylogenetic tree and a pangolin lineage of COVID-19 from the GISAID platform using Nextstrain open source project [25,26], which showed that this virus has quickly mutated and spread across the five continents (Figure 1A and B). According to the frequency of viral mutations and the trend of prevalence between late 2020 and early 2021, nine strains were selected, which account for 74% of the total prevalent strains in the world (Supplementary Table 1). Six strains, B.1 (5%), B.1.1 (5%), B.1.177 (5%), B.1.1.50 (3%), B.1.2 (7%), and B.1.5 (2%), are prevalent in the Americas, giving a total percentage of 27%. For the other three strains, B.1.351 (23%) is prevalent in Africa, B.1.1.7 (17%) is prevalent in Europe, especially in the UK, and B.1.36 (7%) is prevalent in Asia (Figure 1C).
Prevalence of major viral strains of COVID-19 in 2020 and early 2021. A. The phylogenetic tree of COVID-19 from 2019-Dec to 2021-Jan plotted by the GISAID platform. Each viral strain was colored based on the pangolin lineage. B. The geographic distribution of major viral strains between 2020-Oct and 2021-Jan. C. The evolution of COVID-19 from 2020-Jan to 2021-Jan based on the frequency of major viral strains that were sequenced chronologically.(Click on the image to enlarge.)
Calculation of the probability of binding efficiency (PBE) based on single-stranded information of RNA molecules
Although RNA is synthesized into a single strand, its complementary nucleotides, such as A=U and G≡C, can bind to each other via hydrogen bonds to produce different secondary structures, such as loops, stems, and hairpins. Previously, we reported that nucleotides located in the loop structure are in a single-stranded configuration and are easier for siRNAs to target, while nucleotides located in the stem region are in a double-stranded configuration and are harder for siRNA to target . As the binding between siRNA and the target RNA is crucial for RISC to degrade target RNA, these single-strand loop regions can be selected as potential sites for designing siRNAs.
To identify effective target sites on the RNA molecule, we designed a formula to calculate the probability of binding efficiency (PBE) between an siRNA and its target RNA based on the single-strand information obtained from the Mfold web server . PBE calculates the total number of hydrogen bonds between 21 nucleotides of an siRNA and the complementary region of the target RNA based on all the predicted secondary structures (Figure 2A).
We then calculated the PBE values of all possible siRNA candidates against the spike (S), nucleocapsid (N) and membrane (M) viral genes. The distribution curves show that the PBE values of most of the siRNA candidates are below 30, while the PBE values of a small portion of the siRNA candidates are above 30 (Figure 2B). It should be emphasized that for different genes with different length and the number of secondary structures, the values of PBE can differ significantly. Thus, there is no universal cut-off to indicate the efficacy of siRNAs. We recommend choosing the ones with the highest PBE values as candidates. We then selected the 3-4 top-ranked siRNAs with the highest PBE values from individual viral genes for structural and functional assessments.
siRNAs with high PBE values can effectively reduce the expression of the spike gene
We used the RNA sequence of the S gene from the B.1 strain, as this strain is currently the origin of the most prevalent strains in the world . The S gene of the B.1 strain contains a single mutation from A to G at nucleotide 1841, resulting in a change from aspartic acid to glycine at amino acid position 614 . The RNA sequence of the S gene was folded by using the RNAfold web server , and its secondary structure was displayed by using the Forna platform  (Figure 3A). We also compared the sequences of the S gene between the original B.1 strain and the other eight variants and found that there were 17 mutations, most of which could cause changes to the amino acid sequence. These 17 mutation sites are highlighted in red in Figure 3A.
In theory, as the RNA sequence of the S gene contains 3,822 nucleotides (nt), 3,802 siRNA candidates with a length of 21 nt can be generated. We calculated the PBE values of all siRNA candidates against the S gene and found that they ranged from 2.3 to 46.0. The top four candidate siRNAs with the highest PBE values of 31.1-46.0 (Figure 3B) and outside of the 17 mutation sites were selected (Figure 3A).
In the detailed secondary structures of these siRNAs, each nucleotide is colored from steel blue to white, which indicates an increased probability of hydrogen bond formation between the siRNA and target RNA (Figure 3C-F). Among these siRNAs, 3329i and 1878i had higher PBE values, and their targeting regions did not encode any functional domains of the S protein (Figure 3A-D). The siRNA 1104i targets loop 13, which encodes a small region of the receptor binding domain (RBD) (Figure 3A and E). The siRNA 2351i targets loop 5, which partially encodes the fusion peptide (FP) of the S protein (Figure 3A and F).
To determine the gene silencing effect of the siRNA candidates against the S gene, we first showed that this viral gene was exogenously expressed in human HUVECs and A549 cells at the level comparable to the high expression of a housekeeping gene (GAPDH). The qPCR results showed that the expression of S gene was not detectable (ND) in these two types of cells without transfecting the viral gene (Figure 3G and H).
Calculation of the probability of binding efficiency (PBE) based on single-stranded information of RNA molecules. A. PBE was calculated by summing the hydrogen bonds formed between siRNA and the target region of RNA molecules. NHB refers to the number of hydrogen bonds formed by each nucleotide base pair. B. The distribution of PBE values for three viral genes: spike (S), nucleocapsid (N) and membrane (M) genes.(Click on the image to enlarge.)
Selection of high PBE siRNAs for the spike (S) gene. A. The secondary structure of the RNA molecule of the S gene. The nucleotide sequences corresponding to functional protein domains are labeled in different colors. Mutated sites among the nine viral strains are labeled in red. The loops numbered 5, 13, 17, 25 and 26 are highlighted in orange, as these regions were selected for siRNA targeting. B. The sequence, GC content, PBE and RNA structure of the four selected siRNAs. C-F. Detailed structures of targeted regions by four siRNAs, 3329i (C), 1878i (D), 1104i (E) and 2351i (F), on the S gene with the siRNA sequence outlined in orange. The color of each nucleotide reflects the possibility of being single-stranded. G-H. Real-time qPCR was used to measure the relative expression of the S gene in both HUVECs (G) and A549 cells (H) after transfection with the plasmid DNA of the S gene and individual siRNAs at 24 hours. Cells that were not transfected were used as a background control. All values were normalized to the siNC group, in which cells were transfected with the S gene and a negative control siRNA simultaneously. In No siRNA group, only plasmid DNA was transfected.(Click on the image to enlarge.)
To test the silencing efficiency of the designed siRNAs, plasmid DNA expressing the S gene and individual siRNAs (10 nM) were transfected into HUVECs and A549 cells using Lipofectamine 2000 for 6 hours. The cells were then incubated for 18 hours before RNA isolation and quantification. The results showed that although the S gene was expressed at very high levels in the cells, all four siRNAs significantly reduced its mRNA levels in both HUVECs and A549 cells (Figure 3G and H). Specifically, in HUVECs, the 3329i siRNA, which had the highest PBE value, produced the greatest level of gene silencing of 70%, while the 1878i and 2351i siRNAs reduced the levels of S gene expression by over 50% (Figure 3G). In A549 cells, although siRNA 3329i reduced S gene expression by 55%, the other three siRNAs achieved similar silencing effects ranging from 48-53% (Figure 3H). The results suggest that the PBE value is useful to select candidate regions within RNA for designing siRNAs.
Designed siRNAs can effectively reduce the expression of the nucleocapsid gene
The nucleocapsid (N) gene is conserved among coronaviruses, and the encoded protein is important for viral replication and eliciting the immune response for SARS-CoV-2 [31,32]; thus, inhibiting N gene amplification is a good strategy to prevent viral propagation. We compared the sequences of the N gene among the nine prevalent viral strains and found 13 mutations. The locations of these mutated nucleotides are indicated in red in the secondary structure of the N gene with 1,260 nucleotides (Figure 4A). We then chose the top four siRNAs that had the highest PBE values of 33.3-41.0 among the possible 1,240 siRNA candidates (Figure 4B-F) and compared their silencing effects in the two cell lines that overexpressed the N gene (Figure 4G and H). In HUVECs, 418i siRNA, which had the highest PBE value, produced the best silencing effect, with a reduction of N gene expression at mRNA level of 64%, and the other three siRNAs generated slightly lower but significant silencing effects of 55-61% (Figure 4G). In A549 cells, the top two ranked siRNAs in terms of PBE values, 418i and 1068i, achieved a better silencing effect than the other two siRNAs (Figure 4H). Taken together, the top three ranked siRNAs reduced N gene expression by over 50% in both types of cells.
Designed siRNAs can effectively reduce the expression of the membrane gene
The membrane (M) gene encodes the M glycoprotein that is located at the membrane of the viral particle. It is the most abundant structural protein of the SARS-CoV-2 virus . By comparing the viral sequences of the M gene among the nine strains, we found that only four point mutations occurred in three of the viral strains. None of these mutations caused any changes at the amino acid level (Figure 5A).
The secondary structure of the M gene was generated, and ss-count information was used to calculate the PBE values covering the entire length of M gene with 669 nucleotides (Figure 5B). Based on the PBE values of 649 candidates, we selected the top three siRNAs (Figure 5C), and the detailed secondary structures of these siRNA-targeted regions are shown in Figure 5D-F. The results showed that both types of cells expressed high levels of the M gene after the transfection of plasmid DNA for the M gene (Figure 5G and H). However, when molecules of siRNA 607i were transfected together with the plasmid DNA of the M gene, they reduced M gene expression by 39-43% in both HUVECs and A549 cells (Figure 5G and H). The other two siRNAs reduced M gene expression by 39-47% in A549 cells but only decreased M gene expression by 26-30% in HUVECs (Figure 5G and H).
The S and N proteins are the major focuses for developing antiviral drugs and vaccines for SARS-CoV-2 [31,35], especially the D614G form of the S protein, which has been suggested to enhance viral fitness and has become prevalent in major viral strains . In this study, we used siRNA to decrease the levels of viral RNAs, which has several advantages over other methods.
First, siRNAs directly exert silencing effects and act faster than other types of RNAi, such as shRNA or microRNA (miRNA), as the latter two need to be intracellularly expressed and processed into mature forms before they can bind with the target genes.
Second, the design of siRNAs is flexible for highly mutable targets such as the RBD of the S gene . In this study, we avoided mutated regions identified among the nine prevalent strains. In the future, when more mutations emerge in the viral genome, we can modify our siRNAs to target these new viral strains. Furthermore, we can introduce a pool of siRNAs to target multiple variants at the same time. In this manner, the development of siRNA-based preventative or curative medication can readily keep pace with the evolution of SARS-CoV-2.
Third, as a noninvasive approach, inhalation of liposome-encapsulated siRNAs can theoretically reduce the risk of nasal infection to SARS-CoV-2 without bringing about severe side effects . In addition, such an approach more easily reaches a relatively high drug concentration in the cells lining the respiratory tract, and an overall lower drug dosage is required to achieve a protective purpose.
Although there has been no consensus on how many viral particles are needed for an effective infection, it is widely observed that a higher viral load positively correlates with a higher risk of transmission  and increased disease severity . A cohort study reported that in the respiratory tract of patients with severe COVID-19 there were approximately 7×106 viral copies per milliliter of specimens . In our study, we transfected cells with 30 ng of plasmid DNA corresponding to 3.3-4.9×109 DNA molecules per milliliter of culture medium in a culture well containing 1×105 cells; thus, more than 1×104 DNA molecules were estimated to enter each transfected cell. Under the cytomegalovirus promoter known to have a strong transcriptional effect, many more RNA molecules were synthesized, which resulted in a high abundance of expressed viral genes with a threshold cycle (Ct) number comparable to that of a highly expressed housing keeping gene, GAPDH, in the host cells. Even at such high concentrations of viral genes, the siRNAs selected based on high PBE values significantly reduced the expression of the viral genes within 24 hours. From this point of view, we expect our siRNAs to produce better effects when applied clinically, as the RNA copy number is usually lower in real infections .
Selection of high PBE siRNAs for the nucleocapsid (N) gene. A. The secondary structure of the RNA molecule of the N gene where the color of each nucleotide reflects the possibility of being single-stranded. The mutated sites among the nine viral strains are labeled in red. The loops numbered 5, 8, 10 and 13 are highlighted in orange, as these regions were selected for siRNA targeting. B. The sequence, GC content, PBE and RNA structure of the four selected siRNAs. C-F. Detailed structures of targeted regions by four siRNAs, 418i (C), 881i (D), 214i (E) and 1068i (F), on the N gene where the siRNA sequences are outlined in orange. G-H. Real-time qPCR was used to measure the relative expression of the N gene in both HUVECs (G) and A549 cells (H) after transfection with the plasmid DNA of the N gene and individual siRNAs at 24 hours. Cells that were not transfected were used as a background control. All values were normalized to the siNC group, in which cells were transfected with the N gene and a negative control siRNA simultaneously.(Click on the image to enlarge.)
Selection of high PBE siRNAs for the membrane (M) gene. A. Mutation information of the M gene among the nine viral strains. B. The secondary structure of the RNA molecule of the M gene where the mutated sites are labeled in red. The loops numbered 2, 5 and 9 are highlighted in orange, as these regions were selected for siRNA targeting. C. The sequence, GC content, PBE and RNA structure of the three selected siRNAs. D-F. Detailed structures of targeted regions by three siRNAs, 607i (D), 344i (E) and 79i (F), on the M gene with the siRNA sequences outlined in orange. G-H. Real-time qPCR was used to measure the relative expression of the M gene in both HUVECs (G) and A549 cells (H) after transfection with the plasmid DNA of the M gene and individual siRNAs at 24 hours. Cells that were not transfected were used as a background control. All values were normalized to the siNC group in which cells were transfected with the M gene and a negative control siRNA simultaneously.(Click on the image to enlarge.)
In summary, we used the PBE value and the secondary structure of three key viral genes to guide the design of candidate siRNAs. By using an in vitro viral expression and targeting model, we showed that most of the selected siRNAs achieved a silencing effect of over 50%. In response to the progression of the COVID-19 pandemic and the constant evolution of SARS-CoV-2, we proposed a flexible and promising siRNA-based approach for viral prevention as well as early treatment. For the application prospects in the future, this approach can complement the preventative medication for populations that are unrecommended for vaccination, as well as persons who need a quick and reinforced prevention to viral infection.
This work was supported by grants from the Science and Technology Development Fund (FDCT), Macao SAR, China (No. 083/2016/A2 and 068/2017/A2). We gratefully acknowledge the authors listed in Supplementary Table 1 for sharing the genetic sequence of viral strains via the GISAID Initiative. We would like to thank Dr. Hao Jia from our laboratory for his careful reading of this manuscript.
R. Wu and K.Q. Luo designed the research study. R. Wu conducted the experiments. R. Wu and K.Q. Luo analyzed the data and wrote the manuscript.
The authors have declared that no competing interest exists.
1. Roser M, Ritchie H, Esteban O-O, Joe H. Coronavirus Pandemic (COVID-19). Our World Data. 2020 https://ourworldindata.org/coronavirus
2. Machhi J, Herskovitz J, Senan AM, Dutta D, Nath B, Oleynikov MD. et al. The Natural History, Pathobiology, and Clinical Manifestations of SARS-CoV-2 Infections. Journal of Neuroimmune Pharmacology. 2020;15:359-386
3. Hamming I, Timens W, Bulthuis MLC, Lely AT, Navis GJ, van Goor H. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. Journal of Pathology. 2004;203:631-637
4. Kim D, Lee JY, Yang JS, Kim JW, Kim VN, Chang H. The Architecture of SARS-CoV-2 Transcriptome. Cell. 2020;181:914-921.e10
5. Denison MR, Graham RL, Donaldson EF, Eckerle LD, Baric RS. Coronaviruses: an RNA proofreading machine regulates replication fidelity and diversity. RNA Biology. 2011;8:270-279
6. Callaway E. The coronavirus is mutating - does it matter?. Nature. 2020;585:174-177
7. WHO. SARS-CoV-2 variants. 2020. www.who.int/csr/don/31-december-2020-sars-cov2-variants/en (accessed 31 Dec 2020).
8. Asselah T, Durantel D, Pasmant E, Lau G, Schinazi RF. COVID-19: Discovery, diagnostics and drug development. Journal of Hepatology. 2021;74:168-184
9. Callaway E. The race for coronavirus vaccines: a graphical guide. Nature. 2020;580:576-577
10. Davis MM, Zickafoose JS, Halvorson AE, Patrick SW. Parents' Likelihood to Vaccinate Their Children and Themselves against COVID-19. medRxiv. 2020. 2020 11.10.20228759
11. Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC. Potent and specific genetic interference by double-stranded RNA in caenorhabditis elegans. Nature. 1998;391:806-811
12. Hammond SM, Bernstein E, Beach D, Hannon GJ. An RNA-directed nuclease mediates post-transcriptional gene silencing in Drosophila cells. Nature. 2000;404:293-296
13. Kelleher AD, Cortez-Jugo C, Cavalieri F, Qu Y, Glanville AR, Caruso F. et al. RNAi therapeutics: an antiviral strategy for human infections. Current Opinion in Pharmacology. 2020;54:121-129
14. Agrawal N, Dasaradhi PVN, Mohmmed A, Malhotra P, Bhatnagar RK, Mukherjee SK. RNA Interference: Biology, Mechanism, and Applications. Microbiology and Molecular Biology Reviews. 2003;67:657-685
15. Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC. Potent and specific genetic interference by double-stranded RNA in caenorhabditis elegans. Nature. 1998;391:806-811
16. Luo KQ, Chang DC. The gene-silencing efficiency of siRNA is strongly dependent on the local structure of mRNA at the targeted region. Biochemical and Biophysical Research Communications. 2004;318:303-310
17. Fu A, Ma S, Wei N, Tan BXX, Tan EY, Luo KQ. High expression of MnSOD promotes survival of circulating breast cancer cells and increases their resistance to doxorubicin. Oncotarget. 2016;7:50239-50257
18. Ma S, Fu A, Lim S, Chiew GGY, Luo KQ. MnSOD mediates shear stress-promoted tumor cell migration and adhesion. Free Radical Biology and Medicine. 2018;129:46-58
19. Yin Q, Flemington EK. siRNAs against the Epstein Barr virus latency replication factor, EBNA1, inhibit its function and growth of EBV-dependent tumor cells. Virology. 2006;346:385-393
20. Wooddell CI, Yuen MF, Chan HLY, Gish RG, Locarnini SA, Chavez D. et al. Rnai-based treatment of chronically infected patients and chimpanzees reveals that integrated hepatitis b virus DNA is a source of hbsag. Science Translational Medicine. 2017;9:eaan0214
21. Delville M, Touzot F, Couzin C, Hmitou I, Djerroudi L, Ouedrani A. et al. Safety of CD34 + Hematopoietic Stem Cells and CD4 + T Lymphocytes Transduced with LVsh5/C46 in HIV-1 Infected Patients with High-Risk Lymphoma. Molecular Therapy - Methods and Clinical Development. 2019;13:303-309
22. Li BJ, Tang Q, Cheng D, Qin C, Xie FY, Wei Q. et al. Using siRNA in prophylactic and therapeutic regimens against SARS coronavirus in Rhesus macaque. Nature Medicine. 2005;11:944-951
23. Setten RL, Rossi JJ, Han S ping. The current state and future directions of RNAi-based therapeutics. Nature Reviews Drug Discovery. 2019;18:421-446
24. Shu Y, McCauley J. GISAID: Global initiative on sharing all influenza data - from vision to reality. Eurosurveillance. 2017;22:30494
25. Elbe S, Buckland-Merrett G. Data, disease and diplomacy: GISAID's innovative contribution to global health. Global Challenges. 2017;1:33-46
26. Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C. et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. 2018;34:4121-4123
27. Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Research. 2003;31:3406-3415
28. Plante JA, Liu Y, Liu J, Xia H, Johnson BA, Lokugamage KG. et al. Spike mutation D614G alters SARS-CoV-2 fitness. Nature. 2020
29. Gruber AR, Lorenz R, Bernhart SH, Neuböck R, Hofacker IL. The Vienna RNA websuite. Nucleic acids research. 2008;36:W70-W74
30. Kerpedjiev P, Hammer S, Hofacker IL. Forna (force-directed RNA): Simple and effective online RNA secondary structure diagrams. Bioinformatics. 2015;31:3377-3379
31. Dutta NK, Mazumdar K, Gordy JT. The Nucleocapsid Protein of SARS-CoV-2: a Target for Vaccine Development. Journal of Virology. 2020 94
32. Grifoni A, Sidney J, Zhang Y, Scheuermann RH, Peters B, Sette A. A Sequence Homology and Bioinformatic Approach Can Predict Candidate Targets for Immune Responses to SARS-CoV-2. Cell Host and Microbe. 2020;27:671-680.e2
33. Mousavizadeh L, Ghasemi S. Genotype and phenotype of COVID-19: Their roles in pathogenesis. Journal of Microbiology, Immunology and Infection. 2020 in press
34. Thomas S. The structure of the membrane protein of sars-cov-2 resembles the sugar transporter semisweet. Pathogens and Immunity. 2020;5:342-363
35. Huang Y, Yang C, Xu X feng, Xu W, Liu S wen. Structural and functional properties of SARS-CoV-2 spike protein: potential antivirus drug development for COVID-19. Acta Pharmacologica Sinica. 2020;41:1141-1149
36. Li Q, Wu J, Nie J, Zhang L, Hao H, Liu S. et al. The Impact of Mutations in SARS-CoV-2 Spike on Viral Infectivity and Antigenicity. Cell. 2020;182:1284-1294.e9
37. Chow MYT, Qiu Y, Lam JKW. Inhaled RNA Therapy: From Promise to Reality. Trends Pharmacological Sciences. 2020;41:715-729
38. Shrestha NK, Marco Canosa F, Nowacki AS, Procop GW, Vogel S, Fraser TG. et al. Distribution of Transmission Potential During Nonsevere COVID-19 Illness. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America. 2020;71:2927-2932
39. Fajnzylber J, Regan J, Coxen K, Corry H, Wong C, Rosenthal A. et al. SARS-CoV-2 viral load is associated with increased disease severity and mortality. Nature Communications. 2020;11:1-9
40. To KKW, Tsang OTY, Leung WS, Tam AR, Wu TC, Lung DC. et al. Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by SARS-CoV-2: an observational cohort study. The Lancet Infectious Diseases. 2020;20:565-574
41. Chang MC, Hur J, Park D. Interpreting the COVID-19 Test Results: A Guide for Physiatrists. American Journal of Physical Medicine and Rehabilitation. 2020;99:583-585
Corresponding author: kluoedu.mo. Faculty of Health Sciences, University of Macau, Macao SAR, China