Cutting Edge: PDGF-DD Binding to NKp44 Costimulates TLR9 Signaling and Proinflammatory Cytokine Secretion in Human Plasmacytoid Dendritic Cells.

NKp44 is a human receptor originally found on activated NK cells, group 1 and group 3 innate lymphoid cells that binds dimers of platelet-derived growth factor D (PDGF-DD). NKp44 is also expressed on tissue plasmacytoid dendritic cells (PDCs), but NKp44-PDGF-DD interaction on PDCs remains unstudied. Engagement of NKp44 with PDGF-DD in vitro enhanced PDC secretion of IFN-α, TNF, and IL-6 in response to the TLR9 ligand CpG-ODN, but not TLR7/8 ligands. In tissues, PDCs were found in close contact with PDGF-DD-expressing cells in the high endothelial venules and epithelium of tonsils, melanomas, and skin lesions infected with Molluscum contagiosum. Recombinant PDGF-DD enhanced the serum IFN-α response to systemic HSV-1 infection in a humanized mouse model. We conclude that NKp44 integrates with TLR9 signaling to enhance PDC cytokine production. These findings may have bearings for immune responses to TLR9-based adjuvants, therapy for tumors expressing PDGF-DD, and infections with DNA viruses that induce PDGF-DD expression to enhance viral spread.

IL17RB and IL17REL Expression Are Associated with Improved Prognosis in HPV-Infected Head and Neck Squamous Cell Carcinomas.

Changes in the cellular secretome are implicated in virus infection, malignancy, and anti-tumor immunity. We analyzed the association between transcriptional signatures (TS) from 24 different immune and stromal cell types on the prognosis of HPV-infected and HPV-free head and neck squamous carcinoma (HNSCC) patients from The Cancer Genome Atlas (TCGA) cohort. We found that HPV-positive HNSCC patients have tumors with elevated immune cell TS and improved prognosis, which was specifically associated with an increased tumor abundance of memory B and activated natural killer (NK) cell TS, compared to HPV-free HNSCC patients. HPV-infected patients upregulated many transcripts encoding secreted factors, such as growth factors, hormones, chemokines and cytokines, and their cognate receptors. Analysis of secretome transcripts and cognate receptors revealed that tumor expression of IL17RB and IL17REL are associated with a higher viral load and memory B and activated NK cell TS, as well as improved prognosis in HPV-infected HNSCC patients. The transcriptional parameters that we describe may be optimized to improve prognosis and risk stratification in the clinic and provide insights into gene and cellular targets that may potentially enhance anti-tumor immunity mediated by NK cells and memory B cells in HPV-infected HNSCC patients.

A Transcriptional Signature of IL-2 Expanded Natural Killer Cells Predicts More Favorable Prognosis in Bladder Cancer.

Activation of natural killer (NK) cell function is regulated by cytokines, such as IL-2, and secreted factors upregulated in the tumor microenvironment, such as platelet-derived growth factor D (PDGF-DD). In order to elucidate a clinical role for these important regulators of NK cell function in antitumor immunity, we generated transcriptional signatures representing resting, IL-2-expanded, and PDGF-DD-activated, NK cell phenotypes and established their abundance in The Cancer Genome Atlas bladder cancer (BLCA) dataset using CIBERSORT. The IL-2-expanded NK cell phenotype was the most abundant in low and high grades of BLCA tumors and was associated with improved prognosis. In contrast, PDGFD expression was associated with numerous cancer hallmark pathways in BLCA tumors compared with normal bladder tissue, and a high tumor abundance of PDGFD transcripts and the PDGF-DD-activated NK cell phenotype were associated with a poor BLCA prognosis. Finally, high tumor expression of transcripts encoding the activating NK cell receptors, KLRK1 and the CD160-TNFRSF14 receptor-ligand pair, was strongly correlated with the IL-2-expanded NK cell phenotype and improved BLCA prognosis. The transcriptional parameters we describe may be optimized to improve BLCA patient prognosis and risk stratification in the clinic and potentially provide gene targets of therapeutic significance for enhancing NK cell antitumor immunity in BLCA.

Bidirectional promoters: an enigmatic genome architecture and their roles in cancers.

