key: cord-291719-1ku6cmwj
authors: Hajjo, Rima; Tropsha, Alexander
title: A Systems Biology Workflow for Drug and Vaccine Repurposing: Identifying Small-Molecule BCG Mimics to Reduce or Prevent COVID-19 Mortality
date: 2020-10-06
journal: Pharm Res
DOI: 10.1007/s11095-020-02930-9
sha: 
doc_id: 291719
cord_uid: 1ku6cmwj

PURPOSE: Coronavirus disease 2019 (COVID-19) is expected to continue to cause worldwide fatalities until the World population develops ‘herd immunity’, or until a vaccine is developed and used as a prevention. Meanwhile, there is an urgent need to identify alternative means of antiviral defense. Bacillus Calmette–Guérin (BCG) vaccine that has been recognized for its off-target beneficial effects on the immune system can be exploited to boast immunity and protect from emerging novel viruses. METHODS: We developed and employed a systems biology workflow capable of identifying small-molecule antiviral drugs and vaccines that can boast immunity and affect a wide variety of viral disease pathways to protect from the fatal consequences of emerging viruses. RESULTS: Our analysis demonstrates that BCG vaccine affects the production and maturation of naïve T cells resulting in enhanced, long-lasting trained innate immune responses that can provide protection against novel viruses. We have identified small-molecule BCG mimics, including antiviral drugs such as raltegravir and lopinavir as high confidence hits. Strikingly, our top hits emetine and lopinavir were independently validated by recent experimental findings that these compounds inhibit the growth of SARS-CoV-2 in vitro. CONCLUSIONS: Our results provide systems biology support for using BCG and small-molecule BCG mimics as putative vaccine and drug candidates against emergent viruses including SARS-CoV-2. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s11095-020-02930-9) contains supplementary material, which is available to authorized users.

Few months after the declaration of COVID-19 pandemic by the World Health Organization (WHO), the disease-causing virus is still sweeping the globe, causing more fatalities, failing health care systems, and resulting in severe economic losses. Currently there are no approved drugs to treat COVID-19, and new vaccine development is expected to take at least 12-18 months (1,2), with growing fears of possible failure associated with changes in viral antigenic determinants (3) or shortlived immunity (4) . Additionally, the highly specific virusneutralizing antibodies in recovered patients may be short lived and ineffective in preventing the disease caused by the emerging variable strains of the virus (4) . With these uncertainties regarding an eminent specific SARS-CoV-2 vaccine, there is a need to search for current alternatives, such as agents that can stimulate or emulate the unique capabilities of our innate immune system. Recent immuno-oncology success stories indicate that the best cancer-fighting strategies results from unleashing the patients' immune power (5) (6) (7) (8) . There is an increased awareness that harnessing innate immune responses, opens up new possibilities for long-term, multifaceted tumor control (9, 10) and infectious disease prevention (11) (12) (13) . Therefore, next generation antiviral vaccines should be capable of boosting innate immune responses to tackle a wide range of novel pathogens very early after exposure, as single treatments or adjuvants to traditional vaccines targeting the adaptive immune system.

Accumulating evidence from the biomedical literature indicates that SARS-CoV-mediated pathology, a very similar pathology to SARS-CoV-2, was mainly caused by ineffective innate immune responses, associated with a severe reduction in the number of T cells in the blood (14) . Recent evidence indicated that SARS-CoV-2 and Mycobacterium tuberculosis (Mtb) share unique similarities in terms of the host protein interaction partners, and both pathogens infect lung tissues (15) . On the other hand, old 'polypharmacological' vaccines, such as the BCG vaccine for tuberculosis (TB), has shown promising therapeutic effects for a wide range of infectious and non-infectious diseases including bladder cancer (16) (17) (18) . Studies showed that BCG's polypharmacological effects were not limited to memory T cell immunity, but promoted strong, beneficial, and long-lasting effects on innate immunity. The WHO also recognized these beneficial 'offtarget' effects of BCG, calling for a further investigation to repurpose this vaccine for other orphan life-threatening diseases (19) . Indeed, there are multiple clinical trials testing BCG for 216 conditions other than TB including 19 studies for COVID-19 as reported on clinicaltrials.gov (20) . Additionally, few recent peer-reviewed reports have pointed to an epidemiological relationship between BCG and COVID-19 without providing substantial evidence (21) (22) (23) (24) . One should expect that the results of the randomized clinical trials (RCTs) will help establish the value of the BCG vaccine as a treatment or prophylactic against the disease.

Herein, we describe a unique drug and vaccine repurposing workflow, and list high confidence proteins and pharmacological classes of compounds, that work as BCG mimics at the system level by inducing beneficial long lasting trained immune response. We also propose that BCG mimics can be used as alternatives to BCG in protecting from COVID-19 and other emergent infectious diseases.

