A combined computational and microarray-based approach identifies novel microRNAs encoded by human gamma-herpesviruses

Epigenetic Mechanisms in Latent Epstein-Barr Virus Infection and Associated Cancers

The Epstein-Barr Virus (EBV) is a double-stranded DNA-based human tumor virus that was first isolated in 1964 from lymphoma biopsies. Since its initial discovery, EBV has been identified as a major contributor to numerous cancers and chronic autoimmune disorders. The virus is particularly efficient at infecting B-cells but can also infect epithelial cells, utilizing an array of epigenetic strategies to establish long-term latent infection. The association with histone modifications, alteration of DNA methylation patterns in host and viral genomes, and microRNA targeting of host cell factors are core epigenetic strategies that drive interactions between host and virus, which are necessary for viral persistence and progression of EBV-associated diseases. Therefore, understanding epigenetic regulation and its role in post-entry viral dynamics is an elusive area of EBV research. Here, we present current outlooks of EBV epigenetic regulation as it pertains to viral interactions with its host during latent infection and its propensity to induce tumorigenesis. We review the important epigenetic regulators of EBV latency and explore how the strategies involved during latent infection drive differential epigenetic profiles and host-virus interactions in EBV-associated cancers.

Application of Integrated Computational Approaches in Prediction of Plant Virus Encoded miRNAs and Their Targeted Plant Genes

This chapter presents a comprehensive approach to predict novel miRNAs encoded by plant viruses and identify their target plant genes, through integration of various ab initio computational approaches. The predictive process begins with the analysis of plant viral sequences using the VMir Analyzer software. VMir Viewer software is then used to extract primary hairpins from these sequences. To distinguish real miRNA precursors from pseudo miRNA precursors, MiPred web-based software is employed. Verified real pre-miRNA sequences with a minimum free energy of < -20 Kcal/mol, are further analyzed using the RNAshapes software. Validation of predictions involves comparing them with available Expressed Sequence Tags (ESTs) from the relevant plant using BlastN. Short sequences with lengths ranging from 19 to 25 nucleotides and exhibiting <5 mismatches are prioritized for miRNA prediction. The precise locations of these short sequences within pre-miRNA structures generated using RNAshapes are meticulously identified, with a focus on those situated on the 5' and 3' arms of the structures, indicating potential miRNAs. Sequences within the arms of pre-miRNA structures are used to predict target sites within the ESTs of the specific plant, facilitated by psRNA Target software, revealing genes with potential regulatory roles in the plant. To confirm the outcome of target prediction, results are individually submitted to the RNAhybrid web-based software. For practical demonstration, this approach is applied to analyze African cassava mosaic virus (ACMV) and East African cassava mosaic virus-Uganda (EACMV-UG) viruses, as well as the ESTs of Jatropha and cassava.

Hepatitis B Virus and microRNAs: A Bioinformatics Approach

In recent decades, microRNAs (miRNAs) have emerged as key regulators of gene expression, and the identification of viral miRNAs (v-miRNAs) within some viruses, including hepatitis B virus (HBV), has attracted significant attention. HBV infections often progress to chronic states (CHB) and may induce fibrosis/cirrhosis and hepatocellular carcinoma (HCC). The presence of HBV can dysregulate host miRNA expression, influencing several biological pathways, such as apoptosis, innate and immune response, viral replication, and pathogenesis. Consequently, miRNAs are considered a promising biomarker for diagnostic, prognostic, and treatment response. The dynamics of miRNAs during HBV infection are multifaceted, influenced by host variability and miRNA interactions. Given the ability of miRNAs to target multiple messenger RNA (mRNA), understanding the viral-host (human) interplay is complex but essential to develop novel clinical applications. Therefore, bioinformatics can help to analyze, identify, and interpret a vast amount of miRNA data. This review explores the bioinformatics tools available for viral and host miRNA research. Moreover, we introduce a brief overview focusing on the role of miRNAs during HBV infection. In this way, this review aims to help the selection of the most appropriate bioinformatics tools based on requirements and research goals.

Insight into the Epigenetics of Kaposi's Sarcoma-Associated Herpesvirus

Epigenetic reprogramming represents a series of essential events during many cellular processes including oncogenesis. The genome of Kaposi's sarcoma-associated herpesvirus (KSHV), an oncogenic herpesvirus, is predetermined for a well-orchestrated epigenetic reprogramming once it enters into the host cell. The initial epigenetic reprogramming of the KSHV genome allows restricted expression of encoded genes and helps to hide from host immune recognition. Infection with KSHV is associated with Kaposi's sarcoma, multicentric Castleman's disease, KSHV inflammatory cytokine syndrome, and primary effusion lymphoma. The major epigenetic modifications associated with KSHV can be labeled under three broad categories: DNA methylation, histone modifications, and the role of noncoding RNAs. These epigenetic modifications significantly contribute toward the latent-lytic switch of the KSHV lifecycle. This review gives a brief account of the major epigenetic modifications affiliated with the KSHV genome in infected cells and their impact on pathogenesis.

MicroRNAs as powerful tool against COVID-19: Computational perspective

Severe acute respiratory syndrome coronavirus 2 is the virus that is responsible for the current pandemic, COVID-19 (SARS-CoV-2). MiRNAs, a component of RNAi technology, belong to the family of short, noncoding ssRNAs, and may be crucial in the battle against this global threat since they are involved in regulating complex biochemical pathways and may prevent viral proliferation, translation, and host expression. The complicated metabolic pathways are modulated by the activity of many proteins, mRNAs, and miRNAs working together in miRNA-mediated genetic control. The amount of omics data has increased dramatically in recent years. This massive, linked, yet complex metabolic regulatory network data offers a wealth of opportunity for iterative analysis; hence, extensive, in-depth, but time-efficient screening is necessary to acquire fresh discoveries; this is readily performed with the use of bioinformatics. We have reviewed the literature on microRNAs, bioinformatics, and COVID-19 infection to summarize (1) the function of miRNAs in combating COVID-19, and (2) the use of computational methods in combating COVID-19 in certain noteworthy studies, and (3) computational tools used by these studies against COVID-19 in several purposes. This article is categorized under: Infectious Diseases > Computational Models.

Virus-Derived Small RNAs and microRNAs in Health and Disease

MicroRNAs (miRNAs) are short noncoding RNAs that can regulate all steps of gene expression (induction, transcription, and translation). Several virus families, primarily double-stranded DNA viruses, encode small RNAs (sRNAs), including miRNAs. These virus-derived miRNAs (v-miRNAs) help the virus evade the host's innate and adaptive immune system and maintain an environment of chronic latent infection. In this review, the functions of the sRNA-mediated virus-host interactions are highlighted, delineating their implication in chronic stress, inflammation, immunopathology, and disease. We provide insights into the latest viral RNA-based research-in silico approaches for functional characterization of v-miRNAs and other RNA types. The latest research can assist toward the identification of therapeutic targets to combat viral infections.

Mir-155 expression is downregulated in hematopoietic stem cell transplantation patients with Epstein-Barr virus infection

BACKGROUND

MicroRNAs (miRNAs) are a class of short length double strand genome encoded RNAs that are produced to repress post-transcriptionally the expression of cellular mRNAs. 2578 unique mature miRNAs are currently annotated in the human genome and participate in the regulation of multiple events, such as cellular proliferation or apoptosis. The over-expression of miR-155 of cellular origin might play a key role in the life cycle of EBV. In this study 24 pediatric patients undergoing HSCT seropositive and seronegative to EBV were enrolled. Thirty-one peripheral blood samples were collected from these patients. The mir-155 expression profile has been evaluated by a stem-loop Real Time PCR in all these conditions.

METHODS

Of 24 patients, 4 were seronegative to EBV and EBV negative to PCR (Group I), 10 were seropositive to EBV and EBV negative to PCR (Group II) and 10 were seropositive to EBV and EBV positive to PCR (Group III).

