key: cord-0887988-l5wtb9lo authors: Patra, Sinjini; Saxena, Shivam; Sahu, Nilanjan; Pradhan, Biswaranjan; Roychowdhury, Anasuya title: Systematic Network and Meta-analysis on the Antiviral Mechanisms of Probiotics: A Preventive and Treatment Strategy to Mitigate SARS-CoV-2 Infection date: 2021-02-03 journal: Probiotics Antimicrob Proteins DOI: 10.1007/s12602-021-09748-w sha: c7dfb6b243d6573b08da6e43946fe34bde50ccb8 doc_id: 887988 cord_uid: l5wtb9lo With the alarming rise of infected cases and deaths, COVID-19 is a pandemic, affecting 220 countries worldwide. Until now, no specific treatment is available against SARS-CoV-2. The causal virus SARS-CoV-2 primarily infects lung cells, leading to respiratory illness ranging in severity from the common cold to deadly pneumonia. This, with comorbidities, worsens the clinical outcome, particularly for immunosuppressed individuals with COVID-19. Interestingly, the commensal gut microbiota has been shown to improve lung infections by modulating the immune system. Therefore, fine-tuning of the gut microbiome with probiotics could be an alternative strategy for boosting immunity and treating COVID-19. Here, we present a systematic biological network and meta-analysis to provide a rationale for the implementation of probiotics in preventing and/or treating COVID-19. We have identified 90 training genes from the literature analysis (according to PRISMA guidelines) and generated an association network concerning the candidate genes linked with COVID-19 and probiotic treatment. The functional modules and pathway enrichment analysis of the association network clearly show that the application of probiotics could have therapeutic effects on ACE2-mediated virus entry, activation of the systemic immune response, nlrp3-mediated immunomodulatory pathways, immune cell migration resulting in lung tissue damage and cardiovascular difficulties, and altered glucose/lipid metabolic pathways in the disease prognosis. We also demonstrate the potential mechanistic domains as molecular targets for probiotic applications to combat the viral infection. Our study, therefore, offers probiotics-mediated novel preventive and therapeutic strategies for COVID-19 warfare. The coronavirus disease-2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has been declared as a pandemic by the World Health Organization (WHO). To date, more than 67 million confirmed cases of COVID-19 have been reported in 220 countries, with more than 1.5 million confirmed deaths (as of 4:22 pm CET, 9 December 2020, https :// covid 19.who.int/). The death toll has climbed up since no specific treatment for COVID-19 is currently available. The global emergency of the pandemic therefore urgently demands the investigation of the novel therapeutic strategy effective for SARS-CoV-2. SARS-CoV-2 is a positive-sense RNA virus and belongs to pathogenic beta-coronavirus. The other two betacoronavirus that caused earlier major outbreaks of lethal pneumonia are SARS-CoV and MERS-CoV. SARS-CoV-2 shares almost 80% of the genome with SARS-CoV. However, SARS-CoV-2 is more distant from MERS-CoV, with only 50% identity [1] . Epidemiology indicates SARS-CoV-2 is more infectious but less fatal (2-5%) than SARS-CoV (9.4%) and MERS-CoV (34.5%) [2] . Nevertheless, similar to them, SARS-CoV-2 primarily infects alveolar epithelial cells of the lung, leading to a respiratory syndrome with a variable degree of severity, ranging from upper respiratory tract infection (URTI) to severe interstitial pneumonia and acute respiratory distress syndrome (ARDS) [3] . Older people and immunocompromised individuals with existing medical conditions like diabetes, hypertension, and cardiovascular complications are critically affected by the disease [4] . Interestingly, respiratory infections, sepsis, and ARDS are associated with a change in gut microbiota composition, indicating their possible role in pulmonary health [5] . It is also shown that the macrophage response to the respiratory viruses during the viral infection depends on the composition of gut microbes [6] . Therefore, the gut and lung are interconnected organs and influence their homeostasis through immunological communication [7] . Strikingly, similar cross-talk among gut and lungs appears in COVID-19 instances as well [7] . A clear decrease in cell number of two common gut bacteria, Lactobacillus and Bifidobacterium spp., has also been found in COVID-19 patients [8] . Next, it is also reported that bronchoalveolar lavage fluid samples of COVID-19 and community-acquired-pneumonia patients are dominated by bacteria that are common in the oral cavity and upper respiratory tract (URT) [7] . Hence, the crosstalk between the gut and the lung through the gut-lung axis may impart a significant role in SARS-CoV-2 infection [5] . Therefore, manipulation of the intestinal microbiota by application of probiotics could be a potential therapeutic strategy or adjuvant therapeutic option for maintaining health and preventing and/or treating the disease. Consumption of probiotics has also been shown to improve the URTI in numerous human clinical trials [9] (Table 1) . Since the specific strains of probiotics, applied with appropriate dosage and mode of administration, can efficiently treat respiratory complications, we hypothesize that the application of probiotics could be effective against SARS-CoV-2 infection. Probiotics build up a protective mucus barrier, impart healthy immune response to all age group individuals minimizing the probability of infectious disease occurrence [10] . Probiotics also maintain the balance of the gut microbial community and protect the respiratory system by preventing secondary bacterial infections [7] . Recently, Tiwari et al. (2020) have reviewed the beneficial strains of lactic acid bacteria (LAB) and their potential antiviral properties [11] . Metabolites like shortchain fatty acids (SCFAs), lactic acid, hydrogen peroxide, and antimicrobial peptides (AMPs) like bacteriocins produced by the LAB have been suggested to be beneficial in restricting the viral entry and/or reducing the viral load [11] . Additionally, such AMPs are also important in the colonization of probiotics in the GI tract, cell signaling, and host defense [11] . The activities of bacteriocins in reinforcing the immunity with their immunomodulatory effects further suggest the prophylactic use of probiotics against viral infections, including COVID-19 [11] . Network analysis offers an effective approach to identify molecular mechanisms and connections between genes and their pathways from dynamic networks [12] . Therefore, in the present network and meta-analysis study, a biological association network is generated concerning training genes of SARS-CoV-2 infection and probiotic treatment. The meta-analysis shows the rational justification for the implementation of the probiotics approach in respiratory infections. The analysis of functional modules, pathway enrichment, and topological network parameters reveal that probiotics could have tremendous therapeutic potential during the pathophysiological events of COVID-19. The comprehensive network analysis study indicates that the application of probiotics at the significant domains of the infection mechanism could be beneficial in the prevention and treatment of SARS-CoV-2 infection. Despite the genomic similarity with the SARS-CoV, highthroughput data and an unambiguous model for SARS-CoV-2 immunopathology are unavailable. The literature search for the study was performed according to the PRISMA guideline (Fig. 1) . The rationale for probiotics as a cure for the global pandemic was justified by the meta-analysis study (Fig. 2) . Only a few clinical trials concerning the direct application of other probiotics like yeast (Saccharomyces boulardii), Bacilli spp. (Bacillus coagulans), and Enterococcus spp. (Enterococcus faecalis) on URTI patients are available. Therefore, the study was restricted to the most common probiotic strains Lactobacillus spp. and Bifidobacterium spp. to increase the robustness and accuracy of the network metaanalysis. Further, a network-analysis study was considered to evaluate the probable bioprotective mechanisms of probiotics against SARS-CoV-2 infection. The network module analysis revealed functionally related genes and provided key domains for the potential mechanism of probiotics against SARS-CoV-2 infection. The entire methodology is summarized in Fig. S1 . The meta-analysis was conducted to analyze similar studies that are already available in the literature. In our study concerning probiotic treatment on URTI patients, Q-statistics is used to assess whether all effect sizes in the sample of single studies are homogeneous (belong to the same population). The I 2 index indicates the degree of heterogeneity. Further, the heterogeneity on the pooled estimates of the individual outcomes of the meta-analysis was assessed by the forest plot. The effects of probiotics depend on specific strains, age-group, clinical dosages, and the mode of administration. Clinical studies with randomized, double-blind placebo-controlled human trials were screened. The overall effects of specific probiotic strains were measured as mean difference (MD) or standardized mean difference (SMD), effects size, 95% confidence interval (CI), and weight percentage for studies by using a random-effects model. The mean change, i.e., the standard deviation in the number of healthy individuals and patients, was used to calculate the MD between the probiotics (or synbiotics) and the placebo groups. The statistical analyses were performed using MetaPython. The p value of Q statistic < 0.05 was defined as an indicator of heterogeneity, and data were considered heterogeneous for I 2 values higher than 40%. Begg's funnel plot was performed to examine the publication bias. All the reported p values were two-sided, and p values ≤ 0.05 were regarded as statistically significant for all included studies. The literature study was conducted following the PRISMA guidelines ( Fig. 1 ). The comparatively unexplored functional genes (training genes) relevant to the viral infection were detected from the structured databases by data-mining strategy. The primary search strategy included manual analysis of the published literature indexed in the electronic search database PubMed archive (http://www.ncbi.nlm. nih.gov/pubme d/) using reference management software EndNote X7 (Bld 7072) (endnote.com). The search terms or MeSH terms used in EndNote included "probiotics" AND "coronavirus." In the scarcity of the available studies regarding SARS-CoV-2 immunopathology and antiviral probiotic application against them, the literature search was continued using the NCBI PubMed database. The search terms or keywords used were "probiotics" ∩ "COVID-19", "probiotics" ∩ "severe acute respiratory syndrome coronavirus 2", "probiotics" ∩ "SARS-CoV-2," and "probiotics" ∩ "coronavirus." The search was limited to the studies conducted on humans, and available literature from April 18 to June 8, 2020 was considered for the present study. Therefore, data-mining considered the immunopathology of SARS-CoV-2 infection and the application of probiotics to mitigate the disease severity. This approach allowed us to pull out the training genes associated with novel coronavirus infection and probiotic treatment. The top-ranked genes were thus selected and used for text-mining. An efficient and integrated combination of data-mining and text-mining approaches was used to retrieve documents eliminating biases toward vigorously investigated disease phenotypes. The unstructured biomedical texts were regained with simultaneous construction of interaction networks by text-mining strategy. Thus, the candidate genes associated with SARS-CoV-2 infection and antiviral probiotic treatment were constructed in an association network by screening the training genes from literature based on the text-mining approach. An open-source bioinformatics tool, Cytoscape, and the metasearch plug-in Agilent Literature Search was used to visualize the molecular interaction networks, integrating the gene expression profiles of the respective training genes. Agilent Literature Search software fetched documents based on the entered query using multiple text-based search engines, parsed into sentences. The software generated the association network based on a lexicon set definition which defines the gene names of the parsed sentences as "concepts" and the interaction terms of interests as "verbs." The network was generated by extracting an association for every parsed sentence containing at least two "concepts" and one "verb" or known interaction term. The gene names and association terms represented nodes and edges respectively of the textmining-based network, generated through Cytoscape. In this work, Cytoscape 3.7.1 App and Agilent Literature Search 3.1.1 beta (LitSearch version 2.69) were used to analyze the existing data source from the published literature available in the PubMed database. The training genes were used as "search terms"; Max Engine Matches was set at 10; the "use context" and the "concept lexicon restrict search" options were selected as search controls; "Homo sapiens" was used as "extraction control". Any disease pathogenesis manifests the complex interaction of biological events modulated by the disease-related association network. Identification of tightly interconnected nodes or genes from a very densely connected network was useful for understanding significant biological events in terms of modules and their interconnections. The association network generated with SARS-CoV-2 pathogenesis candidate genes was subdivided into modules by MCODE plug-in of Cytoscape according to local neighborhood density. The extracted modules were graphically displayed as an isolated, more considerable dense region with functionally similar genes. The network view of the clusters enabled us to understand the local topology and functional features concerning the whole network consisting of several other candidate genes and edges. Network modules with MCODE score of more than three and a minimum of four nodes were considered as significant and carried forward for further analysis of functional annotation. The functional annotations and the pathways associated with candidate genes of the selected MCODE clusters were determined by GeneCards and NCBI. The network topology parameter details of the association network were obtained by analyzing the association network with the NetworkAnalyzer tool in Cytoscape. The node degree and clustering coefficient of the individual node genes together represented the number of connections associated with a node and the degree of involvement of a node in the participating clusters. This tool showed the degree of participation of individual genes in the whole network and participating clusters. The functional annotations of the highest MCODE score with clustered candidate genes were determined by the integrative human gene database GeneCards (version 4.14). It provided gene-centric data of the annotated and predicted human genes with the functions and pathways associated with the candidate genes. To identify enriched biological processes that are affected by COVID-19 immunopathology, a gene ontology (GO) functional enrichment analysis was performed using the Biological Networks Gene Ontology tool (BiNGO, version 3.0.3) in Cytoscape 2.8.0 (http://www.cytos cape.org/), with a threshold of p < 0.001. BiNGO is a tool to determine statistically overrepresented GO terms for a set of genes associated with any biological processes as recorded in the Gene Ontology database. The BiNGO analysis was done with the MCODE-derived gene clusters of the association network. The statistical test was set to "hypergeometric test," multiple testing correction as "Benjamini & Hochberg False Discovery Rate (FDR) correction," the significance level was set to 0.05, the categories to be visualized were "overrepresented categories after correction," reference set was "use the whole annotation as reference set," ontology file was "GO_Biological_Process," and finally organism/ annotation was selected as "Homo sapiens." URTI is one of the principal symptoms of SARS-CoV-2 infection. Therefore, in all the meta-analysis studies, the duration and severity of URTI were compared between individual probiotic or probiotic formulation (or synbiotics) treatment and placebo-controlled groups, including children and adult individuals . The forest plot showed the effect of probiotic treatment on URTI (Fig. 2a) . As indicated by an I 2 value of 31.47% and Q (chi-square) statistics of 42.32 (p = 0.05), there was no statistical heterogeneity. The overall effect size of the study was 2.75, with p = 0.006. The effect of probiotic treatment on URTI was plotted according to specific strains of probiotics. The effects size (ES), confidence level (95% CI), and the weight percentage of clinical trials concerning the specific strains were also depicted. Thus, the forest plot of the meta-analysis study showed statistically significant effects of probiotic treatment on URTI. The funnel plot showed minimal evidence of publication bias among the selected studies associated with randomized, double-blind, placebo-controlled human clinical trials of probiotic application on URTI (Fig. 2b) . The funnel plot with strain-specific probiotic treatment on URTI patients was symmetrical, which indicated that no significant asymmetry was detected in any of the analyses. The data-mining exploratory technique found out about 36 data, including clinical trials, in vitro experiments, bioinformatics analysis, reviews, minireviews, and editorial letters associated with COVID-19 pathophysiology and probiotic treatment from the literature [1, [41] [42] [43] [44] [45] [46] [47] . The use of specific MeSH terms in EndNote and keywords used in the PubMed database excluded non-specific search results. This confirmed the accuracy of the mining study performed with the reports exclusively on SARS-CoV-2 infection and the application of probiotics. Manual analysis of the literature shortlisted 90 training genes of interest associated with COVID-19 pathogenesis (Table 2 ). These genes were used for the generation of the association network and further analysis. All the 90 shortlisted training genes were used for the development of a text-mining-based association network. The primary association network obtained through Cytoscape with training genes showed 453 nodes representing the candidate genes and 1273 edges. [44] crp, hmgcoa, app, scfa, gpr43, hdac Akour [45] g-csf, ip10/cxcl10, mcp-1, thp-1, il22, il8, tnfα, mpo, inos, cox-2, nf-κb, il6 Renzo et al. [46] ace2, il10, il17, mtor Shi et al. [47] hla, has2 indicated the topology of the central association network, highlighting the listed training genes involved in the SARS-CoV-2 infection. Table 3 shows the network parameter statistics. According to the network topology, "the number of connected components" represented the network connectivity, which reflected a pairwise connection among all the nodes. Hence, a lower number of connected components of 27, obtained by the network statistics, suggested robust connectivity of the association network ( Table 3 ). The clustering coefficient 0.633 designated the average clustering coefficient or degree involvement of respective nodes in the participating cluster of the network, with 5.620 average number of neighbors. The tendency to form high degree nodes or hub nodes was described by high "network heterogeneity" of 1.137 and the absence of any isolated nodes. The presence of hub nodes indicated a real functional network with multiple biologically relevant pathways associated with the novel coronavirus infection. The text-mining approach in conjunction with Cytoscape has been successfully adopted to illustrate actin dynamics during the post-ejaculatory life of spermatozoa [48] , to demonstrate the synergistic mechanisms of therapeutic herbs for rheumatic arthritis [49] , and to find out the effects of psychological stress on innate immunity and metabolism in human [50] . MCODE tool of Cytoscape derived 38 densely bridged modules from the first network (Fig. S3 ) by an algorithm called "top overlap." It made groups of "genes of interest" according to their correlation coefficient and highest-scoring edge among all other candidate genes of the network. The algorithm also minimized the chances of false-positive results, increasing the robustness and accuracy of the analysis [51] . The individual clusters were further selected based on their scores for functional studies with the removal of lower confidence data points. The details of the modules are given in table S1. Among the 38 derived modules, 11 resultant clusters with MCODE score ≥ 3, nodes ≥ 4, edges ≥ 6 were chosen for functional annotation (Fig. 3 ). Module 1 with MCODE score 16, 16 nodes, and 120 edges was the highest-scoring module, whereas module 2 with MCODE score 14, 14 nodes, and 91 edges was the second-highest, and module 3 with score 8.457, 36 nodes, and 147 edges was the third-highest. The network matrices like the node degree, the clustering coefficient of the primary network, and MCODE-derived modules provided the topological parameters to understand the importance of nodes for pathway analysis. The highestscoring module 1 (MCODE score 16) consisted of 16 nodes (src, limk1, rps6ka3, aak1, mapk1, mknk2, map3k1, gak, map2k1, fgfr1, mapkapk5, map3k7, bmp2k , zak, gsk3, yes1) (Fig. 3a) A higher-scoring MCODE module with a high number of nodes and edges implied a possible multi-functional role. Such modules could be assigned in the essential regulatory functions of the disease pathogenesis and thereby could be the target for the probiotic treatment. The genomic information of the nodes showed that MCODE cluster 1 candidate genes src, limk1, rps6ka3, aak1, mapk1, mknk2, map3k1, gak, map2k1, fgfr1, mapkapk5, map3k7, bmp2k, zak, gsk3, yes1 were involved in the receptormediated endocytosis and phagocytosis, stress-mediated cellular metabolic pathways (Fig. 3a) . Therefore, probiotic supplementation could heal the disrupted intestinal barrier and prevent virus entry. Table 3 Network statistics of the association network of the probiotics-COVID-19 axis The table contains details of the primary association network obtained by text-mining results using an Agilent Literature search (ALS) plugin. "Network clustering coefficient" is the average of the clustering coefficients for all nodes in the network. The "average number of neighbors" indicates the average connectivity of a node in the network (Fig. 3b) response to stress stimuli, induction of the NF-κb-mediated inflammatory pathways, and subsequent apoptosis as well as the production of cytokines. Several probiotic strains are reported to balance the immune activations and inflammatory cytokines, which could have immense implications in the SARS-CoV-2-induced immunological complications. Although MCODE cluster 3 was the third in rank based on score, it contained the highest number of nodes and edges that might indicate multiple cellular pathways (Fig. 3c) . The GeneCards-derived functional profile of nodes expectedly showed the interconnections of four cellular events upon the viral infection: (a) activation of principal antiviral interferons (IFNs) and innate immune cells (stat3, stat5, ifna1, soat1, cd68, il9, itgam, ptprc, cd14); (b) activation of proinflammatory cytokines and their downstream signaling through JAK-STAT pathway (ccl2, chi3l1, csf2, il5, il25, klk15, lcn2, il18) , (c) activation of proinflammatory cytokines and chemokines resulting in the immune cell (neutrophil, monocyte) migration and infiltration in lung tissues, and (d) synthesis of lung-damaging hyaluronan by has2 or hyaluronan synthase-2 (cxcl5, cxcl12, tnc, has2, cxcl6, mmrn1, mmp13, smox, rela) , a proinflammatory cytokine-induced growth factor, responsible for vascular permeability and pulmonary dysfunction (vegfa, eng, plau, pgf, prl, angpt2, igfbp1, hgf, erbb2, fgf2) . Although limited studies are available on the ability of probiotics to degrade glycosaminoglycan molecules, probiotics could have an impact on the degradation of limk1, rps6ka3, aak1, mapk1, mknk2, map3k1, gak, map2k1, fgfr1, mapkapk5, map3k7, bmp2k, zak, gsk3, yes1) , and edges 120 (a). MCODE-derived cluster 2: score 14, nodes 14 (tnf, il6, map1lc3b, ros1, cd38, bax, fas, sqstm1, il10, fcgr3a, cybb, il32, bcl2, icam1) , and edges 91 (b). MCODE-derived cluster 3: score 8.457, nodes 36 (cxcl5, cxcl12, tnc, has2, cxcl6, mmrn1, mmp13, smox, rela, ccl2, chi3l1, csf2, il5, il25, klk15, lcn2, il18, vegfa, eng, plau, pgf, prl, angpt2, igfbp1, hgf, erbb2, fgf2, stat3, stat5, ifna1, soat1, cd68, il9 , itgam, ptprc, cd14), and edges 148 (c). MCODE-derived cluster 5: score 6.769, nodes 14 (tlr2, tlr4, nlrp3, myd88, hmgb1, mapk14, creb1, rps6ka5, jun, ephb2, mapk8, atf2, fos, maa) , and edges 44 (d). MCODE-derived cluster 6: score 5, nodes 5 (dgat2, adipor2, dgat1, pgc, adipoq), and edges 10 (e). MCODE-derived cluster 9: score 4.818, nodes 23 (il27ra, il27, cxcl10, il17d, il6st, il6r, jak2, ptpn18, mapk3, epo, akt1, egfr, reg1a, ptpn1, cat, frap1, malat1, acan, eif-4ebp1, smad4, smad6, smad2, smad7) , and edges 53 (f). MCODEderived cluster 10: score 4.167 and 4, nodes 13 (ace2, mas1, ace, ang, agtr1, il22, il2, gli2, cd4, cd40lg, fus, th1l, pdgfb) (g). MCODEderived cluster 11: score 4, nodes 5 (rab18, rab13, mtg1, rce1), edges 25 (h). MCODE-derived cluster 13: score of 4, nodes 4 (hif1a, dlk1, adam17, epas1), and edges 6 (i). MCODE-derived cluster 35: score of 3, nodes 7 (pparg, twist1, lpa, slc12a, mbtps1, srebf1, mbtps2), and edges 9 (j). MCODE-derived cluster 38: score of 2.857, nodes 8 (hk2, nr1i2, foxa2, pik3ca, f10, inpp5d, hspb1, mcl1), and edges 10 (k) hyaluronan substances synthesized by SARS-CoV-2-induced proinflammatory cytokines. The candidate genes of cluster 5 (tlr2, tlr4, nlrp3, myd88, hmgb1, mapk14, creb1, rps6ka5, jun, ephb2, mapk8, atf2 , fos, maa) represented the genes responsible for TLRmediated innate immune responses and subsequent activation of inflammasome complex, with stress-induced MAPK signaling pathway (Fig. 3d) . The virus entry might activate the host's innate immune system through TLR-myd88-mediated pathways and eventually activate the cellular inflammasome complex. The functional annotations of cluster 6 candidate genes (dgat2, adipor2, dgat1, pgc, adipoq) indicated glucose and lipid metabolic pathways and innate immune responses (Fig. 3e) . Host defense mechanisms by the proliferation of Th1 cell suppressing Th2 and Th17 cells, isotype switching and anti-inflammatory cytokine production, and induction of growth regulatory pathways were represented by MCODE cluster 9 candidate genes (il27ra, il27, cxcl10, il17d, il6st, il6r, jak2, ptpn18, mapk3, epo, akt1, egfr, reg1a, ptpn1, cat, frap1, malat1, acan, eif4ebp1, smad4, smad6, smad2, smad7) (Fig. 3f) . The genomic functional annotations of cluster 10 (ace2, mas1, ace, ang, agtr1, il22, il2, gli2, cd4, cd40lg, fus, th1l, pdgfb) and cluster 11 (rab18, rab13, mtg1, rce1) candidate genes showed the reported ACE2-ADAM17mediated entry of SARS-CoV-2 in the host cell with the resulting dysregulation of the renin-angiotensin system (RAS) involving ace2, ace, mas1, ang genes ( Fig. 3g and h) . Moreover, cluster 13 candidate genes (hif1a, dlk1, adam17, epas1) indicated the receptor-mediated endocytosis of the virus entry with the consequent activation of oxidative stressresponsive pathways mediated by hif1a (Fig. 3i) . Probiotic supplementation could alter the adverse conditions of the infection by regulating host metabolic pathways through the secretion of SCFAs. Additionally, probiotics could balance the inflammatory cytokines and chemokine levels in the body by regulating oxidative stress, altering the SARS-CoV-2 infection condition. The effects of probiotics in controlling blood pressure may significantly alter ace2mediated dysregulation of RAS. Finally, the candidate genes of MCODE cluster 35 (pparg, twist1, lpa, slc12a, mbtps1, srebf1, mbtps2) and cluster 38 (hk2, nr1i2, foxa2, pik3ca, f10, inpp5d, hspb1, mcl1 ) designated the involvement of glucose and lipid metabolic pathways, adipocyte differentiation, VEGF signaling pathway, and induction of atherosclerosis ( Fig. 3j and k) . Although the MCODE score of cluster 38 was <3, the module was selected for analysis due to the 8 nodes and 10 edges that might reflect significant function and pathways in the disease. Treatment with probiotic formulations might reverse any metabolic abnormality upon SARS-CoV-2 infections by releasing SCFAs. The molecular events and pathways indicated by MCODE clusters were supported by the published scientific shreds of evidence available for SARS-CoV-2 infection and pathogenesis, which elucidated the promising yet unexplored functional pathways of these genes. The enrichment or overrepresented GO terms associated with the MCODE cluster candidate genes, designated by p values in BiNGO enrichment analysis, are detailed in table S2. The analysis showed that cluster 1, cluster 2, cluster 3, and cluster 5 involved the pathways related to viral infection, pathogenassociated molecular patterns (PAMPs)-mediated immune response to extracellular signals, stress, leukocyte migration, activation of cellular oxidative stress and inflammatory processes through phosphate metabolic processes, regulation of MAPKKK cascade, and jun kinase activity. Therefore, these physiological events associated with SARS-CoV-2 infections could be modulated by probiotic supplementation. Probiotics could heal the intestinal barriers and minimize viral entry, as well as balance the immune responses and inflammatory conditions. Cluster 6, cluster 9, cluster 10, cluster 11, and cluster 13 were involved in the regulation of innate immune cell differentiation, regulation of glucose and lipid transport, cholesterol efflux, glycerol, glycerolipid, acylglycerol, triglyceride, TNF, and cytokine signaling by JAK-STAT pathway in addition to regulation of renin-angiotensin-mediated blood pressure. Probiotic supplementation could play an immense role in immune cell development and differentiation, regulation of inflammatory processes, metabolic activity, and the regulation of blood pressure through RAS. Cluster 35 and cluster 38 genes were involved in lipid, cholesterol, steroid, glucose, and monosaccharide metabolism, xenobiotic drug transport, regulation of monocyte and neutrophil differentiation, and response to cytokine IL6. This enrichment analysis of the GO terms of candidate genes enlightened us on the major biological processes associated with SARS-CoV-2 infection pathogenesis. It provided the cellular domains for antiviral probiotic mechanisms to defeat the viral infection. Several research groups have used similar BiNGO enrichment analysis to understand breast cancer susceptibility or to detect key pathways responsible for gastric cancer or analyze the protein-protein network [12] . The successful application of probiotics in respiratory infections of human clinical trials encourages us to study the probiotics application to alleviate SARS-CoV-2 infection. The meta-analysis is an epidemiological study design that quantitatively examines the available outcomes and estimates the effect of treatment much precisely than a single study contributing to the pooled analysis. In this study, by using a systematic network and meta-analysis, we provide a promising probiotic mechanism that can reinforce the immunity and mitigate the SARS-CoV-2 infection. ACE2 receptor is a master regulator of RAS and blood pressure. However, SARS-CoV-2 enters the host enterocytes and lung tissues through ADAM17-mediated proteolytic cleavage of the ACE2 receptor [53] . The disruption of the protective "gut-lung axis" increases the propensity of the infection in the underlying ACE2 containing epithelial cells. SARS-CoV-2 infection may downregulate cellular ACE2 expression with the resulting accumulation of angiotensin II (Ang II). Therefore, infected lung alveolar and small intestinal epithelial cells, heart, kidney, vascular endothelial, and smooth muscle cells (highly expressing ACE2) cause vasoconstriction, tissue inflammation, and oxidative stress, which worsen the health conditions of COVID-19 patients [53] . Our functional annotations and enrichment analysis of candidate genes of MCODE cluster 10 (ace2, mas1, ace, ang, agtr1, il22, il2, gli2, cd4, cd40lg, fus, th1l, pdgfb) and MCODE cluster 13 (hif1a, dlk1, adam17, epas1) ( Fig. 3g and i; Table S3 ) indicate that probiotics could have potential roles in limiting viral entry through ADAM17-ACE2-mediated receptor endocytosis. Further enrichment analysis shows that the application of probiotics could mitigate the adverse effects of dysregulated RAS system, hif1a-mediated oxidative stress [54] , activated immune cells (monocytes and NK cells), and elevated proinflammatory cytokines. Probiotics are known to heal the damaged epithelial barrier and thereby protect the underlying ACE2 receptorexpressing cells [55] . Bacteriocin-producing probiotic strains can colonize in the GI tract better than bacteriocin nonproducers. Therefore, bacteriocins also play a crucial role in microbial colonization [11] . Probiotics have been reported to reduce nitric oxide (NO) production [56] , hypertension [57] and oxidative stress [58] , and secrete ACE-inhibitory peptides [59] , SCFAs (acetate, propionate, and butyrate) to induce anti-inflammation and control the blood pressure [60] . Recombinant probiotic L. paracasei in conjugation with ACE2 has been applied for the cardiovascular protective role of ACE2 in mice [61] . Additionally, cinnamycin-like lantibiotics (bacteriocins) mediate inflammatory responses by inactivating phospholipase A2. One such lantibiotic ancovenin is reported to inhibit ACE, thereby preventing the conversion of Ang I to Ang II and regulate the blood pressure [11] . Frequent interaction with the external environment in terms of food and air alters the microbial composition of intestinal microbiota as well as lungs. Several GI disorders also have manifestations of respiratory complications. The high expression of ACE2 receptor in gut and lung epithelial cells with the prevalence of respiratory and GI tract infections in COVID-19 patients strongly suggests the bi-directional cross-talk of the "gut-lung axis" [7] . After entry, the viral RNAs activate the cellular innate immune system (TLRs) [62] and inflammatory pathways (NLRP3 and NF-κB) as a protective mechanism. Activated TLRs promote first-line antiviral response through MYD88 and IRF3/7-mediated type-I IFN production [63, 64] . The NLRP3 inflammasome complex assists in the secretion of proinflammatory IL1b and IL18, which activates T-cells or macrophages to secrete IL6 and TNFα. The released proinflammatory cytokines (IL1B, IL18, IL6, TNFα) further convert naive T-cells to Th1/CTLs/CD8+ or Th17 and triggers the secretion of proinflammatory IFNγ and IL17. Therefore, the gut microbiome controls innate and adaptive immune responses against the respiratory pathogens [7] . However, the NF-κB pathway can be turned on in either of two ways: activated NLRP3 or TLR4 and stress-induced MAPK signaling (ERK1/2, JNKs, and p38/MAPK14) pathway [65, 66] . The activated NF-κB pathway contributes to proinflammatory cytokine secretion and apoptosis in enterocytes and lung tissues. The resulting tissue injury activates local circulatory innate immune cells and establishes a proinflammatory feedforward loop of cytokines called cytokine storm [65] . The surge of cytokines and chemokines induces VEGF, IL8, and additional IL6 and reduces E-cadherin expression on endothelial cells ensuing vascular permeability, while it also elicits leukocyte trafficking and migration of monocytes, neutrophils, NK cells, macrophages, and dendritic cells in the lung cells. The resulting immune cells and cytokines induce hyaluronan synthesis, which participates in the pathophysiology of ARDS, the hallmark of SARS-CoV-2 infection [67] . The "gut-lung axis" facilitates the migration of immune cells from the gut to the lung through circulation and enforces immunity in health and disease. Therefore, the immunological coordination between the gut and lung affects each other's homeostasis and influences the acquired immunity in health and disease of the host. Our functional annotations and enrichment analysis of MCODE cluster 5 candidate genes (tlr2, tlr4, nlrp3, myd88, hmgb1, mapk14, creb1, rps6ka5, jun, ephb2, mapk8, atf2, fos, maa) indicate that probiotics could be associated with TLR-mediated innate immune response during SARS-CoV-2 infection. TLR-3, 7/8, reported to recognize the SARS-CoV-2 RNA, was also present in the association network but not in the MCODE cluster ( Fig. 3d; Table S4 ). The presence of nlrp3 in cluster 5 hints the activation of inflammatory and antiviral IFNs which is further supported by cluster 2 (tnf, il6, map1lc3b, ros1, cd38, bax, fas, sqstm1, il10, fcgr3a, cybb, il32, bcl2, icam1) and cluster 3 (cxcl5, cxcl12, tnc, has2, cxcl6, mmrn1, mmp13, smoxlus, rela, ccl2, chi3l1, csf2, il5, il25, klk15, lcn2, il18, vegfa, eng, plau, pgf, prl, angpt2, igfbp1, hgf, erbb2, fgf2, stat3, stat5, ifna1, soat1, cd68, il9, itgam, ptprc, cd14) candidate genes (Fig. 3b and c; Table S5 ). The highest node degree of IL6 and all the cytokine and chemokine genes confirm their significant participation indicating reported elevation in the immunopathology of COVID-19. Additionally, apoptotic genes like bax, bcl2, fas, and cybb may indicate TNFα, IL6-mediated programmed cell death of T-cells depicting the observed global lymphocytopenia in the COVID-19 patients [68] . However, receptors like cd68 (NK cells), cd14 (monocytes), and has2 gene may depict the rapid migration and infiltration of innate immune cells in the lung tissue with injury in response to chemoattractant chemokines. The candidate genes of cluster 9 (il27ra, il27, cxcl10, il17d, il6st, il6r, jak2, ptpn18, mapk3, epo, akt1, egfr, reg1a, ptpn1, cat, frap1, malat1, acan, eif4ebp1, smad4, smad6, smad2 , smad7) and cluster 1 (src, limk1, rps6ka3, aak1, mapk1, mknk2, map3k1, gak, map2k1, fgfr1, mapkapk5, map3k7, bmp2k, zak, gsk3, yes1) show the involvement of MAPK signaling and proinflammatory cytokine-induced growth regulatory pathways (Figs. 3f and S1; Table S6 ). Therefore, our analysis indicates that probiotics could play a protective role in the above-mentioned signaling pathways. Probiotics protect the intestinal barrier by inhibiting cytokine-induced intestinal epithelial cell apoptosis [69] . A probiotics mixture consisting of L. acidophilus, L. casei, L. reuteri, B. bifidium, and Streptococcus thermophilus is reported to induce both T-cell and B-cell hyporesponsiveness and downregulate T helper (Th) Th1, Th2, and Th17 cytokines without inducing apoptosis [70] . Further, probiotic-derived molecules (i.e., bacteriocins) are also known to trigger CD4(+) and CD8(+) T cell activation and modulate interleukin production [11] . Thus, probiotics could have a promising role against lymphocytopenia or reduction of T cell count, which is observed in COVID-19 patients. Moreover, probiotic formulations can balance pro-and antiinflammatory cytokine secretion, which is the key factor for robust immune (adaptive and innate) system. Although few reports are available on the immunomodulatory effects of bacteriocins, the activity of nisin, nisaplin, and pediocin is well documented [11] . Several clinical trials on URTI individuals have reported the decreased proinflammatory and increased anti-inflammatory cytokine upon probiotic application. Limited search has been done, which reveals the ability of probiotics to degrade glycosaminoglycan molecules [71] that may result in reduced expression of has2 and hyaluronan breakdown. Probiotics are also proven to downregulate the NF-κB signaling pathway by regulating MAPK and ERK pathways reducing systemic inflammation [72] . COVID-19 comorbidities like obesity, metabolic syndrome (e.g., IBD), are associated with an inadequate gut-microbial composition linked to disruption of the bi-directional communication network in the gut-heart axis. Such dysbiosis reduces gut-barrier integrity, accelerates intestinal inflammation, and contributes as an extragenomic factor to cardiovascular complications [73] . Dysregulation of RAS due to SARS-CoV-2 entry promotes inflammation, vasoconstriction, hypertrophy, proliferation, and fibrosis, all factors that contribute to the development and progression of cardiopulmonary diseases [74] . Chronic elevation of IL6 and cytokine storm promotes macrophages to release MCP-1, which aids atherogenesis, expression of cell adhesion molecules, and proliferation and migration of vascular smooth muscle cells resulting in cardiovascular diseases (CVDs) (e.g., coronary atherosclerosis, inflammation in the vascular system, and thrombosis) [75] . Enhanced angiotensin II also causes EGFR transactivationinduced vascular remodeling [76] . The induction of CVDs further functionally damages intestinal epithelial cells (IECs), which, together with barrier dysfunction, allows translocation of gut bacteria-derived endotoxins (lipopolysaccharide or LPS) and metabolites into the circulation. The recognition of LPS by receptors like TLR4 on cardiomyocytes induces tissue inflammation and increases cytokine secretion (TNFα). The pathophysiology of coronary atherosclerosis is also associated with aberrant production of metabolites like betaine, choline, and TMAO by microbial digestion. These metabolites and LPS affect cholesterol and lipoprotein (LDL) metabolic pathways, which amplifies inflammation and worsen SARS-CoV-2 infection [73] . The candidate genes encoding growth factors of MCODE cluster 3 (vegfa, eng, plau, pgf, prl, angpt2, igfbp1, hgf, erbb2, fgf2, stat3, stat5) along with overrepresentation of GO terms indicates the role of probiotics in the mitigation of elevated angiotensin II-induced cardiovascular complications ( Fig. 3c ; Table S5 ). The functional annotations and enrichment analysis of the candidate genes of MCODE cluster 6 (dgat2, adipor2, dgat1, pgc, adipoq), cluster 35 (pparg, twist1, lpa, slc12a, mbtps1, srebf1, mbtps2), and cluster 38 (hk2, nr1i2, foxa2, pik3ca, f10, inpp5d, hspb1, mcl1) indicates that probiotics could have potential roles in glucose and cholesterol or lipoprotein metabolism (Fig. 3e , j, k; Table S7 ). The lipidomic and cholesterol metabolic abnormalities due to SARS-CoV-2 infection are poorly documented. The limited reports indicate enhanced glucose and lipid need for viral replication and metabolism since viruses hijack the host's metabolic processes. The overrepresentation of GO terms associated with lipid and cholesterol metabolism reflects the lipiddemanding processes (viral replication, endocytosis, and exocytosis) involved in SARS-CoV-2 infection [77] . They also indicate cardiovascular complications due to gut-dysbiosis and COVID-19. Probiotic supplementation has been proved to reduce total cholesterol, LDL, and triglycerides and increase HDL count. The antihypertension capability of the probiotics makes them an affordable and adjunctive treatment option in hypertension, diabetes, and cardiovascular diseases, and other dyslipidemiaassociated health issues [78] . SCFAs released by gut bacteria and probiotics are reported to regulate blood pressure (BP). Receptors like GPR41 and Olfr78 are responsible for hypo-and hypertensive effects of SCFAs, respectively [73] . A number of LAB probiotic strains are reported to exhibit antiviral activity by producing bacteriocins, which mostly inhibit viral replication that indicates their potential role in reducing SARS-CoV-2 infection [11] . Our systematic network and meta-analysis study aid us to propose a mechanistic model of probiotic actions in the alleviation of COVID-19 (Fig. 4) . Probiotic consumption could reduce the propensity of viral entry by healing the ACE2 containing epithelial barrier. The antimicrobial peptides or bacteriocins, SCFAs, and ACE inhibitory peptides released by beneficial bacteria could Fig. 4 The mechanistic model of probiotics action against COVID-19. Probiotics induce a stronger epithelial barrier that prevents viral entry through the gut (i), probiotics modulate gut microbiota and induce the synthesis of SCFAs that regulate blood pressure and inflammation (ii). Probiotics also release ACE-inhibitory peptides that could reduce angiotensin II (Ang II) expression, thereby inhibiting viral entry into the cell (iii). Probiotics induce anti-inflammation by suppressing NF-κb signaling and reducing the levels of IL1β, IL18, NO, and TNF (iv). Bacteriocin and other anti-, as well as proinflammatory cytokines produced by the effects of probiotics modulate Th1, Th2, and Th17 cells (v). Which, in turn, help in the production of more anti-inflammatory cytokines (vi). The anti-inflammatory cytokines regulate monocytes, macrophages, dendritic cells, and neutrophils (vii) to downregulate SARS-CoV-2 infection-mediated cytokine storm (viii), resulting in decreased total cholesterol, LDL, triglycerides, VEGF, EGF, PDGF, TNF, and CRP level in the blood stream (ix). The reduced cytokine storm and inflammation exerted by probiotics cause the reduction in hyaluronan synthesis, which eventually could improve the ARDS condition in SARS-CoV-2 infection (x) balance the dysregulated RAS. Hence, the blood pressure or cardiovascular complications can be mitigated. The reduced NO production could decrease cellular oxidative stress. This can lead to the downregulation of inflammatory (NLRP3 and NF-κB) pathways and eventually, the secretion of proinflammatory cytokines or chemokines. Probiotics and bacteriocins might balance pro-and anti-inflammatory cytokine levels and increase the T-cell count in the SARS-CoV-2-infected patients. Finally, probiotics might also degrade the hyaluronan and hence could improve ARDS. Therefore, probiotics could be considered a potential preventive and alternative treatment strategy for both mild and severe stages of COVID-19. However, further research is required for elucidation. Immune response in COVID-19: addressing a pharmacological challenge by targeting pathways triggered by SARS-CoV-2 The outbreak of COVID-19: an overview A novel coronavirus from patients with pneumonia in China Repurposing antiviral protease inhibitors using extracellular vesicles for potential therapy of COVID-19 Gut microbiota and Covid-19-possible link and implications Respiratory viral infection-induced microbiome alterations and secondary bacterial pneumonia Immunological co-ordination between gut and lungs in SARS-CoV-2 infection Management of corona virus disease-19 (COVID-19): the Zhejiang experience Probiotics and paraprobiotics in viral infection: clinical application and effects on the innate and acquired immune systems Using probiotics to flatten the curve of coronavirus disease COVID-2019 pandemic Probiotics at war against viruses: what is missing from the picture? A systems biology approach to detect key pathways and interaction networks in gastric cancer on the basis of microarray analysis Probiotics and child care absence due to infections: a randomized controlled trial Probiotic effects on cold and influenza-like symptom incidence and duration in children A new chance of preventing winter diseases by the administration of synbiotic formulations Bifidobacterium animalis subsp. lactis fails to prevent common infections in hospitalized children: a randomized, double-blind, placebo-controlled study Bifidobacterium animalis subsp. lactis in prevention of common infections in healthy children attending day care centersrandomized, double blind, placebo-controlled study Alterations of intestinal flora and the effects of probiotics in children with recurrent respiratory tract infection Efficacy of a synbiotic supplementation in the prevention of common winter diseases in children: a randomized, double-blind, placebo-controlled pilot study Probiotics and vitamin C for the prevention of respiratory tract infections in children attending preschool: a randomised controlled pilot study Role of short-term use of L. acidophilus DDS-1 and B. lactis UABLA-12 in acute respiratory infections in children: a randomized controlled trial Probiotic supplements beneficially affect tryptophankynurenine metabolism and reduce the incidence of upper respiratory tract infections in trained athletes: a randomized, double-blinded, placebo-controlled trial Effect of Lactobacillus paracasei subsp. paracasei, L. casei 431 on immune response to influenza vaccination and upper respiratory tract infections in healthy adult volunteers: a randomized, double-blind, placebo-controlled, parallel-group study Yogurt supplemented with probiotics can protect the healthy elderly from respiratory infections: a randomized controlled open-label trial Preventive effect of cow's milk fermented with Lactobacillus paracasei CBA L74 on common infectious diseases in children: a multicenter randomized controlled trial Cow's milk and rice fermented with Lactobacillus paracasei CBA L74 prevent infectious diseases in children: a randomized controlled trial Randomised, double-blind and placebo-controlled study using new probiotic lactobacilli for strengthening the body immune defence against viral infections Bifidobacterium longum PL03, Lactobacillus rhamnosus KL53A, and Lactobacillus plantarum PL02 in the prevention of antibiotic-associated diarrhea in children: a randomized controlled pilot trial Oral intake of heat-killed Lactobacillus plantarum L-137 decreases the incidence of upper respiratory tract infection in healthy subjects with high levels of psychological stress Increasing work-place healthiness with the probiotic Lactobacillus reuteri: a randomised, double-blind placebo-controlled study Long-term safety and impact on infection rates of postnatal probiotic and prebiotic (synbiotic) treatment: randomized, double-blind, placebo-controlled trial Lactobacillus GG in the prevention of gastrointestinal and respiratory tract infections in children who attend day care centers: a randomized, double-blind, placebo-controlled trial Lactobacillus GG in the prevention of nosocomial gastrointestinal and respiratory tract infections Milk containing probiotic Lactobacillus rhamnosus GG and respiratory illness in children: a randomized, double-blind, placebo-controlled trial Effect of Lactobacillus rhamnosus LGG® and Bifidobacterium animalis ssp. lactis BB-12® on health-related quality of life in college students affected by upper respiratory infections Probiotics to prevent respiratory infections in nursing homes: a pilot randomized controlled trial Daily intake of fermented milk with Lactobacillus casei strain Shirota reduces the incidence and duration of upper respiratory tract infections in healthy middle-aged office workers Decreased duration of acute upper respiratory tract infections with daily intake of fermented milk: a multicenter, double-blinded, randomized comparative study in users of day care facilities for the elderly population Consumption of a fermented dairy product containing the probiotic Lactobacillus casei DN-114 001 reduces the duration of respiratory infections in the elderly in a randomised controlled trial Effects of consumption of a fermented dairy product containing the probiotic Lactobacillus casei DN-114 001 on common respiratory and gastrointestinal infections in shift workers in a randomized controlled trial Oral delivery of probiotics expressing dendritic cell-targeting peptide fused with porcine epidemic diarrhea virus COE antigen: a promising vaccine strategy against PEDV Antiviral effects of a probiotic Enterococcus faecium strain against transmissible gastroenteritis coronavirus Short report-medical nutrition therapy for critically ill patients with COVID-19 The gut-liver-lung axis. Modulation of the innate immune response and its possible role in chronic obstructive pulmonary disease Probiotics and COVID-19: is there any link? Are probiotics effective adjuvant therapeutic choice in patients with COVID-19? COVID-19 infection: the perspectives on immune responses Networks models of actin dynamics during spermatozoa postejaculatory life: a comparison among humanmade and text mining-based models A network pharmacology approach to determine the synergetic mechanisms of herb couple for treating rheumatic arthritis Effects of psychological stress on innate immunity and metabolism in humans: a systematic analysis An automated method for finding molecular complexes in large protein interaction networks Global H3. 3 dynamic deposition defines its bimodal role in cell fate transition Angiotensin-converting enzyme-2 (ACE2), SARS-CoV-2 and pathophysiology of coronavirus disease 2019 (COVID-19) Silent hypoxia: a harbinger of clinical deterioration in patients with COVID-19 Probiotics or pro-healers: the role of beneficial bacteria in tissue repair Inhibition of nitric oxide production, oxidative stress prevention, and probiotic activity of lactic acid bacteria isolated from the human vagina and fermented food Effect of probiotic fermented milk on blood pressure: a metaanalysis of randomised controlled trials Lactobacillus fermentum ME-3-an antimicrobial and antioxidative probiotic Proteolytic profiles and angiotensin-I converting enzyme and α-glucosidase inhibitory activities of selected lactic acid bacteria Current perspectives on antihypertensive probiotics Expression of human ACE2 in Lactobacillus and beneficial effects in diabetic retinopathy in mice TLR7 and TLR8 activate distinct pathways in monocytes during RNA virus infection NLRP3 inflammasome-a key player in antiviral responses Toll-like receptor 4 in acute viral infection: too much of a good thing Should we stimulate or suppress immune responses in COVID-19? cytokine and anti-cytokine interventions MAPK phosphatasesregulating the immune response Defective lung function following influenza virus is due to prolonged, reversible hyaluronan synthesis Paradigm shift for COVID-19 response: identifying high-risk individuals and treating inflammation Probiotic bacterium prevents cytokineinduced apoptosis in intestinal epithelial cells Probiotics and immune health Probiotics in human gut microbiota can degrade host glycosaminoglycans VSL# 3 probiotics regulate the intestinal epithelial barrier in vivo and in vitro via the p38 and ERK signaling pathways Gut-organ axis: a microbial outreach and networking ACE2 and microbiota: emerging targets for cardiopulmonary disease therapy Elevated serum interleukin-6 is predictive of coronary artery disease in intermediate risk overweight patients referred for coronary angiography Understanding angiotensin II type 1 receptor signaling in vascular pathophysiology The role of lipid metabolism in COVID-19 virus infection and as a drug target A mini-review of human studies on cholesterol-lowering properties of probiotics The online version contains supplementary material available at https ://doi.org/10.1007/s1260 2-021-09748 -w. Coronavirus Data Availability Data sharing not applicable to this article as no datasets were generated or analyzed during the current study. Competing Interests The authors declare that there is no conflict of interest.