We have developed and applied a systems biology workflow to study the BCG network pharmacology and prioritize smallmolecule BCG mimics and antivirals. This workflow is based on our original chemocentric informatics workflow described thoroughly in a previous report (25) . Our current workflow ( Fig. 1) incorporates three major components: (1) a module for mining and prioritizing gene signatures representative of a condition or a biological state; (2) a network-mining module to identify genetic perturbations that induce gene expression profiles that are highly enriched with the genes constituting the condition gene signature; and (3) a pathway enrichment module to understand the biological processes involved in the mechanism of action of BCG and highly correlated genetic perturbagens.

A consensus gene signature for BCG vaccine was derived from gene expression profiles in peripheral blood mononuclear cells (PBMCs) in response to a BCG challenge test reported by Matsumiya et al (76) , GSE58636 dataset on NCBI Gene Expression Omnibus (GEO) (27) . We used the data collected from whole blood samples taken from healthy human subjects enrolled in phase 1 trial (clinical trials registration: NCT01194180). For the purposes of this study we used the gene expression profiles generated from two human subject groups included in the above trial: group 1 (BCG naive), and group 2 (BCG vaccinated; median time since vaccination, 10 years). To study network pharmacology and query the connectivity map, we developed a consensus gene signature using genes that showed significant differential gene expression in response to a BCG challenge test (stimulated) in comparison with controls (unstimulated) on days 0 and 14 in both groups 1 and 2.

A systematic search, for nearest neighbor (NN) genes/proteins of the upregulated and downregulated genes in BCG's gene signature, was conducted in Cytoscape (77) version 3.8.0 using the STRING (78) protein query application. All retrieved protein-protein interactions (PPIs), including both physical and functional interactions were retrieved from widely used and reliable databases such as MINT (79) , HPRD (80), BIND (81) , DIP (82), BioGRID (83) , KEGG (84), Reactome (60), EcoCyc (85) , NCI-Nature Pathway Interaction Database (86) , and Gene Ontology (GO) (87) protein complexes. Network building tools in Cytoscape version 3.7.2 were used to generate PPI networks for BCG-CGS.

Enrichment analysis was conducted in Cytoscape (77) and MetaCore to identify pathways and biological processes associated with BCG-CGS and CMap genetic connections. The significance of the enrichment was determined by the hypergeometric test (88) . All terms from the ontology were ranked based on their calculated p values. Ontology terms with p values less than the p value threshold 0.05 are defined as statistically significant and therefore relevant to the studied list of genes. All terms from the ontology were ranked according to their calculated p values.

The CMap (27, 89) is a chemogenomics database that catalogs 1.3 Million profiles of transcriptional responses of human cells to chemical and genetic perturbations. Currently, there are 27,927 perturbagens (19,911 small molecules, and 7494 genetic perturbagens) producing 476,251 expression signatures in 9 human cell lines: PC3, VCAP, A375, A459, HA1E, HCC515, HT29, MCF7, HEPG2. This database of cellular signatures has been produced using the L1000 platform (27) ; a high-throughput gene expression assay that measures the mRNA transcript abundance of 978 "landmark" genes from human cells.

Causal reasoning (90) analysis identifies genes and proteins of a 'topological significance' in order to make decisions whether these genes/proteins are eligible for targeting in the studied phenotype. In this study, we applied causal reasoning to identify molecular regulators that most likely directly cause the observed expression changes in transcriptional profiles in response to BCG. In this approach, changes in gene expression in both directions as well as the effect of edges in the network are taken into account. For each node (i.e., gene) in the causal reasoning network, observed changes in expression are matched with the expected changes inferred from the network structure given the hypothesis that the observed gene expression is decreased or increased due to its activity. Each node has an outgoing activation or inhibition effects on other objects in the knowledge database, and a key hub with a predicted increase in activity shows increased expression for those genes that the hub is known to activate, and it shows decreased expression for genes it is known to inhibit. Each predicted key hub has a prediction P value which is produced as a result of a binomial test used to assess the probability of making a given number of supportive data out of all defined differentially expressed genes (DEGs) in the examined data. It is noteworthy that causal reasoning examines both direct neighbors of differentially expressed genes, and remote (several steps away) regulators. All causal reasoning predictions were performed in Key Pathway Advisor from Clarivate Analytics, using the Pollard method (91).

Gplots (92) v3.0.1.2 was used for plotting enhanced heatmaps for transcriptional data (e.g., heatmap representing BCG-CGS in Fig. 2 ). Heat maps were generated using the heatmap.2 function included in this package.