RESULTS

Based on relative quantification, the mir-155 expression was compared among the groups. The comparison between HSCT patients without EBV infection seronegative to EBV (Group I) showed higher levels of mir-155 expression than patients seropositive to EBV (P=0.1419). The mir-155 expression levels in seronegative to EBV were not significantly different compared with the patients seropositive to EBV (P=0.6504). The mir-155 expression levels in seropositive to EBV without and with EBV infection (positive viral load), were not significantly (P=0.7667). Also, when we evaluated the mir-155 expression levels comparing all EBV negative patients with an active EBV infection, we did not observe a statistically significant difference (P=0.9782).

CONCLUSIONS

Our results are controversial, showing a higher production of mir-155 levels during EBV infection.

Cross-Kingdom Interaction of miRNAs and Gut Microbiota with Non-Invasive Diagnostic and Therapeutic Implications in Colorectal Cancer

Colorectal cancer (CRC) has one of the highest incidences among all types of malignant diseases, affecting millions of people worldwide. It shows slow progression, making it preventable. However, this is not the case due to shortcomings in its diagnostic and management procedure and a lack of effective non-invasive biomarkers for screening. Here, we discuss CRC-associated microRNAs (miRNAs) and gut microbial species with potential as CRC diagnostic and therapy biomarkers. We provide rich evidence of cross-kingdom miRNA-mediated interactions between the host and gut microbiome. miRNAs have emerged with the ability to shape the composition and dynamics of gut microbiota. Intestinal microbes can uptake miRNAs, which in turn influence microbial growth and provide the ability to regulate the abundance of various microbial species. In the context of CRC, targeting miRNAs could aid in manipulating the balance of the microbiota. Our findings suggest the need for correlation analysis between the composition of the gut microbiome and the miRNA expression profile.

A microRNA Arising from the Negative Strand of SARS-CoV-2 Genome Targets FOS to Reduce AP-1 Activity

Virus-encoded microRNAs were first reported in the Epstein-Barr virus in 2004. Subsequently, a few hundred viral miRNAs have been identified, mainly in DNA viruses belonging to the herpesviridae family. To date, only 30 viral miRNAs encoded by RNA viruses are reported by miRBase. Since the outbreak of the SARS-CoV-2 pandemic, several studies have predicted and, in some cases, experimentally validated miRNAs originating from the positive strand of the SARS-CoV-2 genome. By integrating NGS data analysis and qRT-PCR approaches, we found that SARS-CoV-2 also encodes for a viral miRNA arising from the minus (antisense) strand of the viral genome, in the region encoding for ORF1ab, herein referred to as SARS-CoV-2-miR-AS1. Our data show that the expression of this microRNA increases in a time course analysis of SARS-CoV-2 infected cells. Furthermore, enoxacin treatment enhances the accumulation of the mature SARS-CoV-2-miR-AS1 in SARS-CoV-2 infected cells, arguing for a Dicer-dependent processing of this small RNA. In silico analysis suggests that SARS-CoV-2-miR-AS1 targets a set of genes which are translationally repressed during SARS-CoV-2 infection. We experimentally validated that SARS-CoV-2-miR-AS1 targets FOS, thus repressing the AP-1 transcription factor activity in human cells.

DIANA-microT 2023: including predicted targets of virally encoded miRNAs

DIANA-microT-CDS is a state-of-the-art miRNA target prediction algorithm catering the scientific community since 2009. It is one of the first algorithms to predict miRNA binding sites in both the 3' Untranslated Region (3'-UTR) and the coding sequence (CDS) of transcripts, with increased performance. Its current version, DIANA-microT 2023 (www.microrna.gr/microt_webserver/), brings forward a significantly updated set of interactions. DIANA-microT-CDS has been executed utilizing annotation information from Ensembl v102, miRBase 22.1 and, for the first time, MirGeneDB 2.1, yielding more than 83 million interactions in human, mouse, rat, chicken, fly and worm species. Additionally, this version delivers predicted interactions of miRNAs encoded from 20 viruses against host transcripts from human, mouse and chicken species. Numerous resources have been interconnected into DIANA-microT, including DIANA-TarBase, plasmiR, HMDD, UCSC, dbSNP, ClinVar, as well as miRNA/gene abundance values for 369 distinct cell-lines/tissues. The server interface has been redesigned allowing users to use smart filtering options, identify abundance patterns of interest, pinpoint known SNPs residing on binding sites and obtain miRNA-disease information. The contents of DIANA-microT webserver are freely accessible and can also be locally downloaded without any login requirements.

In silico analysis of genomic landscape of SARS-CoV-2 and its variant of concerns (Delta and Omicron) reveals changes in the coding potential of miRNAs and their target genes

COVID-19 related morbidities and mortalities are still continued due to the emergence of new variants of SARS-CoV-2. In the last few years, viral miRNAs have been the centre of study to understand the disease pathophysiology. In this work, we aimed to predict the change in coding potential of the viral miRNAs in SARS-CoV-2's VOCs, Delta and Omicron compared to the Reference (Wuhan origin) strain using bioinformatics tools. After ab-intio based screening by the Vmir tool and validation, we retrieved 22, 6, and 6 pre-miRNAs for Reference, Delta, and Omicron. Most of the predicted unique pre-miRNAs of Delta and Omicron were found to be encoded from the terminal and origin of the genomic sequence, respectively. Mature miRNAs identified by MatureBayes from the unique pre-miRNAs were used for target identification using miRDB. A total of 1786, 216, and 143 high-confidence target genes were captured for GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis. The GO and KEGG pathways terms analysis revealed the involvement of Delta miRNAs targeted genes in the pathways such as Human cytomegalovirus infection, Breast cancer, Apoptosis, Neurotrophin signaling, and Axon guidance whereas the Sphingolipid signaling pathway was found for the Omicron. Furthermore, we focussed our analysis on target genes that were validated through GEO's (Gene Expression Omnibus) DEGs (Differentially Expressed Genes) dataset, in which FGL2, TNSF12, OGN, GDF11, and BMP11 target genes were found to be down-regulated by Reference miRNAs and YAE1 and RSU1 by Delta. Few genes were also observed to be validated among in up-regulated gene set of the GEO dataset, in which MMP14, TNFRSF21, SGMS1, and TMEM192 were related to Reference whereas ZEB2 was detected in all three strains. This study thus provides an in-silico based analysis that deciphered the unique pre-miRNAs in Delta and Omicron compared to Reference. However, the findings need future wet lab studies for validation.

MicroRNA Regulation of Human Herpesvirus Latency

Herpesviruses are ubiquitous human pathogens. After productive (lytic) infection, all human herpesviruses are able to establish life-long latent infection and reactivate from it. Latent infection entails suppression of viral replication, maintenance of the viral genome in infected cells, and the ability to reactivate. Most human herpesviruses encode microRNAs (miRNAs) that regulate these processes during latency. Meanwhile, cellular miRNAs are hijacked by herpesviruses to participate in these processes. The viral or cellular miRNAs either directly target viral transcripts or indirectly affect viral infection through host pathways. These findings shed light on the molecular determinants that control the lytic-latent switch and may lead to novel therapeutics targeting latent infection. We discuss the multiple mechanisms by which miRNAs regulate herpesvirus latency, focusing on the patterns in these mechanisms.

KSHV/HHV8-mediated hematologic diseases

