Coronaviruses are a large family of viruses that were first described over 50 years ago (1) . Since the turn of the millennium, there have been two major global outbreaks caused by coronaviruses, namely SARS-CoV in 2003 and MERS-CoV in 2012 (2) . The ongoing COVID-19 pandemic caused by SARS-CoV-2 represents the third and most devastating of these outbreaks. These outbreaks, notably the COVID-19 pandemic, are harsh reminders of the challenges posed by emerging infectious diseases. The global impact of the COVID-19 pandemic has brought to the forefront the need to rapidly develop and deploy an effective vaccine. However, despite ongoing concerted research efforts, it is accepted that such a vaccine will not be available in time to contribute to the containment of the ongoing pandemic. Current management guidelines include the use of repurposed drugs such as chloroquine and its analog hydroxychloroquine as well as antiviral agents (3) . However, the need for well-designed clinical trials to validate their efficacy continues to be highlighted. To effectively address the ongoing COVID-19 pandemic, there is a recognized need for a framework for rapid identification of novel targets for diagnostic and therapeutic interventions as well as determine clinically approved drugs with high potential for repurposed use against SARS-CoV-2. Publicly available transcriptomic datasets generated from SARS-CoV infected humans, and mammalian cells represent a wealth of data that could be used to identify consistent differentially expressed genes, which could then be validated against SARS-CoV-2 infected epithelial cells transcriptomic datasets. A comprehensive toxicogenomic analysis of the identified genes could potentially identify possible interactions with clinically proven drugs. This simple approach can be used for the rapid identification of novel targets and drugs for further validation. In this study, we have applied this approach, and our findings have identified IFITM3 as an early upregulated gene and indicate that valproic acid enhances IFITM3 mRNA expression and antiviral action.
Publicly available transcriptomic datasets were retrieved from Gene Expression Omnibus (GEO) 
Mouse lung tissue transcriptome response to a mouse-adapted strain of SARS-CoV in wild type C57BL6/NJ mice and TLR3-/-mice C57BL6/NJ MA15 GSE68820 (11) (https://www.ncbi.nlm.nih.gov/geo/). Only microarray gene expression datasets with the word "SARS-COV, " virus, or modified strain infected vs. mock-infected and no more than 48 h after the infection. Twelve datasets fulfilled the criteria, as detailed in Table 1 . We used GEOquery and limma R packages through the GEO2R tool for each dataset (12) . After sorting the genes according to the False Discovery Rate (FDR), the top 2,000 differentially expressed probes with FDR <0.05 were selected from each dataset. The annotated genes of the 5,000 probes in each dataset were intersected with differentially expressed genes (DEGs) from all other datasets. The DEGs that were common in at least 9 out of the 12 (75%) datasets were identified as shared genes that are consistently DEG in the first 48 h of SARS-COV infection. Enriched Ontology Clustering for the identified genes was performed to explore using the Metascape (http://metascape.org/gp/index.html#/main/step1). The shortlisted genes expression was then explored in another dataset (GSE147507), where RNA-Sequencing of primary human lung epithelium (NHBE) mock-treated or infected with SARS-CoV-2 was done to examine whether there is a difference in the response of SARS-CoV-2 from other strains in terms of DEGs (13) .
In total, 9,692 genes were differentially expressed genes (DEGs) between mock-infected and virally infected models in the 12 studies. Thirty-eight genes that were DEGs in at least 9 out of 12 studies (75%) were considered common DEGs due to SARS-COV infection of the lungs in the first 48 h post-infection. These genes are listed in Table 2 .
In order to identify DEG in SARS-CoV infected lung tissuespecific to each of the models used and those which are shared, we intersected the DEGs from datasets that use the same model. Human epithelial cells datasets (GSE17400, GSE33267, GSE37827, GSE47960, and GSE48142), mice datasets (GSE33266, GSE50000, GSE52920, GSE64660, and GSE68820), Ferret (GSE11704) and Cynomolgus macaques (GSE23955) were all intersected with the COVID-19 infected epithelial cells dataset as shown in Figure 1 . The number of DEG intersected between different species is listed in the Table 3 . Epithelial cells infected with SARS-CoV-2 shared 9 DEGs (MX1, OAS3, XAF1, IFI44, MX2, IRF7, STAT1, IFIT3, and IFIT1) with Human Lung Cells, Mice, and Cynomolgus maca.
The Identified Genes Are Involved in the Immune Response Against RNA Viruses As expected, the top genes identified are involved in innate immune responses against RNA viruses. These include the cytosolic DNA-sensing pathway, Toll-like receptor signaling pathway, and negative regulation of binding. Interferon (IFN) response to viral infections such as type I interferon signaling pathway, defense response to the virus, the antiviral mechanism by IFN-stimulated genes, regulation of type I interferon production, response to interferon-alpha, and regulation of defense response to virus and Influenza A, were also upregulated. Genes that play significant roles in activating immune systems such as regulation of response to cytokine stimulus, negative regulation of immune response, myeloid cell homeostasis, and positive regulation of the multi-organism process are also upregulated (Figure 2 and Table 4 ).
The identified genes expression levels were higher in human bronchial epithelium infected with SARS-CoV-2 compared to those mock-infected (Figure 3) . However, only IFITM3 showed a significant difference (p < 0.05), while two other genes OAS2 and MX1 showed a trend of enhancement, although it was not statistically significant (p = 0.06 using the two-stage linear stepup procedure of Benjamini, Krieger, and Yekutieli). IFITM3 mRNA levels were one of the highly expressed genes compared to the other identified genes at baseline in mock-infected HBE and were further induced by the virus, which results in overall high mRNA levels.
Valproic Acid Can Upregulate the IFITM3 mRNA Expression
Next, we searched the Comparative Toxicogenomics Database (http://ctdbase.org/) to identify drugs/chemicals that might affect the mRNA expression of IFITM3 in at least two reference studies (14) . Interestingly valproic acid, carbon nanotubes, nickel, and tert-butylhydroperoxide were shown to upregulate IFITM3 expression while pirinixic acid, acetaminophen, and Ethinyl estradiol decreased such an expression ( Table 5) .
In order to examine the effect of valproic acid therapy on the IFITM3 mRNA expression in immune cells of the blood, a publicly available transcriptomics dataset (GSE143272) was extracted and reanalyzed. Healthy controls were compared to responders and non-responders patients on valproic acid therapy. We found upregulation of the mRNA expression of IFITM3 in patients, and the difference was significant in the responder group only (p < 0.05) compared to healthy controls (Figure 4) . 
In response to viral RNAs, like in the case of SARS-CoV-2, the innate immune system will unleash interferon (IFN), to activate antiviral mechanisms and effector cells like natural killers (15) . In mice infected with SARS-CoV, a delayed and prolonged type I interferon (IFN-I) signaling leads to lung immunopathology as it promotes the accumulation of pathogenic inflammatory cells with increased lung cytokine/chemokine levels and vascular leakage (16) . This prolonged IFN-I and virally induced IL-10 set the scene for secondary bacterial infection, which can add a strong IL1β and TNFα-mediated inflammatory response to magnify lung damage (17) . Understanding how SARS-CoV-2 can manipulate IFN is vital in deciphering the battle of the body against the viral spread and consequence. Our reanalysis of transcriptomic data showed that although the IFN pathway is upregulated consistently in SARS-CoV related infection, SARS-CoV-2 showed specific upregulation of the gene for a unique interferon-induced protein, namely IFITM3. IFITM3 is a 15-kDa protein that localizes to endosomes and lysosomes and is possibly acquired by mammalian ancestral cells via horizontal gene transfer (18) .
Interferon-induced transmembrane proteins (IFITMs 1, 2, and 3) are innate immune responders to virus infections as they regulate the fusion of invading virus and endocytic vesicles and direct it to the lysosomes (19, 20) . IFITM3 can further alter membrane rigidity and curvature to inhibit virus membrane fusion (21) . Such action is important to prevent the release of viral particles into the cytoplasm, which controls viral spread (22) . During influenza A infection of human airway epithelial cells, IFITM3 was shown to clusters on viruscontaining endosomes and lysosomes within few hours postinfection, indicating its role in the early phase of viral entry (23) . Even platelets and megakaryocytes were shown to remarkably upregulate IFITM3 to prevent viral progression during influenza infection (24) . BST2, IFI35, IFIT2, IFIT1, IFIT3, IRF7, ISG20,  MX1, MX2, OAS2, OAS3, SP100, STAT1,  ISG15, IFITM3, USP18, XAF1, ZBP1, RSAD2,  TRIM21, DDX58, IL6, CXCL10, CXCL11 The epithelial cell and resident leukocytes in lung upper and lower airways that constitutively express IFITM3 can withstand viral infections, and this is vital to decide viral tropism as viruses favor cells with low IFITM3 expression (25) . IFITM3 enhances the accumulation of CD8+ T cells in airways to promote mucosal immune cell persistence (26) . Lung and circulating immune cells were reported to express less IFITM3 than other tissues, and this was a suggestive reason for COVID-19 severity and cytokine release syndrome (27) .
Interestingly, IFITM3-rs12252-C/C SNP prevalence in the Chinese population is 26.5%, and recent research confirmed that SNPs in IFITM3 could change the severity of influenza infection, as was shown in one case with COVID-19 (28) . IFITM3 polymorphisms have been linked with hospitalization and mortality during influenza virus infection (29) .
Expressing the gene is not the only prerequisite to the antiviral action of IFITM3, as it was found recently that within the protein, an amphipathic helix is critical for its blocking effect of viral fusion of similar pathogenic viruses like influenza A virus and Zika virus (30). Another factor that regulates the IFITM3 trafficking specificity to such viruses is that it requires S-palmitoylation (19, 20) . S-palmitoylation (S-PALM) is the reversible process of linking a fatty acid chain to cysteine residues of the substrate protein (31) . Multiple zinc finger DHHC domain-containing palmitoyltransferases (ZDHHCs) can palmitoylate IFITM3 to make it a fully functional antiviral protein (32) . It seems that bats (order Chiroptera), which act as natural hosts for many viral infections, use IFITM3 as an antiviral mechanism if there is S-palmitoylation of the protein; however, if this modification is disturbed, the bat can develop viral infection (33) . Based on that, we can suggest that severe COVID-19 cases might be due to either non-functional IFITM3 by SNP, failure of lung cells to upregulate IFITM3 in response to interferon, a mutation in amphipathic helix sequence or modification in S-palmitoylation. Further examination and screening for the IFITM3 dynamics in COVID-19 might explain the possible therapeutic and diagnostic options.
Our toxicogenomic analysis showed that valproic acid increased the mRNA expression of IFITM3, supporting a new report that the SARS-CoV-2-human protein-protein interaction map showed that valproic acid might be a potential repurposing drug for COVID-19 (34) . Virtual screening, docking, binding energy calculation, and simulation show that valproic acid forms stable interaction with nsP12 of CoV and can inhibit its function (14) . Valproic acid is currently used for the treatment of epilepsy and known to target histone deacetylases (HDACs) that modify the gene expression epigenetically (35) . Valproic acid was shown to inhibit mature and fully infectious enveloped viruses release as it alters cellular membrane composition (36) . The modest and broad antiviral activity of valproic acid made the drug an attractive possibility due to its availability, and limited side effects for a short term of use during acute viral disease (37) . The reported side effects like hepatotoxicity and teratogenicity are mainly associated with the parental compound valproate and can be avoided by the use of its derivatives like Valpromide (VPD) and valnoctamide (VCD) . A recent open-label proof-of-concept trial of 10 days Intravenous Valproic Acid for Severe COVID-19 showed a 50% reduction in the case fatality rate and length of stay (38) . More studies are needed to explore the promising potential of valproic acid in the treatment of COVID-19.
One limitation of the study is that it is based on the publicly available transcriptome dataset, which is limited in number, partly because this is a novel disease, but also because ongoing lockdowns have made it challenging for scientists to carry out the extensive laboratory work required.
Our evaluation showed that the analysis of publicly available transcriptomic data could be a reasonable approach to identify the novel target and suggest drugs that can modify the action of such targets during the early phases of emerging infections like COVID-19 until a complete understanding of the disease become clear. This can justify the experimental use of clinically approved drugs and guide the clinicians in their limited options against such lethal disease.
All datasets presented in this study are included in the article/supplementary material.