To study the BCG polypharmacology and potential beneficial effects of this vaccine in preventing the fatal consequences of COVID-19, we have devised and implemented a 'network biology' workflow ( Fig. 1) to interrogate the hypothesis that BCG vaccination may protect from COVID-19 fatalities. This workflow is based on our drug repurposing chemocentric informatics workflow which has been validated previously for small-molecule drug repurposing (25) . The current workflow is tweaked towards vaccine repurposing by employing novel bioinformatic approaches to computationally model and connect molecular networks in an effort to understand the underlying 'network' biology of vaccines, and pinpoint the regulatory genes and proteins responsible for causing the observed beneficial multitherapeutic effects. Although we are not the first group to use network biology approaches to study the transcriptional changes of vaccines, to our knowledge, this is the first study that uses these approaches both to support vaccine repurposing, specifically for COVID-19, as well as identify putative small molecule drugs that can mimic the vaccine effects.

Our workflow starts with the prioritization of a gene signature to study the BCG network pharmacology. First, we derived a consensus gene signature (CGS) for BCG based on GEO's dataset GSE58636 (26) . Details on BCG-CGS signature are found in Table S1 (Supporting Information). Twenty-two differentially-expressed genes across all 4 experiments (2 Groups × 2 time points discussed in Methods) formed BCG's consensus gene signature (BCG-CGS) shown in Fig.  2a .

All 22 genes in BCG-CGS were used as seed nodes to build a protein-protein interaction network for signature genes (Fig.  2b ). Interactions were extracted from STRING database, and high confidence interactions included physical interactions (e.g., binding), functional interactions (e.g., activation, inhibition, catalysis), or gene co-expression. Two types of networks were generated: 1) high-confidence 'core' network restricted to BCG signature genes as network nodes and high confidence (≥0.70) interactions as network edges, and 2) mediumconfidence interaction network obtained from expanding the core network by 20 additional nodes (Fig. 3 ).

Enrichment analysis results performed in Cytoscape, using STRING's protein-protein interactions, indicated that BCG-CGS is enriched in inflammatory cytokines and immune response modulators (Fig. 2b) . Some signature genes are also involved in the negative control of important viral processes (e.g., (FCN1, TNF and CCL3), and others are involved in the response to viral infections (e.g., IFNG, RNASE6, IL6 and TNF). The complete lists of enriched pathways are included in Tables S2 and S3 (Supporting Information).

We identified 291 key hubs using the causal reasoning method, which seeks to identify molecular regulators that can directly cause the observed transcriptional changes in response to BCG vaccination. Key regulators can be transcriptional factors and proteins with potentially altered activity that explains the transcriptional changes. Top five statistically significant inhibited key hubs included HEY1, DSIPI (GILZ), Jagged1, HAND1 and miR-129-1-3p, whereas top five statistically significant activated key hubs were PHF20, TAFII70, Glutaredoxin, RUNX2 and NOTCH1 (NICD). Top 30 causal key hubs are shown in Table I and all identified 291 key hubs are included in Table S4 (Supporting Information) .

In order to identify experimentally validated upstream regulators that cause transcriptional changes similar to those induced by BCG, we queried the Connectivity Map (CMap) (27) database of the Broad Institute with BCG-CGS and identified proteins and small-molecule drugs that have strong connectivity scores with BCG (Fig. 1) . The CMap approach enabled us to compare BCG-CGS with 'experimentally' predefined signatures of therapeutic compounds and genetic perturbations (i.e., over expression or knockdown) included in the CMap and ranked according to a connectivity scores (ranging from +100 to −100), representing relative similarity to BCG-CGS. The connectivity score itself is derived using a nonparametric, rankbased, pattern-matching strategy based on t he Kolmogorov-Smirnov statistics (28) . All instances in the database are then ranked according to their connectivity scores with BCG-CGS; those at the top (+) are most strongly correlated to the query signature and looked at as BCG mimics, and those at the bottom (−) are most strongly anticorrelated and can reverse BCG's gene signature.

Our analysis identified three highly enriched classes of genetic knockdown (KD) perturbagens and one pharmacological class of drugs that have positive connectivity scores in alveolar A549 cells (i.e., caused similar transcriptional changes to those induced by BCG in alveolar A549 cells). These hits can be considered as BCG mimics capable of inducing transcriptional changes similar to those caused by the BCG vaccine. Therefore, we suggest that BCG mimics can be used as alternatives to BCG vaccination to promote long-lasting b a beneficial effects on immune cells. The three enriched protein classes are: protein phosphatases (with best positive connection for PPP4C KD), histone deacetylases (with best positive connection for HDAC10 KD followed by HDAC11 KD), and mediator complex proteins (with best positive connection for MED6 KD followed by MED7 KD). Additionally, protein kinase C (PKC) activators were enriched as a drug class; and top three PKC activators with highest CMap connectivity scores to BCG-CGS were prostratin, phorbol-12-myristate-13-acetate, and ingenol. It is evident that all of the above four classes of proteins share one common feature: they participate in the transcriptional and metabolic regulation of immune cells in response to environmental cues including responses to pathogens (29) (30) (31) (32) . All top-scoring PKC activators from the CMap, are also known to have antiviral effects or affect T cell activation (33) (34) (35) (36) (37) .