Kaposi sarcoma (KS) herpesvirus (KSHV), also known as human herpesvirus 8, is the causal agent of KS but is also pathogenetically related to several lymphoproliferative disorders, including primary effusion lymphoma (PEL)/extracavitary (EC) PEL, KSHV-associated multicentric Castleman disease (MCD), KSHV+ diffuse large B-cell lymphoma, and germinotropic lymphoproliferative disorder. These different KSHV-associated diseases may co-occur and may have overlapping features. KSHV, similar to Epstein-Barr virus (EBV), is a lymphotropic gammaherpesvirus that is preferentially present in abnormal lymphoid proliferations occurring in immunecompromised individuals. Notably, both KSHV and EBV can infect and transform the same B cell, which is frequently seen in KSHV+ EBV+ PEL/EC-PEL. The mechanisms by which KSHV leads to lymphoproliferative disorders is thought to be related to the expression of a few transforming viral genes that can affect cellular proliferation and survival. There are critical differences between KSHV-MCD and PEL/EC-PEL, the 2 most common KSHV-associated lymphoid proliferations, including viral associations, patterns of viral gene expression, and cellular differentiation stage reflected by the phenotype and genotype of the infected abnormal B cells. Advances in treatment have improved outcomes, but mortality rates remain high. Our deepening understanding of KSHV biology, clinical features of KSHV-associated diseases, and newer clinical interventions should lead to improved and increasingly targeted therapeutic interventions.

EBV-associated diseases: Current therapeutics and emerging technologies

EBV is a prevalent virus, infecting >90% of the world's population. This is an oncogenic virus that causes ~200,000 cancer-related deaths annually. It is, in addition, a significant contributor to the burden of autoimmune diseases. Thus, EBV represents a significant public health burden. Upon infection, EBV remains dormant in host cells for long periods of time. However, the presence or episodic reactivation of the virus increases the risk of transforming healthy cells to malignant cells that routinely escape host immune surveillance or of producing pathogenic autoantibodies. Cancers caused by EBV display distinct molecular behaviors compared to those of the same tissue type that are not caused by EBV, presenting opportunities for targeted treatments. Despite some encouraging results from exploration of vaccines, antiviral agents and immune- and cell-based treatments, the efficacy and safety of most therapeutics remain unclear. Here, we provide an up-to-date review focusing on underlying immune and environmental mechanisms, current therapeutics and vaccines, animal models and emerging technologies to study EBV-associated diseases that may help provide insights for the development of novel effective treatments.

Are Argonaute-Associated Tiny RNAs Junk, Inferior miRNAs, or a New Type of Functional RNAs?

The biosynthesis pathways of microRNAs (miRNAs) have been well characterized with the identification of the required components. miRNAs are synthesized from the transcripts of miRNA genes and other RNAs, such as introns, transfer RNAs, ribosomal RNAs, small nucleolar RNAs, and even viral miRNAs. These small RNAs are loaded into Argonaute (AGO) proteins and recruit the effector complexes to target mRNAs, repressing their gene expression post-transcriptionally. While mature miRNAs were defined as 19–23 nucleotides (nt), tiny RNAs (tyRNAs) shorter than 19 nt have been found to bind AGOs as equivalent or lesser miRNAs compared to their full-length mature miRNAs. In contrast, my recent study revealed that when human AGO3 loads 14 nt cleavage-inducing tyRNAs (cityRNAs), comprised of the first 14 nt of their corresponding mature miRNA, it can become a comparable slicer to AGO2. This observation raises the possibility that tyRNAs play distinct roles from their mature form. This minireview focuses on human AGO-associated tyRNAs shorter than 19 nt and discusses their possible biosynthesis pathways and physiological benefits, including how tyRNAs could avoid target-directed miRNA degradation accompanied by AGO polyubiquitination.

CircEAF2 counteracts Epstein-Barr virus-positive diffuse large B-cell lymphoma progression via miR-BART19-3p/APC/β-catenin axis

Background

Epstein-Barr virus (EBV) represents an important pathogenic factor of lymphoma and is significantly associated with poor clinical outcome of diffuse large B-cell lymphoma (DLBCL). Circular RNAs (circRNAs) play an essential role in lymphoma progression. However, the underlying mechanism of circRNA on DLBCL progression related to EBV remains largely unknown.

Methods

CircRNA was screened by high-throughput sequencing in tumor samples of 12 patients with DLBCL according to EBV infection status. Expression of circEAF2, as well as the relationship with clinical characteristics and prognosis, were further analyzed in tumor samples of 100 DLBCL patients using quantitative real-time PCR. Gain- and loss-of-function experiments were conducted to investigate the biological functions of circEAF2 both in vitro and in vivo. The underlying mechanism of circRNA on DLBCL progression were further determined by RNA sequencing, RNA pull down assay, dual-luciferase reporter assay, rescue experiments and western blotting.

Results

We identified a novel circRNA circEAF2, which was downregulated in EBV + DLBCL and negatively correlated with EBV infection and DLBCL progression. In EBV-positive B lymphoma cells, circEAF2 overexpression induced lymphoma cell apoptosis and sensitized lymphoma cells to epirubicin. As mechanism of action, circEAF2 specifically targeted EBV-encoded miR-BART19-3p, upregulated APC, and suppressed downstream β-catenin expression, resulting in inactivation of Wnt signaling pathway and inhibition of EBV + DLBCL cell proliferation. In EBV-positive B-lymphoma murine models, xenografted tumors with circEAF2 overexpression presented decreased Ki-67 positivity, increased cell apoptosis and retarded tumor growth.

Conclusions

CircEAF2 counteracted EBV + DLBCL progression via miR-BART19-3p/APC/β-catenin axis, referring circEAF2 as a potential prognostic biomarker. Therapeutic targeting EBV-encoded miRNA may be a promising strategy in treating EBV-associated lymphoid malignancies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01458-9.

Epstein-Barr Virus miR-BART1-3p Regulates the miR-17-92 Cluster by Targeting E2F3

Epstein-Barr virus (EBV) is associated with several tumors and generates BamHI A rightward transcript (BART) microRNAs (miRNAs) from BART transcript introns. These BART miRNAs are expressed at higher levels in EBV-associated epithelial malignancies than in EBV-infected B lymphomas. To test the effects of EBV miRNA on the cell cycle and cell growth, we transfected miR-BART1-3p, a highly expressed EBV-associated miRNA, into gastric carcinoma cells. We found that miR-BART1-3p induced G0/G1 arrest and suppressed cell growth in gastric carcinoma cells. As our microarray analyses showed that E2F3, a cell cycle regulator, was inhibited by EBV infection, we hypothesized that miR-BART1-3p regulates E2F3. Luciferase assays revealed that miR-BART1-3p directly targeted the 3′-UTR of E2F3 mRNA. Both E2F3 mRNA and encoded protein levels were reduced following miR-BART1-3p transfection. In contrast, E2F3 expression in AGS-EBV cells transfected with a miR-BART1-3p inhibitor was enhanced. As E2F3 has been shown to regulate the expression of highly conserved miR-17-92 clusters in vertebrates, we examined whether this expression is affected by miR-BART1-3p, which can downregulate E2F3. The expression of E2F3, miR-17-92a-1 cluster host gene (MIR17HG), and miR-17-92 cluster miRNAs was significantly reduced in EBV-associated gastric carcinoma (EBVaGC) patients compared with EBV-negative gastric carcinoma (EBVnGC) patients. Further, miR-BART1-3p as well as the siRNA specific to E2F3 inhibited the expression of the miR-17-92 cluster, while inhibition of miR-BART1-3p enhanced the expression of the miR-17-92 cluster in cultured GC cells. Our results suggest a possible role of miR-BART1-3p in cell cycle regulation and in regulation of the miR-17-92 cluster through E2F3 suppression.

Inter-regulatory role of microRNAs in interaction between viruses and stem cells

MicroRNAs (miRNAs) are well known for post-transcriptional regulatory ability over specific mRNA targets. miRNAs exhibit temporal or tissue-specific expression patterns and regulate the cell and tissue developmental pathways. They also have determinative roles in production and differentiation of multiple lineages of stem cells and might have therapeutic advantages. miRNAs are a part of some viruses’ regulatory machinery, not a byproduct. The trace of miRNAs was detected in the genomes of viruses and regulation of cell reprograming and viral pathogenesis. Combination of inter-regulatory systems has been detected for miRNAs during viral infections in stem cells. Contraction between viruses and stem cells may be helpful in therapeutic tactics, pathogenesis, controlling viral infections and defining stem cell developmental strategies that is programmed by miRNAs as a tool. Therefore, in this review we intended to study the inter-regulatory role of miRNAs in the interaction between viruses and stem cells and tried to explain the advantages of miRNA regulatory potentials, which make a new landscape for future studies.