Bidirectional promoters are the transcription regulatory regions of genes positioned head-to-head on opposite strands. Specific sequence signals, chromatin modifications and three-dimensional structures of the transcription site facilitate the unconventional yet tightly regulated transcription proceeding in both directions from these promoters. Mutations or aberrant epigenetic changes can lead to abnormal enhanced or reduced expression from either of the bidirectionally transcribed genes resulting in tumorigenesis. Moreover, bidirectionally transcribed genes might also contribute towards the immune regulation in tumor microenvironment. In this review, we aimed to expound the characteristic features of bidirectional promoters alongside their transcriptional regulations, and ultimately, the association of these enigmatic genomic elements in different cancers.

SARS-CoV-2 mutations altering regulatory properties: Deciphering host's and virus's perspectives.

Since the first recorded case of the SARS-CoV-2, it has acquired several mutations in its genome while spreading throughout the globe. In this study, we investigated the significance of these mutations by analyzing the host miRNA binding and virus's internal ribosome entry site (IRES). Strikingly, we observed that due to the acquired mutations, five host miRNAs lost their affinity for targeting the viral genome, and another five can target the mutated viral genome. Moreover, functional enrichment analysis suggests that targets of both of these miRNAs might be involved in various host immune signaling pathways. Remarkably, we detected that three particular mutations in the IRES can disrupt its secondary structure which can consequently make the virus less functional. These results could be valuable in exploring the functional importance of the mutations of SARS-CoV-2 and could provide novel insights into the differences observed different parts of the world.

Mutational profiles of marker genes of cervical carcinoma in Bangladeshi patients.

BACKGROUND: Cervical cancer is a gynecologic cancer type that develops in the cervix, accounting for 8% mortality of all female cancer patients. Infection with specific human papillomavirus (HPV) types is considered the most severe risk factor for cervical cancer. In the context of our socioeconomic conditions, an increasing burden of this disease and high mortality rate prevail in Bangladesh. Although several researches related to the epidemiology, HPV vaccination, and treatment modalities were conducted, researches on the mutation profiles of marker genes in cervical cancer in Bangladesh remain unexplored. METHODS: In this study, five different genomic regions within the top three most frequently mutated genes (EGFR, KRAS and PIK3CA) in COSMIC database with a key role in the development of cervical cancers were selected to study the mutation frequency in Bangladeshi patients. In silico analysis was done in two steps: nucleotide sequence analysis and its corresponding amino acid analysis. RESULTS: DNA from 46 cervical cancer tissue samples were extracted and amplified by PCR, using 1 set of primers designed for EGFR and 2 sets of primers designed for two different regions of both PIK3CA and KRAS gene. In total, 39 mutations were found in 26 patient samples. Eleven different mutations (23.91%), twenty-four different mutations (52.17%) and four mutations (8.7%) were found in amplified EGFR, PIK3CA and KRAS gene fragments, respectively; among which 1 (EGFR) was common in seven patient samples and 2 (PIKCA) were found in more than 1 patient. Our study shows that except for KRAS, the frequency of observed mutations in our patients is higher than those reported earlier in other parts of the world. Most of the exonic mutations were found only in the PIK3CA and EGFR genes. CONCLUSIONS: The study can be used as a basis to build a mutation database for cervical cancer in Bangladesh with the possibility of targetable oncogenic mutations. Further explorations are needed to establish future diagnostics, personalized medicine decisions, and other pharmaceutical applications for specific cancer subtypes.

Computational analysis revealed miRNAs produced by Chikungunya virus target genes associated with antiviral immune responses and cell cycle regulation.

Chikungunya virus (CHIKV) that causes chikungunya fever, is an alphavirus that belongs to the Togaviridae family containing a single-stranded RNA genome. Mosquitoes of the Aedes species act as the vectors for this virus and can be found in the blood, which can be passed from an infected person to a mosquito through mosquito bites. CHIKV has drawn much attention recently because of its potential of causing an epidemic. As the detailed mechanism of its pathogenesis inside the host system is still lacking, in this in silico research we have hypothesized that CHIKV might create miRNAs, which would target the genes associated with host cellular regulatory pathways, thereby providing the virus with prolonged refuge. Using bioinformatics approaches we found several putative miRNAs produced by CHIKV. Then we predicted the genes of the host targeted by these miRNAs. Functional enrichment analysis of these targeted genes shows the involvement of several biological pathways regulating antiviral immune stimulation, cellular proliferation, and cell cycle, thereby provide themselves with prolonged refuge and facilitate their pathogenesis, which in turn may lead to disease conditions. Finally, we analyzed a publicly available microarray dataset (GSE49985) to determine the altered expression levels of the targeted genes and found genes associated with pathways such as cell differentiation, phagocytosis, T-cell activation, response to cytokine, autophagy, Toll-like receptor signaling, RIG-I like receptor signaling and apoptosis. Our finding presents novel miRNAs and their targeted genes, which upon experimental validation could facilitate in developing new therapeutics to combat CHIKV infection and minimize CHIKV mediated diseases.