Remarkably, analyzing top ten CMap positive connections with BCG-CGS obtained from nine cell lines indicated that two compounds are approved antiviral drugs: raltegravir (top 3rd positive connection, an HIV integrase inhibitor) and lopinavir (top 6th positive connection, an HIV protease inhibitor). More interestingly, emetine (top 4th positive connection) and lopinavir were recently validated to inhibit SARS-CoV-2 replication in vitro (38) . We also found evidence in the biomedical literature indicating that MST-312 (39), narciclasine (40) and verrucarin-a (41) possess antiviral activities. All CMap hits are provided in Tables S5 and S8 (Sup porting Information).

Initial reports from clinical studies evaluating the use of lopinavir in COVID-19 patients showed that the unbound lopinavir concentrations in the lungs were calculated to be subtherapeutic against SARS-CoV-2 (42, 43) . Another study found that the unbound drug concentrations of lopinavir are far from reaching the EC50 of SARS-CoV-2 (16.4 μg/mL), although they clearly suffice to inhibit HIV-1 (44) . The authors mentioned Fig. 3 High-confidence expanded network for BCG-CGS. Nodes are color-coded using a split pie chart coloring scheme indicating pathway/gene set contribution to each node from the top 5 most enriched pathways/gene lists. Core network is composed of genes in the BCG-CGS that are not singletons.

Step 1 expansion, added 10 additional nodes (i.e., genes) to the core network.

Step 2 expansion, added another 10 nodes for the first expansion.

Step 3 expansion, added another 10 nodes to the second expansion. Expansions were performed to see which pathways remained most statistically significant, and therefore are considered high confidence pathways.

that approximately 60-to 120-fold higher concentrations than those found in COVID-19 patients treated with lopinavir-ritonavir, are required to reach the assumed EC50 at trough levels, making effective treatment of COVID-19 with lopinavir and ritonavir at the currently used doses unlikely (44) .

In order to prioritize high confidence BCG genetic mimics, we integrated hypotheses derived independently from the CMap with those predicted by causal reasoning, and accepted common hits only (i.e., CMap positive connections with BCG-CGS that are also predicted as beneficial drug targets by causal reasoning). This analysis resulted in 30 high confidence common hits reported in Table S9 (Supporting Information) .

We tested whether BCG-CGS, CMap positive connections, or predicted key hubs will have any impact on COVID-19 by identifying overlaps with SARS-CoV-2 interactome, i.e., human proteins that were experimentally validated to interact with SARS-CoV-2 and extracted from two recent reports (45, 46) . This analysis (Fig. 4a ) validated 3 protein hits to have physical links to SARS-CoV-2. The three proteins are transcribed by BRD4, PRKACA and SIRT5; they all were positive connections from the CMap, predicted as statistically 3 † Predicted activity of the key hub by causal reasoning is denoted byif the hub is inhibited, and denoted by + if the hub is activated ‡ Correct/total network predictions: correct for the genes in the dataset predicted correctly; total for the total number of genes in the causal reasoning network § Calculation distance: Using causal reasoning one-step key hubs are defined as statistically significant transcriptional factors that are associated with experimental differential expressed genes regulation. Two-step and three-step key hubs are distant key hubs that regulate one-step transcriptional factors *P-value calcualted for the polynomial test

(2020) 37:212

significant key hubs, and were also validated as SARS-CoV-2 interacting proteins (15) . Additionally, 14 high-confidence CMap positive connections, were validated to make physical interactions with SARS-CoV-2 proteins. These proteins are: PSEN2, PABPC1, HMOX1, CIT, PLAT, IGF2R, RIPK1, NDUFS3, NDUFA5, GGH, NEU1, SCARB1, CSNK2B, F2RL1. And two positive connections, MARK2 and MARK3, were reported to have interactions with corona viruses (45) . Predicted causal key hubs, SIGMAR1 and GNB1, were also validated to have physical links to SARS-CoV-2 (15) , and a third key hub PPIA was known as a human protein interacting with proteins from corona viruses (45) .

Additionally, we mined the biomedical literature to identify evidence for linking BCG small molecule mimics with SARS-CoV-2, corona viruses or viral infections in general. We found that two out of ten top positive compound connections (emetine and lopinavir), were recently validated to inhibit SARS-CoV-2 replication in vitro (38) . Other compounds we found to inhibit the growth of corona viruses, or had general antiviral activities (Table II) .