Hydrogel-based hybridization chain reaction (HCR) for detection of urinary exosomal miRNAs as a diagnostic tool of prostate cancer

Although urinary exosomal microRNAs (miRNAs) have recently emerged as potential biomarkers, clinical applications are still limited due to their low concentration in small volumes of clinical samples. Therefore, the development of a non-invasive, specific diagnostic tool, along with profiling exosomal miRNA markers from urine, remains a significant challenge. Here, we present hydrogel-based hybridization chain reaction (HCR) for multiplex signal amplification to detect urinary exosomal miRNAs from human clinical samples. We succeeded in identifying small amounts (~amol) of exosomal miRNAs from 600 μL of urine with up to ~35-fold amplification and enhanced detection limits by over an order of magnitude for two miRNA biomarker candidates, hsa-miR-6090 and hsa-miR-3665. Furthermore, we proposed ratiometric analysis without requiring normalization to a reference miRNA and validated the clinical diagnostic potential toward differentiating prostate cancer patients from healthy controls. Our hydrogel-based HCR could serve as a new diagnostic platform for a non-invasive liquid biopsy before burdensome tissue biopsy of various diseases, including prostate cancer screening, complementing the PSA test.

Identification and host response interaction study of SARS-CoV-2 encoded miRNA-like sequences: an in silico approach

COVID-19, a global pandemic caused by an RNA virus named SARS-CoV-2 has brought the world to a standstill in terms of infectivity, casualty, and commercial plummet. RNA viruses can encode microRNAs (miRNAs) capable of modulating host gene expression, and with that notion, we aimed to predict viral miRNA like sequences of MERS-CoV, SARS-CoV and SARS-CoV-2, analyze sequence reciprocity and investigate SARS-CoV-2 encoded potential miRNA-human genes interaction using bioinformatics tools. In this study, we retrieved 206 SARS-CoV-2 genomes, executed phylogenetic analysis, and the selected reference genome (MT434792.1) exhibited about 99% similarities among the retrieved genomes. We predicted 402, 137, and 85 putative miRNAs of MERS-CoV (NC_019843.3), SARS-CoV (NC_004718.3), and SARS-CoV-2 (MT434792.1) genome, respectively. Sequence similarity was analyzed among 624 miRNAs which revealed that the predicted miRNAs of SARS-CoV-2 share a cluster with the clad of miRNAs from MERS-CoV and SARS-CoV. Only SARS-CoV-2 derived 85 miRNAs were encountered for target prediction and 29 viral miRNAs seemed to target 119 human genes. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis suggested the involvement of respective genes in various pathways and biological processes. Finally, we focused on eight putative miRNAs influencing 14 genes that are involved in the adaptive hypoxic response, neuroinvasion and hormonal regulation, and tumorigenic progression in patients with COVID-19. SARS-CoV-2 encoded miRNAs may cause misexpression of some critical regulators and facilitate viral neuroinvasion, altered hormonal axis, and tumorigenic events in the human host. However, these propositions need validation from future studies.

miRNAs in SARS-CoV 2: A Spoke in the Wheel of Pathogenesis

INTRODUCTION

The rapid emergence of Severe Acute Respiratory Syndrome coronavirus 2 (SARS-- CoV-2) has resulted in an increased mortality rate across the globe. However, the underlying mechanism of SARS-CoV-2 altering human immune response is still elusive. The existing literature on miRNA mediated pathogenesis of RNA virus viz. Dengue virus, West Nile virus, etc. raises a suspicion that miRNA encoded by SARS-CoV-2 might facilitate virus replication and regulate the host's gene expression at the post-transcriptional level.

METHODS

We investigated this possibility via computational prediction of putative miRNAs encoded by the SARS-CoV-2 genome using a novel systematic pipeline that predicts putative mature-miRNA and their targeted genes transcripts. To trace down if viral-miRNAs targeted the genes critical to the immune pathway, we assessed whether mature miRNA transcripts exhibit effective hybridization with the 3'UTR region of human gene transcripts. Conversely, we also tried to study human miRNA-mediated viral gene regulation to get insight into the miRNA mediated offense and defense mechanism of virus and its host organisms in toto.

RESULTS

Our analysis led us to shortlist six putative miRNAs that target, majorly, genes related to cell proliferation/ differentiation/signaling, and senescence. Nonetheless, they also target immune-related genes that directly/ indirectly orchestrate immune pathways like TNF (Tumor Necrosis Factor) signaling and Chemokine signaling pathways putatively serving as the nucleus to cytokine storms.

CONCLUSION

Besides, these six miRNAs were found to be conserved so far across 80 complete genomes of SARS-CoV-2 (NCBI Virus, last assessed 12 April 2020) including Indian strains that are also targeted by 7 human miRNAs and can, therefore, be exploited to develop MicroRNA-Attenuated Vaccines.

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt’s Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi’s sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Burkitt’s Lymphoma and gastric cancer are both associated with EBV, a prolific DNA tumor virus that latently resides in nearly all human beings. Despite mostly restricting viral gene expression to noncoding RNAs, EBV has important influences on the fitness of infected tumor cells. Here, we show that the miRNA class of viral noncoding RNAs are a major viral contributor to remodeling the tumor cell regulatory machinery in patient tumor samples. First, an assessment of miRNA expression in clinical tumor samples showed that EBV miRNAs are expressed at unexpectedly high levels relative to cell miRNAs. Using a highly specific miRNA target identification approach, CLASH, we comprehensively identified both viral and cellular miRNA targets and the relative abundance of each miRNA-mRNA interaction. We also show that viral miRNAs bind to and alter the expression of their mRNA targets more effectively than their cellular miRNA counterparts. Pathway analysis of the most effectively targeted mRNAs revealed enrichment of immune signaling pathways and we show a corresponding inverse correlation between EBV miRNA expression and infiltrating immune cells in EBV positive primary tumors. Altogether, this study shows that EBV miRNAs are key regulators of the tumor cell phenotype and the immune cell microenvironment.

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt's Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi's sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Computational prediction of potential siRNA and human miRNA sequences to silence orf1ab associated genes for future therapeutics against SARS-CoV-2

The coronavirus disease 2019 (COVID-19) is an ongoing pandemic caused by an RNA virus termed as severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). SARS-CoV-2 possesses an almost 30kbp long genome. The genome contains open-reading frame 1ab (ORF1ab) gene, the largest one of SARS-CoV-2, encoding polyprotein PP1ab and PP1a responsible for viral transcription and replication. Several vaccines have already been approved by the respective authorities over the world to develop herd immunity among the population. In consonance with this effort, RNA interference (RNAi) technology holds the possibility to strengthen the fight against this virus. Here, we have implemented a computational approach to predict potential short interfering RNAs including small interfering RNAs (siRNAs) and microRNAs (miRNAs), which are presumed to be intrinsically active against SARS-CoV-2. In doing so, we have screened miRNA library and siRNA library targeting the ORF1ab gene. We predicted the potential miRNA and siRNA candidate molecules utilizing an array of bioinformatic tools. By extending the analysis, out of 24 potential pre-miRNA hairpins and 131 siRNAs, 12 human miRNA and 10 siRNA molecules were sorted as potential therapeutic agents against SARS-CoV-2 based on their GC content, melting temperature (Tm), heat capacity (Cp), hybridization and minimal free energy (MFE) of hybridization. This computational study is focused on lessening the extensive time and labor needed in conventional trial and error based wet lab methods and it has the potential to act as a decent base for future researchers to develop a successful RNAi therapeutic.