Transcriptome of nasopharyngeal samples from COVID-19 patients and a comparative analysis with other SARS-CoV-2 infection models reveal disparate host responses against SARS-CoV-2.

BACKGROUND: Although it is becoming evident that individual's immune system has a decisive influence on SARS-CoV-2 disease progression, pathogenesis is largely unknown. In this study, we aimed to profile the host transcriptome of COVID-19 patients from nasopharyngeal samples along with virus genomic features isolated from respective host, and a comparative analyses of differential host responses in various SARS-CoV-2 infection systems. RESULTS: Unique and rare missense mutations in 3C-like protease observed in all of our reported isolates. Functional enrichment analyses exhibited that the host induced responses are mediated by innate immunity, interferon, and cytokine stimulation. Surprisingly, induction of apoptosis, phagosome, antigen presentation, hypoxia response was lacking within these patients. Upregulation of immune and cytokine signaling genes such as CCL4, TNFA, IL6, IL1A, CCL2, CXCL2, IFN, and CCR1 were observed in lungs. Lungs lacked the overexpression of ACE2 as suspected, however, high ACE2 but low DPP4 expression was observed in nasopharyngeal cells. Interestingly, directly or indirectly, viral proteins specially non-structural protein mediated overexpression of integrins such as ITGAV, ITGA6, ITGB7, ITGB3, ITGA2B, ITGA5, ITGA6, ITGA9, ITGA4, ITGAE, and ITGA8 in lungs compared to nasopharyngeal samples suggesting the possible way of enhanced invasion. Furthermore, we found comparatively highly expressed transcription factors such as CBP, CEBP, NFAT, ATF3, GATA6, HDAC2, TCF12 which have pivotal roles in lung injury. CONCLUSIONS: Even though this study incorporates a limited number of cases, our data will provide valuable insights in developing potential studies to elucidate the differential host responses on the viral pathogenesis in COVID-19, and incorporation of further data will enrich the search of an effective therapeutics.

Genotypic distribution and prevalence of human papillomavirus infection in an apparently healthy female population in Bangladesh.

Objective: Human papillomavirus (HPV) comprises around 120 genotypically related viruses, classified into low- and high-risk HPVs, which are capable of replicating inside the keratinocytes of skin or mucous membranes. Studies suggest that infections with HPV-16 or HPV-18 have a higher rate of developing cancer. The aim of our study was to detect HPV early, and to estimate the genotype-specific prevalence of HPV in apparently healthy and asymptomatic females in Bangladesh. Method: After cervical swab specimen collection, a VIA test was performed to identify any type of abnormality in the cervix. A multiplex PCR amplification of HPV DNA, using L1 consensus primer systems, was performed with type-specific primers, followed by sequencing to detect HPV genotypes. Result: Of the 417 females, 121 were found to be HPV positive. The most prevalent high-risk HPV genotypes were found to be HPV-16 and HPV-18. Different patient demographic parameters, such as age, socioeconomic status, education, and history of first intercourse, were also studied to establish correlations with HPV infection. Conclusion: Our results might provide some insights into factors that influence the development of cervical cancer. They might also help in guiding better patient management, increased public health awareness, further testing, and the implementation of existing vaccines.

Computational engineering the binding affinity of Adalimumab monoclonal antibody for designing potential biosimilar candidate.

Amongst the anti-TNF-α therapy for rheumatoid arthritis and other autoimmune diseases, Adalimumab mAb is one of the best candidates. However, several risk factors are found to be associated with higher doses. Improvement of the binding properties will therefore significantly increase its therapeutic efficacy, reduce the dosage requirements, and ultimately the associated toxicity and treatment cost. Here, we proposed a systematic in silico approach of finding newer mAb variants with improved binding properties. Using various bioinformatics tools, we have identified the significant amino acid residues on Adalimumab mAb. Next, we searched for the suitability of the other residues for mutating the significant residues and from the combinations of suitable mutations, variants were designed. To find the most significant ones, binding properties of the variants were compared with the wild type Adalimumab mAb using molecular docking scrutiny and molecular dynamics simulation. Finally, structural properties between the variant and wild type were analyzed. We have identified the six most significant residues on Adalimumab mAb involved in the antigen-antibody interactions. Using the suitable mutations replacing each of these residues, we have modeled 143 variants. From several docking analyses, we have found five significant variants and after molecular dynamics simulation, one most significant variant with improved binding affinity was identified whose structural properties are similar to the wild type Adalimumab mAb. Designed variant from this study, may provide newer insights on the structure-based affinity improvements of monoclonal antibodies and likewise modifications of the Fc region will also improve the therapeutic effector functions of antibodies too.