Previous peer-reviewed reports indicated that BCG's nonspecific effects on the immune system, can reduce all-cause child mortality (47) , protect individuals from numerous viral infections (48) (49) (50) (51) (52) , and it can even enhance the efficiency of some viral vaccines (53) (54) (55) . Recently, several peer-reviewed studies have pointed to a striking correlation between universal BCG vaccination policies and reduced COVID-19 mortality (23) . However, most epidemiological studies identified this correlation without acknowledging other important study confounders like social, economic, and demographic differences between countries. Lately, Escobar et al. mitigated multiple confounding factors for the first time and still observed several significant associations between BCG vaccination and reduced COVID-19 deaths (21) . The authors of this study highlighted the need for mechanistic studies behind the effect of BCG vaccination on COVID-19, and for clinical evaluation of the effectiveness of BCG vaccination to protect from severe COVID-19. Earlier studies suggested that the documented beneficial off-target effects of BCG in protecting from non-TB infections, including perhaps COVID-19, involve a potentiation of innate immune responses through epigenetic mechanisms (56) (57) (58) .

To our knowledge, we report here on the first study providing a mechanistic insight to explain the relationships between BCG and COVID-19 at the molecular and systems biology levels as well as extend this insight toward proposing several BCG mimetics among known drugs as candidates for repurposing against the disease. Our results indicate that BCG-CGS, key regulatory hubs and BCG-mimics identified from the CMap enrich common biological pathways important for key viral processes such as RNA synthesis and processing, virus-host interactions, positive regulation of viral genome replication, and they are also important for the immune response mounted against the virus. Supporting evidence from the biomedical literature confirms that BCG has many beneficial 'off-target' effects that can protect humans from emerging novel pathogens by boasting their innate immune responses (59) . Our studies suggest that BCG promotes a wide-range of transcriptional and metabolic changes, including beneficial gene commensalism, that have been shown to reduce mortality and morbidity from non-TB infectious diseases (51, 52) . We show that BCG can produce these protective 'off target' effects mainly by increasing the production of thymus-generated short-lived undifferentiated CD4+ cells known as naive T cells (Th 0 ), and triggering their differentiation into the long-lived mature naive T cells (MNTs), such as CD4+ and CD8+ T cells (60) . Interestingly, a very recent study published in Science (61) showed that many unexposed patients (20-50%) carry selective and cross-reactive SARS-CoV-2 T cell epitopes protecting patients against severe infection. Another recent study in Cell (62) reported similar observation of strong SARS-CoV-2 selective memory T cell immunity (reminiscent of the functional patterns observed after successful vaccine immunizations) in patients with asymptomatic or mild infections. Although these studies make no connection to any previous BCG vaccination as a source of selective epitopes, we observe that these observations are consistent with our mechanistic hypothesis concerning the protective effect of BCG against COVID-19. These conclusions are supported by the enrichment results produced using the 'Compare Experiment' algorithm in MetaCore from Clarivate Analytics, which looks for significant coordinated gene expression effects across all experiments to test whether the pathway is being up-or down-regulated in a manner that is unlikely to be accounted for by random chance. The top enriched pathway map, with upregulated genes in response to BCG, is 'Immune response T cell subsets: secreted signals' (Fig. 4b) .

A recent study showing that SARS-CoV-2 reshapes central cellular pathways, such as translation, splicing, carbon metabolism and nucleic acid metabolism (53), provides further support for this observation.

Naturally, bioinformatics techniques relying on gene expression, pathway over-representation and network biology have some limitations and biases: 1) results are impacted by the user-selected cut-off thresholds used to determine significant genes, which could make the results userdependent (63) ; 2) all components in the pathway are given equal weights without paying attention to the nature of the interactions between the different components (64) ; 3) there are underlying assumptions that pathways are independent of each other, contrary to the fact that pathways cross-talk and overlap (65); 4) assuming independence between genes may result in false positive predictions of highly enriched genes or gene sets (66); 5) these methods are incapable of modeling an organism's biology as a dynamic system, and cannot predict changes in the system due genetic mutations or environmental changes (67); 6) most pathway knowledge databases are built by curating experiments performed in different cell types at different time points under different conditions, so they are missing condition-and cell-specific information (64) .