Computationally predicted SARS-COV-2 encoded microRNAs target NFKB, JAK/STAT and TGFB signaling pathways

Recently an outbreak that emerged in Wuhan, China in December 2019, spread to the whole world in a short time and killed >1,410,000 people. It was determined that a new type of beta coronavirus called severe acute respiratory disease coronavirus type 2 (SARS-CoV-2) was causative agent of this outbreak and the disease caused by the virus was named as coronavirus disease 19 (COVID19). Despite the information obtained from the viral genome structure, many aspects of the virus-host interactions during infection is still unknown. In this study we aimed to identify SARS-CoV-2 encoded microRNAs and their cellular targets. We applied a computational method to predict miRNAs encoded by SARS-CoV-2 along with their putative targets in humans. Targets of predicted miRNAs were clustered into groups based on their biological processes, molecular function, and cellular compartments using GO and PANTHER. By using KEGG pathway enrichment analysis top pathways were identified. Finally, we have constructed an integrative pathway network analysis with target genes. We identified 40 SARS-CoV-2 miRNAs and their regulated targets. Our analysis showed that targeted genes including NFKB1, NFKBIE, JAK1-2, STAT3-4, STAT5B, STAT6, SOCS1-6, IL2, IL8, IL10, IL17, TGFBR1-2, SMAD2-4, HDAC1-6 and JARID1A-C, JARID2 play important roles in NFKB, JAK/STAT and TGFB signaling pathways as well as cells' epigenetic regulation pathways. Our results may help to understand virus-host interaction and the role of viral miRNAs during SARS-CoV-2 infection. As there is no current drug and effective treatment available for COVID19, it may also help to develop new treatment strategies.

Azadirachta indica MicroRNAs: Genome-Wide Identification, Target Transcript Prediction, and Expression Analyses

MicroRNAs are short, endogenous, non-coding RNAs, liable for essential regulatory function. Numerous miRNAs have been identified and studied in plants with known genomic or small RNA resources. Despite the availability of genomic and transcriptomic resources, the miRNAs have not been reported in the medicinal tree Azadirachta indica (Neem) till date. Here for the first time, we report extensive identification of miRNAs and their possible targets in A. indica which might help to unravel their therapeutic potential. A comprehensive search of miRNAs in the A. indica genome by C-mii tool was performed. Overall, 123 miRNAs classified into 63 families and their stem-loop hairpin structures were predicted. The size of the A. indica (ain)-miRNAs ranged between 19 and 23 nt in length, and their corresponding ain-miRNA precursor sequence MFEI value averaged as -1.147 kcal/mol. The targets of ain-miRNAs were predicted in A. indica as well as Arabidopsis thaliana plant. The gene ontology (GO) annotation revealed the involvement of ain-miRNA targets in developmental processes, transport, stress, and metabolic processes including secondary metabolism. Stem-loop qRT-PCR was carried out for 25 randomly selected ain-miRNAs and differential expression patterns were observed in different A. indica tissues. Expression of miRNAs and its targets shows negative correlation in a dependent manner.

Critical review of Epstein-Barr virus microRNAs relation with EBV-associated gastric cancer

Epstein-Barr virus (EBV)-associated gastric cancer (EBVaGC) is regarded as the most prevalent malignant tumor triggered by EBV infection. In recent years, increasing attention has been considered to recognize more about the disease process's exact mechanisms. There is accumulating evidence that showing epigenetic modifications play critical roles in the EBVaGC pathogenesis. MicroRNAs (miRNAs), as critical epigenetic modulators, are single-strand short noncoding RNA (length ~ <200 bp), which regulate gene expression through binding to the 3'-untranslated region (3'-UTR) of target RNA transcripts and either degrade or repress their activities. In the latest research on EBV, it was found that this virus could encode miRNAs. Mechanistically, EBV-encoded miRNAs are involved in carcinogenesis and the progression of EBV-associated malignancies. Moreover, these miRNAs implicated in immune evasion, identification of pattern recognition receptors, regulation of lymphocyte activation and lethality, modulation of infected host cell antigen, maintain of EBV infection status, promotion of cell proliferation, invasion and migration, and reduction of apoptosis. As good news, not only has recent data demonstrated the crucial function of EBV-encoded miRNAs in the pathogenesis of EBVaGC, but it has also been revealed that aberrant expression of exosomal miRNAs in EBVaGC has made them biomarkers for detection of EBVaGC. Regarding these substantial characterizes, the critical role of EBV-encoded miRNAs has been a hot topic in research. In this review, we will focus on the multiple mechanisms involved in EBVaGC caused by EBV-encoded miRNAs and briefly discuss their potential application in the clinic as a diagnostic biomarker.

EBV miR-BARTs and human lncRNAs: Shifting the balance in competing endogenous RNA networks in EBV-associated gastric cancer

Non-coding RNAs (ncRNAs) contribute to the regulation of gene expression. By acting as competing endogenous RNA (ceRNA), long non-coding RNAs (lncRNAs) hijack microRNAs (miRNAs) and inhibit their ability to bind their coding targets. Viral miRNAs can compete with and target the same transcripts as human miRNAs, shifting the balance in networks associated with multiple cellular processes and diseases. Epstein-Barr virus (EBV) is an example of how a subset of viral coding RNA and non-coding RNAs can cause deregulation of human transcripts and contribute to the development of EBV-associated malignancies. EBV non-coding transforming genes include lncRNAs (i.e circular RNAs), and small ncRNAs (i.e. miRNAs). Among the latter, most ongoing research has focused on miR-BARTs whereas target many genes associated with apoptosis and epithelial-mesenchymal transition, in EBV-associated gastric cancer (GC). In this review, we propose to include the interactions between EBV ncRNAs human transcripts in the hypothesis known as "competitive viral and host RNAs". These interactions may shift the balance in biological pathways such as apoptosis and epithelial-mesenchymal transition in EBV-associated gastric cancer.

EBV-encoded miRNAs can sensitize nasopharyngeal carcinoma to chemotherapeutic drugs by targeting BRCA1

Nasopharyngeal carcinoma (NPC) is an Epstein‐Barr virus (EBV)‐associated epithelial malignancy. The high expression of BART‐miRNAs (miR‐BARTs) during latent EBV infection in NPC strongly supports their pathological importance in cancer progression. Recently, we found that several BART‐miRNAs work co‐operatively to modulate the DNA damage response (DDR) by reducing Ataxia‐telangiectasia‐mutated (ATM) activity. In this study, we further investigated the role of miR‐BARTs on DDR. The immunohistochemical study showed that the DNA repair gene, BRCA1, is consistently down‐regulated in primary NPCs. Using computer prediction programs and a series of reporter assays, we subsequently identified the negative regulatory role of BART2‐3p, BART12, BART17‐5p and BART19‐3p in BRCA1 expression. The ectopic expression of these four miR‐BARTs suppressed endogenous BRCA1 expression in EBV‐negative epithelial cell lines, whereas BRCA1 expression was enhanced by repressing endogenous miR‐BARTs activities in C666‐1 cells. More importantly, suppressing BRCA1 expression in nasopharyngeal epithelial cell lines using miR‐BART17‐5p and miR‐BART19‐3p mimics reduced the DNA repair capability and increased the cell sensitivity to the DNA‐damaging chemotherapeutic drugs, cisplatin and doxorubicin. Our findings suggest that miR‐BARTs play a novel role in DDR and may facilitate the development of effective NPC therapies.

Popular Computational Tools Used for miRNA Prediction and Their Future Development Prospects

MicroRNAs (miRNAs) are 19-24 nucleotide (nt)-long noncoding, single-stranded RNA molecules that play significant roles in regulating the gene expression, growth, and development of plants and animals. From the year that miRNAs were first discovered until the beginning of the twenty-first century, researchers used experimental methods such as cloning and sequencing to identify new miRNAs and their roles in the posttranscriptional regulation of protein synthesis. Later, in the early 2000s, informatics approaches to the discovery of new miRNAs began to be implemented. With increasing knowledge about miRNA, more efficient algorithms have been developed for computational miRNA prediction. The miRNA research community, hoping for greater coverage and faster results, has shifted from cumbersome and expensive traditional experimental approaches to computational approaches. These computational methods started with homology-based comparisons of known miRNAs with orthologs in the genomes of other species; this method could identify a known miRNA in new species. Second-generation sequencing and next-generation sequencing of mRNA at different developmental stages and in specific tissues, in combination with a better search and alignment algorithm, have accelerated the process of predicting novel miRNAs in a particular species. Using the accumulated annotated miRNA sequence information, researchers have been able to design ab initio algorithms for miRNA prediction independent of genome sequence knowledge. Here, the methods recently used for miRNA computational prediction are summarized and classified into the following four categories: homology-based, target-based, scoring-based, and machine-learning-based approaches. Finally, the future developmental directions of miRNA prediction methods are discussed.

Parallel analysis of miRNAs and mRNAs suggests distinct regulatory networks in Crassostrea gigas infected by Ostreid herpesvirus 1

Background

Since 2008, the aquaculture production of Crassostrea gigas was heavily affected by mass mortalities associated to Ostreid herpesvirus 1 (OsHV-1) microvariants worldwide. Transcriptomic studies revealed the major antiviral pathways of the oyster immune response while other findings suggested that also small non-coding RNAs (sncRNA) such as microRNAs might act as key regulators of the oyster response against OsHV-1. To explore the explicit connection between small non-coding and protein-coding transcripts, we performed paired whole transcriptome analysis of sncRNA and messenger RNA (mRNA) in six oysters selected for different intensities of OsHV-1 infection.

Results

The mRNA profiles of the naturally infected oysters were mostly governed by the transcriptional activity of OsHV-1, with several differentially expressed genes mapping to the interferon, toll, apoptosis, and pro-PO pathways. In contrast, miRNA profiles suggested more complex regulatory mechanisms, with 15 differentially expressed miRNAs (DE-miRNA) pointing to a possible modulation of the host response during OsHV-1 infection. We predicted 68 interactions between DE-miRNAs and oyster 3′-UTRs, but only few of them involved antiviral genes. The sncRNA reads assigned to OsHV-1 rather resembled mRNA degradation products, suggesting the absence of genuine viral miRNAs.

Conclusions

We provided data describing the miRNAome during OsHV-1 infection in C. gigas. This information can be used to understand the role of miRNAs in healthy and diseased oysters, to identify new targets for functional studies and, eventually to disentangle cause and effect relationships during viral infections in marine mollusks.

An in-silico approach to study the possible interactions of miRNA between human and SARS-CoV2

•

The rising epidemics of SARS-CoV2 has globally made a serious concern.

•

Thus, understanding the virus-host relationship remains a serious concern.

•

Here, we adopted the in-silico approach to picturized the similarities between the miRNAs of SARS-CoV2 genomes and human.

•

Further, the assessments and prediction of miRNAs helped us to analyze and understand the mechanism of pathogenesis.

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.

Fated for decay: RNA elements targeted by viral endonucleases

For over a decade, studies of messenger RNA regulation have revealed an unprecedented level of connectivity between the RNA pool and global gene expression. These connections are underpinned by a vast array of RNA elements that coordinate RNA-protein and RNA-RNA interactions, each directing mRNA fate from transcription to translation. Consequently, viruses have evolved an arsenal of strategies to target these RNA features and ultimately take control of the pathways they influence, and these strategies contribute to the global shutdown of the host gene expression machinery known as “Host Shutoff”. This takeover of the host cell is mechanistically orchestrated by a number of non-homologous virally encoded endoribonucleases. Recent large-scale screens estimate that over 70 % of the host transcriptome is decimated by the expression of these viral nucleases. While this takeover strategy seems extraordinarily well conserved, each viral endonuclease has evolved to target distinct mRNA elements. Herein, we will explore each of these RNA structures/sequence features that render messenger RNA susceptible or resistant to viral endonuclease cleavage. By further understanding these targeting and escape mechanisms we will continue to unravel untold depths of cellular RNA regulation that further underscores the integral relationship between RNA fate and the fate of the cell.

Genome-wide computational prediction of miRNAs in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) revealed target genes involved in pulmonary vasculature and antiviral innate immunity