Perversely expressed long noncoding RNAs can alter host response and viral proliferation in SARS-CoV-2 infection.

BACKGROUND: Regulatory roles of long noncoding RNAs (lncRNAs) during viral infection has become more evident in last decade, but are yet to be explored for SARS-CoV-2. MATERIALS & METHODS: We analyzed RNA-seq dataset of SARS-CoV-2 infected lung epithelial cells to identify differentially expressed genes. RESULTS: Our analyses uncover 21 differentially expressed lncRNAs broadly involved in cell survival and regulation of gene expression. These lncRNAs can directly interact with six differentially expressed protein-coding genes, and ten host genes that interact with SARS-CoV-2 proteins. Also, they can block the suppressive effect of nine microRNAs induced in viral infections. CONCLUSION: Our investigation determines that deregulated lncRNAs in SARS-CoV-2 infection are involved in viral proliferation, cellular survival, and immune response, ultimately determining disease outcome.

Lung transcriptome of a COVID-19 patient and systems biology predictions suggest impaired surfactant production which may be druggable by surfactant therapy.

An incomplete understanding of the molecular mechanisms behind impairment of lung pathobiology by COVID-19 complicates its clinical management. In this study, we analyzed the gene expression pattern of cells obtained from biopsies of COVID-19-affected patient and compared to the effects observed in typical SARS-CoV-2 and SARS-CoV-infected cell-lines. We then compared gene expression patterns of COVID-19-affected lung tissues and SARS-CoV-2-infected cell-lines and mapped those to known lung-related molecular networks, including hypoxia induced responses, lung development, respiratory processes, cholesterol biosynthesis and surfactant metabolism; all of which are suspected to be downregulated following SARS-CoV-2 infection based on the observed symptomatic impairments. Network analyses suggest that SARS-CoV-2 infection might lead to acute lung injury in COVID-19 by affecting surfactant proteins and their regulators SPD, SPC, and TTF1 through NSP5 and NSP12; thrombosis regulators PLAT, and EGR1 by ORF8 and NSP12; and mitochondrial NDUFA10, NDUFAF5, and SAMM50 through NSP12. Furthermore, hypoxia response through HIF-1 signaling might also be targeted by SARS-CoV-2 proteins. Drug enrichment analysis of dysregulated genes has allowed us to propose novel therapies, including lung surfactants, respiratory stimulants, sargramostim, and oseltamivir. Our study presents a distinct mechanism of probable virus induced lung damage apart from cytokine storm.

Aberration of the modulatory functions of intronic microRNA hsa-miR-933 on its host gene ATF2 results in type II diabetes mellitus and neurodegenerative disease development.

BACKGROUND: MicroRNAs are ~ 22-nucleotide-long biological modifiers that act as the post-transcriptional modulator of gene expression. Some of them are identified to be embedded within the introns of protein-coding genes, these miRNAs are called the intronic miRNAs. Previous findings state that these intronic miRNAs are co-expressed with their host genes. This co-expression is necessary to maintain the robustness of the biological system. Till to date, only a few experiments are performed discretely to elucidate the functional relationship between few co-expressed intronic miRNAs and their associated host genes. RESULTS: In this study, we have interpreted the underlying modulatory mechanisms of intronic miRNA hsa-miR-933 on its target host gene ATF2 and found that aberration can lead to several disease conditions. A protein-protein interaction network-based approach was adopted, and functional enrichment analysis was performed to elucidate the significantly over-represented biological functions and pathways of the common targets. Our approach delineated that hsa-miR-933 might control the hyperglycemic condition and hyperinsulinism by regulating ATF2 target genes MAP4K4, PRKCE, PEA15, BDNF, PRKACB, and GNAS which can otherwise lead to the development of type II diabetes mellitus. Moreover, we showed that hsa-miR-933 can regulate a target of ATF2, brain-derived neurotrophic factor (BDNF), to modulate the optimal expression of ATF2 in neuron cells to render neuroprotection for the inhibition of neurodegenerative diseases. CONCLUSIONS: Our in silico model provides interesting resources for experimentations in a model organism or cell line for further validation. These findings may extend the common perception of gene expression analysis with new regulatory functionality.

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.