In order to mitigate some of the aforementioned limitations, we used a consensus gene signature since it is more stable than other gene expression signatures, we paired overrepresentation pathway analysis with causal reasoning to b a Fig. 4 (a) A venn diagram showing overlaps between BCG genetic mimics and key hubs with SARS-CoV-2 and Corona viruses interactomes. (b) Top "pathway map" with the highest level of enrichment by genes in BCG-CGS. This map is generated using MetaCore from Clarivate Analytics. Red thermometers indicate genes overexpressed in response to BCG treatment, and the hight of the red bars is representative of the differential gene expression level (i.e., log 2 values of the fold change). The numbers under the thermometers 1-5 refer to the experiment number: 1) gene expression on day 0 in response to BCG vaccination to a BCG-naïve population on day 1; 2) gene expression on day 0 in response to BCG re-vaccination to a previously vaccinated population; 3) gene expression on day 14 in response to BCG vaccination to a BCG-naïve population; 4) gene expression on day 14 in response to BCG revaccination to a previously vaccinated population, and 5) positive connections from the connectivity map, and the red bar in the thermometer number 5 represents presence of the gene only.

predict protein activities based on the nature of interactions between upregulated or downregulated genes, and we also integrated results from several bioinformatics methods such as causal reasoning and CMap predictions to prioritize common hypotheses.

A recent publication (68) in Lancet has questioned whether BCG's effects can last for a long time. Our top enriched pathway map (Fig. 4b) indicates that BCG's effects can be longlasting if the effects were exerted on thymus-generated Th 0 cells, which can occur to a greater extent very early in life before reaching thymic involution by puberty (69) . This pathway map indicates that BCG is capable of affecting both the numbers and the types of produced innate immune cells, as well as their maturation to long-lived memory T cells (i.e., what is known as trained immunity). This is very significant in the context of BCG's protective effects from SARS-CoV-2 and other emergent novel viruses where the individual's ability to eradicate such viruses is dictated by the number and diversity of naive T cell reservoir (70, 71) . Our analysis suggest that BCG may protect individuals from novel pathogens by priming their trained immunity to fight such pathogens, including SARS-CoV-2.

Supporting evidence for this hypothesis is found in the literature (72) indicating that the protective effects of the BCG against TB, can last from 15 to 60 years after vaccination (72, 73) , with longer lasting effects when the vaccine is administered during the first year of life (74, 60) . A recent study indicated that "school-aged BCG vaccination offered moderate protection against tuberculosis for at least 20 years, which is much longer than previously thought" (60, 72) . Another 60year follow-up study, showed that BCG vaccine efficacy persisted for 50 to 60 years after a single dose of BCG (60) . Of special interest is a recent study that showed that mucosal vaccination resulted in an increased frequency of antigenspecific lung tissue-resident CD4+ T cells that provide longterm immunity (75) . These studies serve as additional evidence from the literature supporting our claim that a single dose of an 'effective' BCG vaccination to infants can have a very long duration of protection against pathogens including SARS-CoV-2.

Our findings provided systems biology support for using BCG to protect from the severe consequences of COVID-19. BCG is currently on WHO's List of Essential Medicines; it is considered one of the safest and most effective medicines Unknown †Score refer to the CMap score. It represents the level of similarity between transcriptional effects induced by BCG and each of the compounds ‡ Validation refers to the presence of any supporting evidence from the biomedical literature that the predicted BCG mimics have any antiviral activities. Antiviral means there is evidence that the compound is used as or has antiviral activity; SARS-CoV-2 means that the compound should antiviral activity against SARS-CoV-2; Corona viruses means that the compound showed antiviral activity against corona viruses other than SARS-CoV-2 § COVID-19 CT: there is evidence that the compound is being tested in clinical trials for COVID-19. There are 12 Studies found for Ruxolitinib in COVID-19 on clinicaltrials.gov. needed in a health system. There is also evidence indicating that BCG can improve the response to vaccines directed against viral infections (48, (53) (54) (55) , which may prove useful when SARS-CoV-2-specific vaccines become available. Therefore, we suggest that administering the BCG vaccine to all newborns may protect them from the infection by SARS-CoV-2 and other emerging pathogens.

Since this is an approved vaccine for TB, it can directly enter Phase III testing for the protection from COVID-19 caused fatalities. However, we caution that running these experiments during an active COVID-19 outbreak, might expose participants to aggravated immune responses if they contract COVID-19 during the study. We also advise that clinical study design takes into account several factors that are known to affect the performance of BCG vaccine, such as: the age of the participants, geographies, ethnicities, route of administration and the mycobacterium strain used in the vaccine. It is equally important to run experimental validation studies, to evaluate the effects of BCG mimics, in preventing COVID-19 or for treating urological cancers.

Our results provide systems biology support for using BCG and small-molecule BCG mimics as putative vaccine and drug candidates against emergent viruses including SARS-CoV-2. Of course, any practical actions to repurpose this vaccine as a means of protection against SARS-CoV-2, or other novel viruses, should be preceded by the successful in vitro and animal experimentation. We also caution that previous studies showed that the protective effects of BCG were found to be weaker when the vaccine was given after the first year of life and particularly after puberty (68).