The current outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China threatened humankind worldwide. The coronaviruses contains the largest RNA genome among all other known RNA viruses, therefore the disease etiology can be understood by analyzing the genome sequence of SARS-CoV-2. In this study, we used an ab-intio based computational tool VMir to scan the complete genome of SARS-CoV-2 to predict pre-miRNAs. The potential pre-miRNAs were identified by ViralMir and mature miRNAs were recognized by Mature Bayes. Additionally, predicted mature miRNAs were analysed against human genome by miRDB server to retrieve target genes. Besides that we also retrieved GO (Gene Ontology) terms for pathways, functions and cellular components. We predicted 26 mature miRNAs from genome of SARS-CoV-2 that targets human genes involved in pathways like EGF receptor signaling, apoptosis signaling, VEGF signaling, FGF receptor signaling. Gene enrichment tool analysis and substantial literature evidences suggests role of genes like BMPR2 and p53 in pulmonary vasculature and antiviral innate immunity respectively. Our findings may help research community to understand virus pathogenesis.

Epigenetic Regulation of Kaposi's Sarcoma-Associated Herpesvirus Latency

Kaposi’s sarcoma-associated herpesvirus (KSHV) is an oncogenic γ-herpesvirus that infects humans and exhibits a biphasic life cycle consisting of latent and lytic phases. Following entry into host cells, the KSHV genome undergoes circularization and chromatinization into an extrachromosomal episome ultimately leading to the establishment of latency. The KSHV episome is organized into distinct chromatin domains marked by variations in repressive or activating epigenetic modifications, including DNA methylation, histone methylation, and histone acetylation. Thus, the development of KSHV latency is believed to be governed by epigenetic regulation. In the past decade, interrogation of the KSHV epitome by genome-wide approaches has revealed a complex epigenetic mark landscape across KSHV genome and has uncovered the important regulatory roles of epigenetic modifications in governing the development of KSHV latency. Here, we highlight many of the findings regarding the role of DNA methylation, histone modification, post-translational modification (PTM) of chromatin remodeling proteins, the contribution of long non-coding RNAs (lncRNAs) in regulating KSHV latency development, and the role of higher-order episomal chromatin architecture in the maintenance of latency and the latent-to-lytic switch.

Epigenetic control in Kaposi sarcoma-associated herpesvirus infection and associated disease

Kaposi sarcoma-associated herpesvirus (KSHV) is the etiologic agent of several malignancies of endothelial and B-cell origin. The fact that latently infected tumor cells in these malignancies do not express classical viral oncogenes suggests that pathogenesis of KSHV-associated disease results from multistep processes that, in addition to constitutive viral gene expression, may require accumulation of cellular alterations. Heritable changes of the epigenome have emerged as an important co-factor that contributes to the pathogenesis of many non-viral cancers. Since KSHV encodes a number of factors that directly or indirectly manipulate host cell chromatin, it is an intriguing possibility that epigenetic reprogramming also contributes to the pathogenesis of KSHV-associated tumors. The fact that heritable histone modifications have also been shown to regulate viral gene expression programs in KSHV-infected tumor cells underlines the importance of epigenetic control during latency and tumorigenesis. We here review what is presently known about the role of epigenetic regulation of viral and host chromatin in KSHV infection and discuss how viral manipulation of these processes may contribute to the development of KSHV-associated disease.

EBV microRNA-BHRF1-2-5p targets the 3'UTR of immune checkpoint ligands PD-L1 and PD-L2

