key: cord-1056476-33y0zvq7 authors: Sarma, Abhijit; Phukan, Homen; Halder, Neha; Madanan, Madathiparambil Gopalakrishnan title: An in-silico approach to study the possible interactions of miRNA between human and SARS-CoV2 date: 2020-08-01 journal: Comput Biol Chem DOI: 10.1016/j.compbiolchem.2020.107352 sha: 59be5a3cd59ffe6576fe2a44ca7beaa9e4f7e88e doc_id: 1056476 cord_uid: 33y0zvq7 BACKGROUND: The progressive SARS-CoV2 outbreaks worldwide have evoked global investigation. Despite the numerousin-silico approaches, the virus-host relationship remains a serious concern. MicroRNAs are the small non-coding RNAs that help in regulating gene profiling. The current study utilized miRNA prediction tools along with the PANTHER classification system to demonstrate association and sequence similarities shared between miRNAs of SARS-CoV2 and human host. METHOD: An in-silico approach was carried out using Vmir analyzer to predict miRNAs from SARS-CoV2 viral genomes. Predicted miRNAs from SARS-CoV2 viral genomes were used for effective hybridization sequence identification along the nucleotide similarities with human miRNAs from miRbase database. Further, it was proceeded to analyze the gene ontology using miRDB with PANTHER classification. RESULT: Based on the prediction and analysis, we have identified 22 potential miRNAs from five genomes of SARS-CoV2 linked with 12 human miRNAs. Analysis of human miRNAs hsa-mir-1267, hsa-mir-1-3p, hsa-mir-5683 were found shared between all the five viral SARS-CoV2 miRNAs. Further, PANTHER classification analyzed the gene-ontology being carried by these associations showed that 44 genes were involved in biological functions that includes genes specific for signaling pathway, immune complex generation, enzyme binding with effective role in the virus-host relationship. CONCLUSION: Our analysis concludes that the genes identified in this study can be effective in analyzing the virus-host interaction. It also provides a new direction to understand viral pathogenesis with a probable new way to link, that can be used to understand and relate the miRNAs of the virus to the host conditions. The ongoing pandemic of COVID-19 which originated at Wuhan and its progressing rate of transmission into 188 countries and territories has created a havoc in the society declaring it as a global health emergency by WHO on 30 January 2020 (Rasmussen et al., 2020) . Phylogenetically, SARS-CoV2 carries a positive-sense single-stranded RNA genome; included in the family Coronaviridae, order Nidovirales, and is a β Coronavirus of 2B group which approximately 30 kb size and shares 79.5% with SARS-CoV and 96% genome similarity with Bat Coronavirus, respectively . The clinical manifestations concerned with SARS-CoV2 are fever, dry cough, low or normal peripheral white blood cell count, and low lymphocyte commonly termed as novel coronavirus-infected pneumonia (NCIP) (Qin, 2020) . Till date, the total confirmed cases were recorded as approximately 13.2 million including 5,75,540 deaths worldwide. Since then, the situation has been deteriorating in the European provinces and American regions, where the South Asian countries have also been carrying a worse COVID-19 burden. Currently more than 400 genome sequences SARS-CoV2 are available at NCBI databases and that gives a plethora of information to aid the development of the drug and vaccine. Immuno-informatics combined with molecular approaches carried out so far to study the genome characteristics include comparison within various CoV genomes, structural behavior of different proteins in SARS-CoV2 and mutational variations (30) that have been identified along the SARS-CoV2 genomes (Ahmed et al., 2020; Silipo et al., 2015; ul Qamar et al., 2020) . Further, researchers are specifically targeting the spike protein using various predicted B-cell and T-cell epitopes (Baruah and Bose, 2020; Bhattacharya et al., 2020; Feng et al., 2020; Kalita et al., 2020; Program and Sciences, 2018; Shanmugaraj et al., 2020) . These studies were not only useful to understand host-pathogen interaction but also to design valuable antiviral therapeutics. Eventually, we tried to focus on a different way to understand the mechanism of pathogenesis in SARS-CoV2 and found that the studies on miRNAs and RNAi is useful in developing as an alternative approach, where the investigators reported many miRNAs mediating gene silencing activity (Saxena and Dwivedi, 2013) . Even, these miRNAs helps in regulating diverse biological functions like development, apoptosis, tumorigenesis, proliferation, stress response, and fat metabolism (Chen et al., 2013; Mishra et al., 2020) . The miRNAs are short generally 19-26 bp in length, ~22-nt non-coding RNA species that control the post-transcriptional genes level expressions. The RNA polymerase II helps in transcribing the gene profiling and to from primary miRNA in the nucleus (Kincaid et al., 2011; Kumar et al., 2018; Praianantathavorn et al., 2016; Sannigrahi et al., 2017; Yousef et al., 2009) . After the miRNA is formed, the hairpin complex is transcribed and recognized by Drosha (RNase III enzyme) and DGCR8 (dsRNA-binding protein). After recognition, the pre-miRNA is exported from the nucleus to cytoplasm with the help of enzyme exportin-5 and Ran (Ras-related Nuclear protein). Soon in the cytoplasm, the pre-miRNA is converted to mature RNA i.e. ds RNA by RNase III ribonuclease Dicer into duplex mature RNA (Girardi et al., 2018; Mishra et al., 2020; Yadav et al., 2014) . This duplex RNA with the help of RNA-induced silencing complex (RISC) targets the messenger RNA and degrades the translational activity. This complexity emerged between the 3' untranslated regions (UTR) of miRNA and the seed region of miRNA (2-7bp) helps in cleaving and blocking the translations (Duygu et al., 2020; Tahi, 2012) . The mechanism adopted by the virus to generate miRNA helps to altercate the relevant host. These miRNAs generated can be used to find pertinent targets. Subsequently, several studies have conferred and explored the miRNAs as antiviral therapies against HIV 1, HSV, Dengue, Influenza, and Hepatitis C (HCV) (Mallick et al., 2009; Yousef et al., J o u r n a l P r e -p r o o f 2009). The efficacy of miRNA based treatment is demonstrated and found promising in the case of HCV treatment (Mishra et al., 2020) . Currently, the prevention and treatment approaches to SARS-CoV2 are very limited due to complexity in their genome sequences. However, approaches like convalescent plasma therapy (CPT) has provided a promising result in this aspect (Duan et al., 2020) . Further, to determine the viral pathogenesis and the outcome of the infection, it is necessary to understand the host-pathogen interactions. Although the studies between the host innate immune system and the role of miRNAmediated RNA silencing in SARS-CoV2 infection has not been enlightened yet. So, our current approach of using Vmir analyzer in identifying the maximum possible miRNAs among the genome sequences in SARS-Cov2 will help to identify the potent and effective target that triggers the hostpathogen interactions. Further, the in-silico based techniques used here will be helpful to scrutinize the role of the human miRNA on SARS-CoV2. Thus, we collected SARS-CoV2 genomes from different geographical locations mainly from China, India, USA, Italy, South America (Jamaica) to identify the notable host miRNA targets identifications. Our analysis revealed various miRNAs accompanying the host-pathogen relationship which is necessary for its survival and replication. Further, the genes targeted by these miRNAs can provide a definitive insight into the pathways involved in infection. This approach of targeting the genomes will provide an innovative platform as an antiviral strategy and to understand the interactions network. The SARS-CoV-2 miRNA prediction was carried out using the complete genome sequence of SARS-CoV2 carrying accession number (NC_045512.2, MT435086, MT339041, MT066156 and hairpin was visualized in Vmir visualizer (Duygu et al., 2020; Sardar et al., 2020) . The potential hairpin-like structures were extracted as a candidate's miRNA precursor. Human miRNA sequences are available in the miRBase database (http://www.mirbase.org). This is dependent on the average length of microRNA (~22bp). Each of the Candidate precursor miRNAs was searched for the nucleotide similarity with all human miRNAs by using the SSEARCH (it is useful for finding a short sequence within the library of miRNA) menu of the miRbase database. According to the principal, each of the input viral hairpin segments was aligned with all of the microRNAs in the miRBase then the highly similar were identified as target miRNA. The mature duplex microRNA consisted of two strands of microRNA that are complementary to each other. The complementary strand of the target miRNA can be complementary to the input viral sequences and along with the hybridization between the viral gene fragments and complementary template of the potential miRNA, which will analyzed and sorted out using RNA hybrid (Hasan et al., 2014) . Energetically most favorable hybridization between target microRNA and viral miRNA was predicted by the RNA hybrid tool (http://bibiserv.techfak.uni-bielefeld.de/rnahybrid). The results of RNA hybrid were categorized in terms of pairing energy (minimum free energy) and hybridization pattern. Four types of hybridization patterns were obtained from RNA hybrid analysis including 5' canonical, 3' compensatory, 5' seed, and ineffective hybridization (Hasan et al., 2014; Yadav et al., 2014) . Mature miRNAs were identified from pre-miRNAs sequences using Mature Bayes (http://mirna.imbb.forth.gr/MatureBayes.html), an online tool that uses Naive Bayes Classifier (NBC) taking into account of the sequences as well as structural information of experimental predicted miRNA precursors. All the potential pre-miRNAs identified by ViralmiR was used for analysis (Gkirtzou et al., 2010) . According to the microRNA target prediction principle which requires the sufficient base pairing between the miRNA and target miRNAs that can be classified into 5' canonical, 3' compensatory, 5' seed, and ineffective hybridization. The 5' dominant classes of target sites can be divided into 2 subtypes: 5' seed and 5' canonical. Both indicate the effective base pairing within 2 nd to 8 th position from the 5' end of miRNA. For 3' compensatory pattern, the candidate miRNA should show half sequence from middle to 3' end of miRNA that will perfectly match with miRNA. Pairing energy or minimum free energy (mef) indicating the stability of the hybridization. For the selection of potential miRNA, the pairing energy at -10 kcal/mol was utilized as a cut-off score. The miRNA targeting SARS-CoV2 genes with effective hybridization patterns (5' canonical, 3' compensatory, 5' seed) and minimum free energy of -10 kcal/mol were selected as potential miRNA (Hasan et al., 2014; Sardar et al., 2020; Tahi, 2012; Yousef et al., 2009 ). The RNAfold web server (http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi) was used to predict the secondary structure of pre-miRNA. This program is used to predict the most stable secondary structure of SARS-CoV2 hairpin sequences. The sequence applied for prediction analysis included pre-miRNA about 200bp upstream and about 100bp downstream flanking sequences at each end of the precursor. In all the cases, folding structures with minimal free energy were depicted (Lorenz et al., 2016; Wang et al., 2019) . Gene ontology analysis of the retrieved target genes was performed using PANTHER (Protein Analysis through Evolutionary Relationships) classification system (http://www.pantherdb.org) to gain insight into the molecular function, biological process and cellular components of the target genes products. Gene ID's of target genes were used for this analysis to find GO terms related to gene products (Chen and Glover, 2016; Dalmer and Clugston, 2019; Hu et al., 2019) . MT507794.1 (Jamaican) obtained a customized hairpin miRNAs structure which was found to be 52 (China, INDIA, US), 53 (Italy) and 51 (Jamaica) respectively (Hasan et al., 2014) . Additionally, the Vmir analyzer predicted all possible hairpins either in direct orientation or in reverse orientation according to the location (x-axis) and Vmir score (y-axis). The sequences of all the five viral SARS-CoV2 miRNA candidate precursors obtained was searched using SSEARCH menu the miRbase database for finding the nucleotide similarity with all human microRNAs using the default parameters (E-cut off value: 10, Maximum no. of hits: 100, Specific organism: human). This identified 22 sequences as a candidate miRNA precursor among the five viral miRNAs, based on the sequence similarity with human miRNAs. Subsequently, human miRNAs carrying a minimum of 19bp sequence similarity with candidate miRNA precursors were selected as prime target miRNAs as depicted in Table 1 . Similarly, the possible interaction shared between the human miRNAs and SARS-CoV2 viral miRNAs was shown with the help of diverging radial interaction plot Figure 3 . Additionally, miRNAs seed region (2-7) were selected based on 3' untranslated region of candidate miRNA precursor which was used to find the potential miRNA targets. In general; to depict a picture of gene silencing activity, the perfect complementarity between J o u r n a l P r e -p r o o f region 3' untranslated (UTR) of the miRNA and the (2-7 bp) seed region of miRNA is significant. Additionally, the precursor miRNA hairpins were also classified accordingly based on directions respectively as MD (forward direction) and MR (the reverse direction). For effective hybridization between target human miRNAs and precursor miRNAs of SARS-CoV2, we implemented the RNAhybrid (http://bibiserv2.cebitec.uni-bielefeld.de/rnahybrid/) tool. It is a widely used miRNAs target prediction tool based on parameters depicted as the minimum free energy hybridization shared between the long and short RNAs. The minimum free energy (MFE) or pairing energy specifies the stability of hybridization. Along with, the potential miRNA carrying pairing energy of -10 kcal/mol is essential for an effective hybridization. According to the criteria obtained for selection, we found 4 effective pairings hsa-mir-1267, hsa-mir-325, hsa-mir-5683, hsa-mir-3064-5p among the China; 6 effective pairings hsa-mir-5197, hsa-mir-147b, hsa-mir-6874, hsa-mir-138-1, hsa-mir-664b, hsa-mir-1-3p among India, 5 including hsa-mir-1267, hsa-mir-6874, hsa-mir-4425, hsamir-138-1, hsa-mir-1-3p among the US; 4 interactions hsa-mir-1267, has-mir-138-1, has-mir-664b, hsa-mir-1-3p among Italy and finally 3 collective pairing hsa-mir-1267, hsa-mir-325, hsa-mir-5683 among Jamaica particularly the human miRNAs involved with viral SARS-CoV2 miRNAs Figure 5 . After initial identification and validation, pre-miRNAs sequences were subjected to Mature Bayes for retrieving the small mature miRNAs. Large pre-miRNAs sequences were cleaved to short mature miRNAs of 22 nucleotide length. Mature miRNAs obtained were based on 5' and 3' stem location where we retrieved 10 mature miRNA from 5 pre miRNA from China; 12 mature miRNAs from 6 pre-miRNAs sequences of India, 10 mature miRNAs from 5 pre-miRNAs sequences of the US, 8 mature miRNAs from 4 pre-miRNAs of Italy and 6 mature miRNA from 3 pre-miRNA of Jamaica Table 2 . Depending on the assembly of RISC complex, one or both strands can serve as mature miRNA molecules, which was kept for further analysis. The RNAfold web server was used to predict the secondary structure of pre-miRNA. Expanded sequence dependence thermodynamic parameters were used to observe the hairpin structure of SARS-CoV2 (Mathews et al., 1999) . Further, the minimum free energy (mfe) and partition function along with avoid isolated base pairs (isolated base pairs predicted structures contains isolated base pairs helices of length 1 which may be undesirable) were selected as a fold algorithm. Subsequently, the RNAfold program was able to predict the most stable and standardized secondary structure. Target prediction by miRDB for all the viral SARS-CoV2 miRNAs gave detailed and precise data of all the collection of genes that can be targeted in the human genome. Basically, the server uses the MirTarget algorithm, which is based on a 7-mer seeding approach and further predicts the miRNAs targets among 3' UTRs of the human genes. Thus, we selected target genes with miRDB score >80 because a predicted target with a prediction score >80 is most likely to be real and not required any other supporting shreds of evidence. Finally, we proceeded through the PANTHER classification system to analyze the gene ontology and their involvement in different clusters of molecular functions, cellular components, and biological processes. This clustering was helpful to prove the significant determinant target gene involved in molecular functions, cellular components, and biological processes supplementary (file 1 & 2). The SARS-CoV2 has become a serious public threat globally and has drawn much attention because of its associations with mortality and the development of vaccines (Ahmed et al., 2020; Shanmugaraj et al., 2020) . Despite this, much effort has been carried out in studying the detailed mechanism of the SARS-CoV2 pathogenesis. Many in-silico works were reported regarding this aspect but a similar bioinformatics approach was carried out where it involved in identifying the potential miRNAs for MERS-CoV. This study ruled out the possibility of the MERS-CoV's miRNAs in locating significant importance to human miRNAs but the study does not determine the gene involvement with the pathogenic conditions (Hasan et al., 2014) . Further, to gain empathy and understand more about the nature of pathogenesis; identification of the novel miRNAs along with the gene-ontology encoded by SARS-CoV2 became one of the preliminary aims. In recent years, these miRNAs have emerged as an important intervention in biomedical research because of its involvement in several biological phenomena. So, studying virus-encoded miRNAs can therefore become an important platform to develop and understand the viral-host relationship which could be better aimed for therapeutics (Biology et al., 2016; Duygu et al., 2020; Ghosh et al., 2007; Girardi et al., 2018; Grundhoff, n.d.; Tahi, 2012; Yadav et al., 2014; Yousef et al., 2009 ). It was already reported from the studies that the human miRNAs target viral genes and functions as antiviral mediators to suppress the viral pathogenesis. By silencing those genes, human miRNAs ensure the prevention of complicated events. A study was carried out in ZIKA virus infection where it suggested that the miRNAs in astrocytes are dysregulated followed by the dysregulation of host miRNAs during the infection; which in turn halts the host gene functions (Karim, 2019) . On the contrary, one of the most exciting aspects comes from the virus-encoded miRNAs, providing their role in the survival and targeting the host genes where it successfully activates the host immune system (Bernier and Sagan, 2018) . In this regard, the host has developed various defenses against viruses such as the miRNA mediated host gene silencing is one of the strategies. Silencing the host genes provides the virus to invade the defenses adopted by the host, to replicate within the host and to avoid antiviral strategies (Russo and Potenza, 2011) . Further, to investigate whether the SARS-CoV2 J o u r n a l P r e -p r o o f can mediate and target the host genes we adopted this approach to find the potential role of the host genes in virus survival and replications. Initially, 1600 pre-miRNAs were obtained by Vmir (Grundhoff et al., 2006) . Previously, the miRNAs encoded in Epstein Barr Virus (EBV) were reported to be 24, which was obtained by Vmir (Griffiths-Jones et al., 2008) . Even, some ZIKA miRNAs targeted genes were identified which was found to be associated with the cell cycle process, cell communication, immune system regulation, etc. Halting these processes can have an immediate response in the viral transport, immune evasion which was found to be essential in ZIKA virus replications (Cristina et al., 2016; Karim, 2019) . Takin into account, we utilized a series of bioinformatics tools, which predicted 22 viral miRNAs against the human genome and based on our computational investigations, we hypothesize those miRNAs having a significant role in understanding the host-pathogen relationship. The human miRNAs hsa-mir-5197, hsa-mir-147b, hsa-mir-6874, hsa-mir-138-1, hsa-mir-664b, hsa-mir-1267, hsa-mir-4425, hsa-mir-1-3p, hsa-mir-5197, hsa-mir-325, hsa-mir-5683 and hsa-mir-3064-5p were found to show perfect complementarity and identity with the virus SARS-CoV2 miRNAs. To gain more insights into the associate functions with their gene's involvement; PANTHER classification was adopted to understand the relationship between the predicted miRNAs shared between the SARS-CoV2 and human. With respect to molecular functions, cluster target genes products were depicted to play a role in protein binding (GO:0005515), enzyme binding (GO:0019899), protein domain specific binding (GO:0019904), protein kinase activity (GO:0004672) and Rho-GTPase binding (GO: 0017048). These specific functions are required to lead to an abnormal state in the body. Further, the cellular response implicated by the gene (GO:0009987) implies the pathological damage. This has been reported earlier with Hepatitis C virus infection along with fatty liver conditions. However, in SARS-CoV2; only 2-11% infected patients carry liver morbidities but remains indistinct whether the liver damage in COVID-19 patients was caused by a viral infection or drug toxicity (Liu et al., n.d.) . The analysis can be based on the phenomenon that liver dysfunction can be impaired by SARS-CoV2 viral infection through the release of viral miRNA. Subsequently, biological process cluster classifications also predicted involvement of the genes linked to immune J o u r n a l P r e -p r o o f system (GO:0002520), response to cytokine (GO:0034097), biological adhesion (GO:0065007) and the regulation of signaling pathways (GO:0009966). These combined biological functions are necessary to implement a significant defense against viral infection. Finally, the cellular component classification suggested the target gene products i.e. immunoglobulin generation (GO:0019814), bounding to organelle (GO:0098588). This PANTHER classification system predicted the gene functions on a wide-scale but there occurs always a phenomena where some viruses might downregulate host genes expression to increase their gene expression either in the nucleus or in the cytoplasm (Chen et al., 2013; Dombkowski et al., n.d.; Mishra et al., 2020) . Another phenomenon associated is the basal transcription machinery of the SARS-CoV2 (TAF's) of TFIID complex which could prevent the RNA polymerase II to assemble on the promoters of the host genes at the initiation step that can specifically block the transcription. Therefore, the suppression of such genes could lead to a turnover and provide opportunities for viral miRNA to escape degradation. Similarly, the genes predicted in our study also bear some similar features with the genes associated with basal cell carcinoma which was primarily involved in biological adhesions and differentiation (Dombkowski et al., n.d.) . Additionally, a similar study was carried out to predict the interactions along the miRNA of MERS-CoV and human miRNAs hsa-miR-628-5p, hsa-miR-18a-3p, and hsa-miR-332-3p (Hasan et al., 2014) . Even, the hsa-mir-1-3p found in our study was also prominent in EboV ensuring the cellular signaling pathways and immune response dysregulations (Mishra et al., 2020) . Further, the functional annotation also carries some of the genes involved in signaling pathway (GO:0048014) where it has been reported that the signaling gene coding from the region CXCL16 for salmonella has been found to elicit cell-mediated immunity via IFN-ϒ production, further ARRB2 coding for signal gene plays an important role in inflammation evoking T-lymphocytes and cytokine production in lungs (Khan et al., 2020; Liu et al., n.d.) . Simultaneously, genes of highest enrichment (GO:0002520, GO: 0015031) were found to be similar to MERS-CoV and SARS-CoV, which was found to be associated as virus-encoded miRNAs in the metabolic process (Khan et al., 2020) . The miRNA's also act as a mediator against the viral evasion, thereby controlling the apoptosis (Cullen, 2013; Riley et al., 2012; Yang et al., 2013) . Particularly, GO:0042981 a gene predicted from our study was found to be involved in the regulation of apoptosis, where it was originated to be in close association with the J o u r n a l P r e -p r o o f infected cell. Additionally, the SARS-CoV2 also initiates various pathways; among them is the Rho-GTPase pathways which have been found to increase the vulnerability of infection among the infected host (Cullen, 2013; Heidbreder et al., 2010; Yang et al., 2013) . Similar gene GO:0017048 predicted and analysed in our study were found to be involved in the Rho-GTPase binding. Finally, the gene's (GO:0048014, GO:0032879, GO:0007186, GO:0034097) were predicted to be associated with various functions like regulation of Tie signaling pathways, regulation of localization, G-protein-couples receptor signaling and response to cytokine respectively. All these associated functions were found to be mostly involved in the regulation and survival of viral infection in the host. Apart from that these miRNAs encoded genes have been found to be linked in the prognosis assessment of the SARS-CoV2 infections (Lukassen et al., 2020; Schneider et al., 2020; Vigorito et al., 2007; Zhou et al., 2014) . Additionally, the human miRNAs hsa-mir-1267, hsa-mir-1-3p, hsa-mir-5683 were found prominent in all the five genomes of SARS-CoV2 predicted miRNAs; where the human miRNA hsa-mir-1267 shares the similarities with China, Italy, US and Jamaica; human miRNA hsa-mir-1-3p with India, US, Italy and hsa-mir-5683 with China and Jamaica respectively based on nucleotide similarity and gene ontology. These also define and predict the sequence similarities shared between the five SARS-CoV2 genomes being identified. So, our analysis resulted in the prediction of 44 essential genes involved with the biological functions, 13 with the molecular functions most prominently in protein binding and 13 genes with cellular components where one recognizing the immunoglobulin complex. Basically, the study carried out gave outline information about the dependence and involvement in survival and replication of SARS-CoV2 inside the human host. Further the insight analysis of miRNAs shared between the SARS-CoV2 and human also came into limelight. In conclusion, we propose the virus-host relationship in terms of miRNAs and its associated genes which may involve the upregulation or downregulations accompanying the dependent mechanisms in survival and replication of viruses. Based on data obtained in our computational approach and its analysis, we designed a mechanism of SARS-CoV2 pathogenesis through miRNA associated pathways. After the entry into the host cell cytoplasm, the SARS-CoV2 releases its genomic RNA, which utilizes the host machinery for its J o u r n a l P r e -p r o o f survival and replications Figure 5 . Probably the possible interactions which were predicted in our study may be helpful in mediating the association pathways in the host which accelerates the pathogenic conditions. Simultaneously, a phylogenetic tree was constructed to show the sequence similarities among the data obtained from the miRNA sequences of the SARS-CoV2 figure 6. In this study, we portrayed the pathogenesis of SARS-CoV2 through miRNAs thereby mediating a generalized picture of the putative miRNAs shared between the human genome and SARS-CoV2. The analysis was further carried out to identify the role of the essential genes involved in the survival and replication of the virus inside the host cytoplasm. Thereby our analysis predicted several novel miRNAs of SARS-CoV2 and as expected, the genes targeted by the miRNAs were found to be mostly involved in various biological and molecular functions. Additionally, the predicted 22 SARS-CoV2 miRNAs linked or shared with 12 human miRNAs can be utilized to understand the relationship with the virus. A much-needed work has to be carried out in this area but still, our findings provided a new insight to understand viral pathology. Last but not the least, a detailed experimental evidence is required to provide an assessment of these identified miRNAs to understand the viral replications more in-depth. None. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Preliminary identification of potential vaccine targets for the COVID-19 Coronavirus (SARS-CoV-2) Based on SARS-CoV Immunological Studies Immunoinformatics-aided identification of T cell and B cell epitopes in the surface glycoprotein of 2019-nCoV The diverse roles of microRNAs at the host-virus interface Development of epitope-based peptide vaccine against novel coronavirus 2019 (SARS-COV-2): Immunoinformatics approach The International Journal of Biochemistry Different microRNA alterations contribute to diverse outcomes following EV71 and CA16 infections : Insights from high-throughput sequencing in rhesus monkey peripheral blood mononuclear cells Divergent MicroRNA Targetomes of Closely Related Circulating Strains of a Polyomavirus Large-scale gene function analysis with PANTHER Classification System Genome-wide Prediction of microRNAs in Zika virus Genomes Reveals Possible Interactions with Human Genes review MicroRNAs as mediators of viral evasion of the immune system Gene ontology enrichment analysis of congenital diaphragmatic hernia-associated genes Effectiveness of convalescent plasma therapy in severe COVID-19 patients Computational analysis of microRNA-mediated interactions in SARS-CoV-2 infection Multi-epitope vaccine design using an immunoinformatics approach for 2019 novel coronavirus in China (SARS-CoV-2) miRNomics -The bioinformatics of microRNA genes On the Importance of Host MicroRNAs During MatureBayes: A probabilistic algorithm for identifying the mature miRNA within novel precursors miRBase: Tools for microRNA genomics Chapter 8 Computational Prediction of Viral miRNAs 721 A combined computational and microarray-based approach identifies novel microRNAs encoded by human gamma-herpesviruses A Computational Approach for Predicting Role of Human MicroRNAs in MERS-CoV Genome Downregulation of the Rho GTPase signaling pathway is involved in the microRNA-138-mediated inhibition of cell migration and invasion in tongue squamous cell carcinoma 512 In silico integration approach reveals key MicroRNAs and their target genes in follicular thyroid carcinoma Design of a Peptide-Based Subunit Vaccine against In silico analysis revealed Zika virus miRNAs associated with viral pathogenesis through alteration of host genes involved in immune response and neurological functions 1-11 Epigenetic regulator miRNA pattern differences among SARS-CoV , SARS-CoV-2 and SARS-CoV-2 world-wide isolates delineated the mystery behind the epic pathogenicity and distinct clinical characteristics of pandemic COVID-19 * Correspondence : Dr RNA virus microRNA that mimics a B-cell oncomiR crossm MicroRNA hsa-miR-324-5p Suppresses H5N1 Virus Replication Predicting RNA secondary structures from sequence and probing data SARS-CoV-2 receptor ACE 2 and TMPRSS 2 are primarily expressed in bronchial transient secretory cells 1-15 MicroRNome Analysis Unravels the Molecular Basis of SARS Infection in Bronchoalveolar Stem Cells 4 Expanded Sequence Dependence of Thermodynamic Parameters Improves Prediction of RNA Secondary Structure The Interplay Between Viral-Derived miRNAs and Host Immunity During Infection Human microRNAs profiling in response to influenza A viruses 409-420 In silico B-cell and T-cell epitope-based vaccine designing against Chikungunya virus Unique epidemiological and clinical features of the emerging Coronavirus Disease 2019 (COVID-19) and Pregnancy: What obstetricians need to know Antiviral effects of human microRNAs and conservation of their target sites Role of Host miRNA Hsa-miR-139-3p in HPV-16 -Induced Carcinomas 23 Comparative analyses of SAR-CoV2 genomes from different geographical locations and other coronavirus family genomes reveals unique features potentially consequential to host-virus interaction and pathogenesis In silico identification of miRNAs and their target prediction from Japanese encephalitis 5 SARS-CoV-2 receptor ACE2 and TMPRSS2 are predominantly expressed in a transient secretory cell type Perspectives on monoclonal antibody therapy as potential therapeutic intervention for Coronavirus disease-19 (COVID-19) 1 + 1 + 1 + 1 = ? A Case Approach to Perioper. Drug-Drug Interact A fast ab-initio method for predicting miRNA precursors in genomes Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants Article microRNA-155 Regulates the Generation of Immunoglobulin Class-Switched Plasma Cells DMFold: A novel method to predict RNA secondary structure with pseudoknots based on deep learning and improved base pair maximization principle A review of the 2019 Novel Coronavirus (COVID-19) based on current evidence Computational Identification of MicroRNAs for Targeting Long and Short Segments of Lassa Virus Epstein-Barr Virus _ Encoded LMP1 Upregulates MicroRNA-21 to Promote the Resistance of Nasopharyngeal Carcinoma Cells to Cisplatin-Induced Apoptosis by Suppressing PDCD4 and Fas-L MINIREVIEW A study of microRNAs in silico and in vivo : bioinformatics approaches to microRNA discovery and target identification 276 Interplay between microRNAs and host pathogen recognition receptors ( PRRs ) signaling pathways in response to viral infection This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.J o u r n a l P r e -p r o o f