Research at Al-Zaytoonah University of Jordan grant 2020-2019/17/03. AT acknowledges partial support from NIH grant OT2TR003441. We thank Clarivate Analytics for providing access to MetaCore, a specialized pathway and functional genomics analysis product. Reference to commercial products or services does not constitute their endorsement.

R. Hajjo generated the idea, designed the workflow, generated content, performed data analysis and wrote the manuscript. A. Tropsha was engaged in the critical discussion of the study design, recommended some studies, and edited the manuscript. 

Developing covid-19 vaccines at pandemic speed

From COVID-19 research to vaccine application: Why might it take 17 months not 17 years and what are the wider lessons? Health Research Policy and Systems

Vaccination against coronaviruses in domestic animals [Internet]. Vaccine

Deployment of convalescent plasma for the prevention and treatment of COVID-19

The history and advances in cancer immunotherapy: understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications

Springer Nature

Novel cancer immunotherapy agents with survival benefit: Recent successes and next steps

Advances in cancer immunotherapy 2019 -Latest trends

Immuno-oncology combinations: raising the tail of the survival curve

Harnessing innate immunity in cancer therapy

Targeting innate immunity to enhance the efficacy of radiation therapy

Targeting innate immunity for antiviral therapy through small molecule agonists of the RLR pathway

Antiviral innate immunity pathways

Innate immunity to influenza virus: Implications for future therapy. Expert Review of Clinical Immunology

T cell-mediated immune response to respiratory coronaviruses

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing

BCG-Induced Cross-Protection and Development of Trained Immunity: Implication for Vaccine Design

The mechanism of action of BCG therapy for bladder cancer-A current perspective

Autophagy Controls BCG-Induced Trained Immunity and the Response to Intravesical BCG Therapy for Bladder Cancer