Epstein-Barr virus-positive (EBV+) diffuse large B-cell lymphomas (DLBCLs) express high levels of programmed death ligand 1 (PD-L1) and PD-L2. MicroRNA (miR) regulation is an important mechanism for the fine-tuning of gene expression via 3'-untranslated region (3'UTR) targeting, and we have previously demonstrated strong EBV miR expression in EBV+ DLBCL. Whereas the EBV latent membrane protein-1 (LMP1) is known to induce PD-L1/L2, a potential counterregulatory role of EBV miR in the fine-tuning of PD-L1/L2 expression remains to be established. To examine this, a novel in vitro model of EBV+ DLBCL was developed, using the viral strain EBV WIL, which unlike common laboratory strains retains intact noncoding regions where several EBV miRs reside. This enabled interrogation of the relationship among EBV latency genes, cell of origin (COO), PD-L1, PD-L2, and EBV miRs. The model successfully recapitulated the full spectrum of B-cell differentiation, with 4 discrete COO phases: early and late germinal center B cells (GCBs) and early and late activated B cells (ABCs). Interestingly, PD-L1/L2 levels increased markedly during transition from late GCB to early ABC phase, after LMP1 upregulation. EBV miR-BamHI fragment H rightward open reading frame 1 (BHRF1)-2-5p clustered apart from other EBV miRs, rising during late GCB phase. Bioinformatic prediction, together with functional validation, confirmed EBV miR-BHRF1-2-5p bound to PD-L1 and PD-L2 3'UTRs to reduce PD-L1/L2 surface protein expression. Results indicate a novel mechanism by which EBV miR-BHRF1-2-5p plays a context-dependent counterregulatory role to fine-tune the expression of the LMP1-driven amplification of these inhibitory checkpoint ligands. Further identification of immune checkpoint-targeting miRs may enable potential novel RNA-based therapies to emerge.

The Interplay Between Viral-Derived miRNAs and Host Immunity During Infection

MicroRNAs are short non-coding RNAs that play a crucial role in the regulation of gene expression during cellular processes. The host-encoded miRNAs are known to modulate the antiviral defense during viral infection. In the last decade, multiple DNA and RNA viruses have been shown to produce miRNAs known as viral miRNAs (v-miRNAs) so as to evade the host immune response. In this review, we highlight the origin and biogenesis of viral miRNAs during the viral lifecycle. We also explore the role of viral miRNAs in immune evasion and hence in maintaining chronic infection and disease. Finally, we offer insights into the underexplored role of viral miRNAs as potential targets for developing therapeutics for treating complex viral diseases.

6mer Seed Toxicity in Viral microRNAs

MicroRNAs (miRNAs) are short double-stranded noncoding RNAs (19-23 nucleotides) that regulate gene expression by suppressing mRNAs through RNA interference. Targeting is determined by the seed sequence (position 2-7/8) of the mature miRNA. A minimal G-rich seed of just six nucleotides is highly toxic to cells by targeting genes essential for cell survival. A screen of 215 miRNAs encoded by 17 human pathogenic viruses (v-miRNAs) now suggests that a number of v-miRNAs can kill cells through a G-rich 6mer sequence embedded in their seed. Specifically, we demonstrate that miR-K12-6-5p, an oncoviral mimic of the tumor suppressive miR-15/16 family encoded by human Kaposi sarcoma-associated herpes virus, harbors a noncanonical toxic 6mer seed (position 3-8) and that v-miRNAs are more likely than cellular miRNAs to utilize a noncanonical 6mer seed. Our data suggest that during evolution viruses evolved to use 6mer seed toxicity to kill cells.

•

Tumor suppressive miR-15/16-5p with a toxic 6mer seed targets survival genes

•

kshv-miR-K12-6-5p, a paralog of hsa-miR-15/16-5p carries an offset toxic 6mer seed

•

A screen of 215 viral miRNAs identifies miRNAs that contain a toxic 6mer seed

•

Many human viral miRNAs have the capacity to kill through 6mer seed toxicity

A Gammaherpesvirus MicroRNA Targets EWSR1 (Ewing Sarcoma Breakpoint Region 1) In Vivo To Promote Latent Infection of Germinal Center B Cells

Gammaherpesviruses, including the human pathogens Epstein-Barr virus (EBV) and Kaposi's sarcoma-associated herpesvirus (KSHV), directly contribute to the genesis of multiple types of malignancies. In vivo, these viruses infect B cells and manipulate B cell biology to establish lifelong infection. To accomplish this, gammaherpesviruses employ an array of gene products, including miRNAs, short noncoding RNAs that bind to and repress protein synthesis from specific target mRNAs. The in vivo relevance of repression of targets of gammaherpesvirus miRNAs remains highly elusive. Here, we identified a murine gammaherpesvirus miRNA as critical for in vivo infection and validated the host mRNA EWSR1 (Ewing sarcoma breakpoint region 1) as the predominant target for this miRNA. Using a novel technology, we demonstrated that repression of EWSR1 was essential for in vivo infection of the critical B cell reservoir. These findings provide the first in vivo demonstration of the significance of repression of a specific host mRNA by a gammaherpesvirus miRNA.

Gammaherpesviruses, including the human pathogens Epstein-Barr virus (EBV) and Kaposi’s sarcoma-associated herpesvirus (KSHV), directly contribute to the genesis of multiple types of malignancies, including B cell lymphomas. In vivo, these viruses infect B cells and manipulate B cell biology to establish lifelong latent infection. To accomplish this, gammaherpesviruses employ an array of gene products, including microRNAs (miRNAs). Although numerous host mRNA targets of gammaherpesvirus miRNAs have been identified, the in vivo relevance of repression of these targets remains elusive due to species restriction. Murine gammaherpesvirus 68 (MHV68) provides a robust virus-host system to dissect the in vivo function of conserved gammaherpesvirus genetic elements. We identified here MHV68 mghv-miR-M1-7-5p as critical for in vivo infection and then validated host EWSR1 (Ewing sarcoma breakpoint region 1) as the predominant target for this miRNA. Using novel, target-specific shRNA-expressing viruses, we determined that EWSR1 repression in vivo was essential for germinal center B cell infection. These findings provide the first in vivo demonstration of the biological significance of repression of a specific host mRNA by a gammaherpesvirus miRNA.

Short RNA regulators: the past, the present, the future, and implications for precision medicine and health disparities

We herein provide a brief review of the trajectory that the field of short RNA research followed in the last 25 years. We place emphasis on the unexpected discoveries and the ramifications of these discoveries for the field, as well as offer some thoughts about what the next 25 years may bring. Arguably, the uncovered dependence of different types of short RNAs on individual attributes such as a person's sex, population origin, race, and on tissue type, tissue state, and disease was most unexpected. This dependence has important ramifications in that it will provide a boost to our understanding of the molecular mechanisms of health disparities as well as pave the way for novel approaches to designing improved and personalized diagnostics and therapeutics.

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.

Roles of Non-coding RNAs During Herpesvirus Infection

Non-coding RNAs (ncRNAs) play essential roles in multiple aspects of the life cycles of herpesviruses and contribute to lifelong persistence of herpesviruses within their respective hosts. In this chapter, we discuss the types of ncRNAs produced by the different herpesvirus families during infection, some of the cellular ncRNAs manipulated by these viruses, and the overall contributions of ncRNAs to the viral life cycle, influence on the host environment, and pathogenesis.

AcMNPV-miR-3 is a miRNA encoded by Autographa californica nucleopolyhedrovirus and regulates the viral infection by targeting ac101

MicroRNAs (miRNAs), which are small noncoding RNAs found in plants, animals, and many viruses, regulate various biological processes. Our group has previously reported the first miRNA encoded by Autographa californica multiple Nucleopolyhedrovirus (AcMNPV), AcMNPV-miR-1, which regulates the expression of three viral genes. This study characterizes another miRNA encoded by AcMNPV, AcMNPV-miR-3. This miRNA is located on the opposite strand of the viral gene ac101 coding sequence in the AcMNPV genome, and it can be detected at 6 h post-infection and accumulated to a peak around 12 h post-infection in AcMNPV infected Sf9 cells. Five viral genes (ac101, ac23, ac25, ac86, and ac98) were verified to be regulated by AcMNPV-miR-3. Ac101 was markedly down-regulated by AcMNPV-miR-3 that may be via a siRNA-like cleavage mode. Administrating excessive AcMNPV-miR-3 resulted in decreased production of infectious budded virions (BV) and accelerated the formation of occlusion-derived virions (ODV). These results suggest that AcMNPV-miR-3 may play a regulatory role in BV and ODV production.

Towards Better Understanding of KSHV Life Cycle: from Transcription and Posttranscriptional Regulations to Pathogenesis

Kaposi’s sarcoma-associated herpesvirus (KSHV), also known as human herpesvirus-8 (HHV-8), is etiologically linked to the development of Kaposi’s sarcoma, primary effusion lymphoma, and multicentric Castleman’s disease. These malignancies often occur in immunosuppressed individuals, making KSHV infection-associated diseases an increasing global health concern with persistence of the AIDS epidemic. KSHV exhibits biphasic life cycles between latent and lytic infection and extensive transcriptional and posttranscriptional regulation of gene expression. As a member of the herpesvirus family, KSHV has evolved many strategies to evade the host immune response, which help the virus establish a successful lifelong infection. In this review, we summarize the current research status on the biology of latent and lytic viral infection, the regulation of viral life cycles and the related pathogenesis.

The Phylogeographic Diversity of EBV and Admixed Ancestry in the Americas⁻Another Model of Disrupted Human-Pathogen Co-Evolution

Epstein-Barr virus (EBV) is an etiological agent for gastric cancer with significant worldwide variations. Molecular characterizations of EBV have shown phylogeographical variations among healthy populations and in EBV-associated diseases, particularly the cosegregated BamHI-I fragment and XhoI restriction site of exon 1 of the LMP-1 gene. In the Americas, both cosegregated variants are present in EBV carriers, which aligns with the history of Asian and European human migration to this continent. Furthermore, novel recombinant variants have been found, reflecting the genetic makeup of this continent. However, in the case of EBV-associated gastric cancer (EBV-associated GC), the cosegregated European BamHI-“i” fragment and XhoI restriction site strain prevails. Thus, we propose that a disrupted coevolution between viral phylogeographical strains and mixed human ancestry in the Americas might explain the high prevalence of this particular gastric cancer subtype. This cosegregated region contains two relevant transcripts for EBV-associated GC, the BARF-1 and miR-BARTs. Thus, genome-wide association studies (GWAS) or targeted sequencing of both transcripts may be required to clarify their role as a potential source of this disrupted coevolution.

The roles of EBV-encoded microRNAs in EBV-associated tumors

Epstein-Barr virus (EBV) is believed to be a pathogen causing a number of human cancers, but the pathogenic mechanisms remain unclear. An increasing number of studies have indicated that EBV-encoded microRNAs (EBV miRNAs) are expressed in a latency type- and tumor type-dependent manner, playing important roles in the development and progression of EBV-associated tumors. By targeting one or more genes of the virus and the host, EBV miRNAs are responsible for the deregulation of a variety of viral and host cell biological processes, including viral replication, latency maintenance, immune evasion, cell apoptosis and metabolism, and tumor proliferation and metastasis. In addition, some EBV miRNAs can be used as excellent diagnostic, prognostic and treatment efficacy predictive biomarkers for EBV-associated tumors. More importantly, EBV miRNA-targeting therapeutics have emerged and have been developing rapidly, which may open a new era in the treatment of EBV-associated tumors in the near future.

Epstein-Barr Virus Infection of Cell Lines Derived from Diffuse Large B-Cell Lymphomas Alters MicroRNA Loading of the Ago2 Complex

Diffuse large B-cell lymphoma (DLBCL) is an aggressive lymphoid tumor which is occasionally Epstein-Barr virus (EBV) positive and is further subtyped as activated B-cell DLBCL (ABC-DLBCL) and germinal center B-cell DLBCL (GCB-DLBCL), which has implications for prognosis and treatment. We performed Ago2 RNA immunoprecipitation followed by high-throughput RNA sequencing (Ago2-RIP-seq) to capture functionally active microRNAs (miRNAs) in EBV-negative ABC-DLBCL and GCB-DLBCL cell lines and their EBV-infected counterparts. In parallel, total miRNA profiles of these cells were determined to capture the cellular miRNA profile for comparison with the functionally active profile. Selected miRNAs with differential abundances were validated using real-time quantitative PCR (RT-qPCR) and Northern blotting. We found 6 miRNAs with differential abundances (2 upregulated and 4 downregulated miRNAs) between EBV-negative and -positive ABC-DLBCL cells and 12 miRNAs with differential abundances (3 upregulated and 9 downregulated miRNAs) between EBV-negative and -positive GCB-DLBCL cells. Eight and twelve miRNAs were confirmed using RT-qPCR in ABC-DLBCL and GCB-DLBCL cells, respectively. Selected miRNAs were analyzed in additional type I/II versus type III EBV latency DLBCL cell lines. Furthermore, upregulation of miR-221-3p and downregulation of let7c-5p in ABC-DLBCL cells and upregulation of miR-363-3p and downregulation of miR-423-5p in GCB-DLBCL cells were verified using RIP-Northern blotting. Our comprehensive sequence analysis of the DLBCL miRNA profiles identified sets of deregulated miRNAs by Ago2-RIP-seq. Our Ago2-IP-seq miRNA profile could be considered an important data set for the detection of deregulated functionally active miRNAs in DLBCLs and could possibly lead to the identification of miRNAs as biomarkers for the classification of DLBCLs or even as targets for personalized targeted treatment.IMPORTANCE Diffuse large B-cell lymphoma (DLBCL) is a highly aggressive tumor of lymphoid origin which is occasionally Epstein-Barr virus (EBV) positive. MicroRNAs are found in most multicellular organisms and even in viruses such as EBV. They regulate the synthesis of proteins by binding to their cognate mRNA. MicroRNAs are tethered to their target mRNAs by "Argonaute" proteins. Here we compared the overall miRNA content of the Ago2 complex by differential loading to the overall content of miRNAs in two DLBCL cell lines and their EBV-converted counterparts. In all cell lines, the Ago2 load was different from the overall expression of miRNAs. In addition, the loading of the Ago2 complex was changed upon infection with EBV. This indicates that the virus not only changes the overall content of miRNAs but also influences the expression of proteins by affecting the Ago complexes.

Detection of Epstein-Barr virus (EBV)-encoded microRNAs in plasma of patients with nasopharyngeal carcinoma

BACKGROUND

Nasopharyngeal carcinoma (NPC) latently infected by Epstein-Barr virus (EBV) expresses 40 EBV BART microRNAs (miRNAs). Difference in diagnostic efficacy of these miRNAs on NPC detection was observed. Here, we performed a comprehensive evaluation on the efficacy of these miRNAs.

METHODS

Quantitative polymerase chain reaction was performed on plasma nucleic acid isolated from patients with NPC and noncancer donors.

RESULTS

For primary NPC, BART2-5P, BART6-3P, BART7-3P, BART7-5P, BART9-5P, BART11-3P, BART17-5P, and BART19-5P were significantly elevated. For recurrent NPC, plasma levels of BART2-3P, BART2-5P, BART5-3P, BART5-5P, BART6-3P, BART8-3P, BART9-5P, BART17-5P, BART19-3P, and BART20-3P were significantly increased. Area under curve (AUC) analysis showed that BART19-5P had the best performance to identify NPC which was serologically EBV DNA undetectable. For recurrent NPC, BART8-3P and BART10-3P had highest AUC value for identifying cancer in EBV DNA undetectable plasma.

CONCLUSION

Our data supported the use of circulating EBV miRNAs in NPC and recurrent NPC detection.