A detailed understanding of the molecular mechanism of SARS-CoV-2 pathogenesis is still elusive, and there is a need to address its deadly nature and to design effective therapeutics. Here, we present a study that elucidates the interplay between the SARS-CoV and SARS-CoV-2 viruses' and host's miRNAs, an epigenetic regulator, as a mode of pathogenesis; and we explored how the SARS-CoV and SARS-CoV-2 infections differ in terms of their miRNA-mediated interactions with the host and the implications this has in terms of disease complexity. We have utilized computational approaches to predict potential host and viral miRNAs and their possible roles in different important functional pathways. We have identified several putative host antiviral miRNAs that can target the SARS viruses and also predicted SARS viruses-encoded miRNAs targeting host genes. In silico predicted targets were also integrated with SARS-infected human cell microarray and RNA-seq gene expression data. A comparison between the host miRNA binding profiles on 67 different SARS-CoV-2 genomes from 24 different countries with respective country's normalized death count surprisingly uncovered some miRNA clusters, which are associated with increased death rates. We have found that induced cellular miRNAs can be both a boon and a bane to the host immunity, as they have possible roles in neutralizing the viral threat; conversely, they can also function as proviral factors. On the other hand, from over representation analysis, our study revealed that although both SARS-CoV and SARS-CoV-2 viral miRNAs could target broad immune-signaling pathways; only some of the SARS-CoV-2 miRNAs are found to uniquely target some immune-signaling pathways, such as autophagy, IFN-I signaling, etc., which might suggest their immune-escape mechanisms for prolonged latency inside some hosts without any symptoms of COVID-19. Furthermore, SARS-CoV-2 can modulate several important cellular pathways that might lead to the increased anomalies in patients with comorbidities like cardiovascular diseases, diabetes, breathing complications, etc. This might suggest that miRNAs can be a key epigenetic modulator behind the overcomplications amongst the COVID-19 patients. Our results support that miRNAs of host and SARS-CoV-2 can indeed play a role in the pathogenesis which can be further concluded with more experiments. These results will also be useful in designing RNA therapeutics to alleviate the complications from COVID-19.

A switch in bidirectional histone mark leads to differential modulation of lincRNAs involved in neuronal and hematopoietic cell differentiation from their progenitors.

Long intergenic noncoding RNAs (lincRNAs) are more than 200 bases long, transcribed from intergenic genomic regions and do not undergo translation. They have regulatory roles in differentiation and development. However, how their transcription is activated and how their expression is differentially modulated in differentiation is quite unclear. In this study, we explored and analyzed data at the transcriptomic and epigenetic level to address these questions. Here, we identified novel lincRNAs that are differentially expressed in neuronal and hematopoietic differentiation and showed that such differential modulations are achieved under epigenetic regulations. lincRNAs that are upregulated in mature cells than in progenitor are activated from a bivalent poised state where activating H3K4me3/H3K9ac/H3K27ac and suppressive H3K9me3/H3K27me3 marks are colocalized. And, lincRNAs that are downregulated in mature cells after differentiation are suppressed by the addition of H3K9me3/H3K27me3 marks. Moreover, here we show a tissue-specific expression pattern of lincRNAs in various cell lines and normal tissues. The study reveals bidirectional histone marks as an epigenetic means of directing the differential expression of lincRNAs which are found to be involved in the process of cellular differentiation.

SARS-CoV-2 Proteins Exploit Host's Genetic and Epigenetic Mediators for the Annexation of Key Host Signaling Pathways.

The constant rise of the death toll and cases of COVID-19 has made this pandemic a serious threat to human civilization. Understanding of host-SARS-CoV-2 interaction in viral pathogenesis is still in its infancy. In this study, we utilized a blend of computational and knowledgebase approaches to model the putative virus-host interplay in host signaling pathways by integrating the experimentally validated host interactome proteins and differentially expressed host genes in SARS-CoV-2 infection. While searching for the pathways in which viral proteins interact with host proteins, we discovered various antiviral immune response pathways such as hypoxia-inducible factor 1 (HIF-1) signaling, autophagy, retinoic acid-inducible gene I (RIG-I) signaling, Toll-like receptor signaling, fatty acid oxidation/degradation, and IL-17 signaling. All these pathways can be either hijacked or suppressed by the viral proteins, leading to improved viral survival and life cycle. Aberration in pathways such as HIF-1 signaling and relaxin signaling in the lungs suggests the pathogenic lung pathophysiology in COVID-19. From enrichment analysis, it was evident that the deregulated genes in SARS-CoV-2 infection might also be involved in heart development, kidney development, and AGE-RAGE signaling pathway in diabetic complications. Anomalies in these pathways might suggest the increased vulnerability of COVID-19 patients with comorbidities. Moreover, we noticed several presumed infection-induced differentially expressed transcription factors and epigenetic factors, such as miRNAs and several histone modifiers, which can modulate different immune signaling pathways, helping both host and virus. Our modeling suggests that SARS-CoV-2 integrates its proteins in different immune signaling pathways and other cellular signaling pathways for developing efficient immune evasion mechanisms while leading the host to a more complicated disease condition. Our findings would help in designing more targeted therapeutic interventions against SARS-CoV-2.

Genomic and evolutionary features of two AHPND positive Vibrio parahaemolyticus strains isolated from shrimp (Penaeus monodon) of south-west Bangladesh.

BACKGROUND: Due to its rapid lethal effect in the early development stage of shrimp, acute hepatopancreatic necrosis disease (AHPND) has been causing great economic losses, since its first outbreak in southeast China in 2009. Vibrio parahaemolyticus, carrying the pirA and pirB toxin genes is known to cause AHPND in shrimp. The overall objective of this study was to sequence the whole genome of AHPND positive V. parahaemolyticus strains isolated from shrimp (Peneaus monodon) of the south-west region of Bangladesh in 2016 and 2017 and characterize the genomic features and emergence pattern of this marine pathogen. RESULTS: Two targeted AHPND positive V. parahaemolyticus strains were confirmed using PCR with 16S rRNA, ldh, AP3 and AP4 primers. The assembled genomes of strain MSR16 and MSR17 were comprised of a total of 5,393,740 bp and 5,241,592 bp, respectively. From annotation, several virulence genes involved in chemotaxis and motility, EPS type II secretion system, Type III secretion system-1 (T3SS-1) and its secreted effectors, thermolabile hemolysin were found in both strains. Importantly, the ~ 69 kb plasmid was identified in both MSR16 and MSR17 strains containing the two toxin genes pirA and pirB. Antibiotic resistance genes were predicted against ß-lactam, fluoroquinolone, tetracycline and macrolide groups in both MSR16 and MSR17 strains. CONCLUSIONS: The findings of this research may facilitate the tracking of pathogenic and/or antibiotic-resistant V. parahaemolyticus isolates between production sites, and the identification of candidate strains for the production of vaccines as an aid to control of this devastating disease. Also, the emergence pattern of this pathogen can be highlighted to determine the characteristic differences of other strains found all over the world.

In silico analysis revealed Zika virus miRNAs associated with viral pathogenesis through alteration of host genes involved in immune response and neurological functions.

BACKGROUND: The concurrent Zika Virus (ZIKV) outbreaks in the United States and Northeast Brazil have evoked global surveillance. Zika infection has been correlated with severe clinical symptoms, such as microcephaly, Guillain-Barré syndrome, and other congenital brain abnormalities. Recent data suggest that ZIKV predominantly targets neural progenitor cells leading to neurological impairment. Despite the clinical evidence, detailed experimental mechanism of ZIKV neurotropic pathogenesis has not been fully understood yet. Here we hypothesized that ZIKV produces miRNAs, which target essential host genes involved in various cellular pathways facilitating their survival through immune evasion and progression of disease during brain development. METHODS: From genome sequence information using several bioinformatic tools, we predicted pri-miRNAs, pre-miRNAs, and finally the mature miRNAs produced by ZIKV. We also identified their target genes and performed functional enrichment analysis to identify the biological processes associated with these genes. Finally, we analyzed a publicly available RNA-seq data set to determine the altered expression level of the targeted genes. RESULTS: From ZIKV genome sequence, we identified and validated 47 putative novel miRNAs. Functional enrichment of the targeted genes demonstrates the involvement of various biological pathways regulating cellular signaling, neurological functions, cancer, and fetal development. The expression analysis of these genes showed that ZIKV-produced miRNAs downregulate the key genes involved in these pathways, which in turn may lead to impaired brain development. CONCLUSIONS: Our finding proposes novel ZIKV miRNAs and their targets, which upon experimental validation could help developing new therapeutics to combat ZIKV infection and minimize ZIKV-mediated pathologies.