BCG vaccines 1 Report on BCG vaccine use for protection against mycobacterial infections including tuberculosis, leprosy, and other nontuberculous mycobacteria (NTM) infections Prepared by the SAGE Working Group on BCG Vaccines and WHO Secretariat. 2017. 20. Search of: BCG | COVID-19 -List Results -ClinicalTrials.gov [Internet

BCG vaccine protection from severe coronavirus disease 2019 (COVID-19)

Is global BCG vaccination-induced trained immunity relevant to the progression of SARS-CoV-2 pandemic?

Is BCG vaccination affecting the spread and severity of COVID-19? Allergy

BCG vaccination induced protection from COVID-19

Chemocentric informatics approach to drug discovery: identification and experimental validation of selective estrogen receptor modulators as ligands of 5-hydroxytryptamine-6 receptors and as potential cognition enhancers

Gene Expression and Cytokine Profile Correlate With Mycobacterial Growth in a

A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles

Nonparametric statistical methods

Protein phosphatase 4 is an essential positive regulator for Treg development, function, and protective gut immunity

Phosphatase PP4 negatively regulates type I IFN production and antiviral innate immunity by dephosphorylating and deactivating TBK1

Mediator complex interaction partners organize the transcriptional network that defines neural stem cells

Histone deacetylase function in CD4+ T cells

The effect of Ingenol-B on the suppressive capacity of elite suppressor HIV-specific CD8+ T cells

Prostratin as a new therapeutic agent targeting HIV viral reservoirs. Drug News and Perspectives. Drug News Perspect

Effects of prostratin on Cyclin TI/P-TEFb function and the gene expression profile in primary resting CD4+ T cells

Lipopolysaccharide and phorbol 12-myristate 13-acetate both impair monocyte differentiation, relating cellular function to virus susceptibility

LRR and PYD domains-containing protein 3 inflammasome is activated and inhibited by berberine via toll-like receptor 4/ myeloid differentiation primary response gene 88/nuclear factor-κB pathway, in phorbol 12-myristate 13-acetate-induced macrophage

Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro

Herpes simplex virus virucidal activity of MST-312 and epigallocatechin gallate

Antiviral (RNA) activity of selected amaryllidaceae isoquinoline constituents and synthesis of related substances

Antiviral activity of brefeldin A and verrucarin A

Why lopinavir and hydroxychloroquine do not work on COVID-19 -ScienceDaily

Effect of Systemic Inflammatory Response to SARS-CoV-2 on Lopinavir and Hydroxychloroquine Plasma Concentrations

Pharmacokinetics of Lopinavir and Ritonavir in Patients Hospitalized With Coronavirus Disease 2019 (COVID-19)

American College of Physicians

Master Regulator Analysis of the SARS-CoV-2/Human Interactome

A SARS-CoV-2-Human Protein-Protein Interaction Map Reveals Drug Targets and Potential Drug-Repurposing. bioRxiv

Vaccination and All-Cause Child Mortality From 1985 to 2011: Global Evidence From the Demographic and Health Surveys

Adjuvant effect of Bacille Calmette-Guérin on hepatitis B vaccine immunogenicity in the preterm and term newborn. Frontiers in Immunology

Repeated BCG vaccination is more effective than a single dose in preventing diabetes in non-obese diabetic (NOD) mice. undefined

Non-specific Effects of vaccines illustrated through the BCG example: from observations to demonstrations

Benefits of BCG-induced metabolic switch from oxidative phosphorylation to aerobic glycolysis in autoimmune and nervous system diseases

BCG-induced protection: Effects on innate immune memory

BCG vaccination enhances the immunogenicity of subsequent influenza vaccination in healthy volunteers: A randomized, placebo-controlled pilot study

Non-specific effect of Bacille Calmette-Guérin vaccine on the immune response to routine immunisations

Influence of mycobacterium bovis bacillus calmette-guérin on antibody and cytokine responses to human neonatal vaccination

The American Association of Immunologists

Report on BCG vaccine use for protection against mycobacterial infections including tuberculosis, leprosy, and other nontuberculous mycobacteria (NTM) infections Prepared by the SAGE Working Group on BCG Vaccines and WHO Secretariat

Nonspecific effects of BCG vaccine on viral infections

BCG Vaccination Protects against Experimental Viral Infection in Humans through the Induction of Cytokines Associated with Trained Immunity. Cell Host and Microbe

Bacille Calmette-Guérin induces NOD2-dependent nonspecific protection from reinfection via epigenetic reprogramming of monocytes

Long-term Efficacy of BCG Vaccine in American Indians and Alaska Natives: A 60-Year Follow-up Study

Selective and cross-reactive SARS-CoV-2 T cell epitopes in unexposed humans

Robust T cell immunity in convalescent individuals with asymptomatic or mild COVID-19

Using the gene ontology for microarray data mining: A comparison of methods and application to age effects in human prefrontal cortex

Ten Years of pathway analysis: current approaches and outstanding challenges

Network biology: understanding the cell's functional organization. nature reviews | genetics

Gene set analysis: challenges, opportunities, and future research

COVID-19 vaccination clinical trials should consider multiple doses of BCG

Thymus size and agerelated thymic involution: Early programming, sexual dimorphism, progenitors and stroma

Thymic and postthymic regulation of Naïve CD4+ T-cell lineage fates in humans and mice models

Mechanisms shaping the naïve T cell repertoire in the elderly -Thymic involution or peripheral homeostatic proliferation? Experimental Gerontology

The duration of protection of school-aged BCG vaccination in England: a population-based case-control study

Duration of BCG protection against tuberculosis and change in effectiveness with time since vaccination in Norway: a retrospective population-based cohort study. The lancet infectious diseases

Delaying BCG vaccination from birth to 10 weeks of age may result in an enhanced memory CD4 T cell response

Enhanced protection conferred by mucosal BCG vaccination associates with presence of antigenspecific lung tissue-resident PD-1 + KLRG1 − CD4 + T cells

Gene expression and cytokine profile correlate with mycobacterial growth in a human BCG challenge model. The Journal of Infectious Diseases

Cytoscape: a software environment for integrated models of biomolecular interaction networks

STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

MINT: the molecular INTeraction database

Human protein reference database as a discovery resource for proteomics

BIND: the Biomolecular Interaction Network Database

DIP: the database of interacting proteins

The BioGRID interaction database: 2008 update

KEGG for linking genomes to life and the environment

EcoCyc: a comprehensive database resource for Escherichia coli

The NCI-Nature Pathway Interaction Database: A cell signaling resource

The Gene Ontology (GO) database and informatics resource

Functional analysis of OMICs data and small molecule compounds in an integrated "knowledge-based" platform

The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease. NY).: Science

Causal reasoning on biological networks: interpreting transcriptional changes

A computational model to define the molecular causes of type 2 diabetes mellitus

New pharmacological strategies to fight enveloped viruses

A screen of FDA-approved drugs for inhibitors of Zika virus infection

Antiviral activity of brefeldin A and verrucarin A

Middle East respiratory syndrome and severe acute respiratory syndrome: current therapeutic options and potential targets for novel therapies

Novel activities of safe-in-human broadspectrum antiviral agents

Brefeldin A and Cytochalasin B reduce dengue virus replication in cell cultures but do not protect mice against viral challenge

Emetine inhibits Zika and Ebola virus infections through two molecular mechanisms: inhibiting viral replication and decreasing viral entry

The natural compound homoharringtonine presents broad antiviral activity in vitro and in vivo

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations