Search

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. Author summary: 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. [ABSTRACT FROM AUTHOR]

The count of the nucleotides in a cloned, short genomic sequence has become an important criterion to annotate such a sequence as a miRNA molecule. While the majority of human mature miRNA sequences consist of 22 nucleotides, there exists discrepancy in the characteristic lengths of the miRNA sequences. There is also a lack of systematic studies on such length distribution and on the biological factors that are related to or may affect this length. In this paper, we intend to fill this gap by investigating the sequence structure of human miRNA molecules using statistics tools. We demonstrate that the traditional discrete probability distributions do not model the length distribution of the human mature miRNAs well, and we obtain the statistical distribution model with a decent fit. We observe that the four nucleotide bases in a miRNA sequence are not randomly distributed, implying that possible structural patterns such as dinucleotide (trinucleotide or higher order) may exist. Furthermore, we study the relationships of this length distribution to multiple important factors such as evolutionary conservation, tumorigenesis, the length of precursor loop structures, and the number of predicted targets. The association between the miRNA sequence length and the distributions of target site counts in corresponding predicted genes is also presented. This study results in several novel findings worthy of further investigation that include: (1) rapid evolution introduces variation to the miRNA sequence length distribution; (2) miRNAs with extreme sequence lengths are unlikely to be cancer-related; and (3) the miRNA sequence length is positively correlated to the precursor length and the number of predicted target genes. [ABSTRACT FROM AUTHOR]

Recent studies have shown the contribution of miRNAs to cancer pathogenesis. Prostate cancer is the most commonly diagnosed cancer in men. Unlike other major types of cancer, no single gene has been identified as being mutated in the majority of prostate tumors. This implies that the expression profiling of genes, including the non-coding miRNAs, may substantially vary across individual cases of this cancer. The within-class variability makes it possible to reconstruct or infer disease-specific miRNA-mRNA correlation and regulatory modular networks using high-dimensional microarray data of prostate tumor samples. Furthermore, since miRNAs and tumor suppressor genes are usually tissue specific, miRNA-mRNA modules could potentially differ between primary prostate cancer (PPC) and metastatic prostate cancer (MPC). We herein performed an in silico analysis to explore the miRNA-mRNA correlation network modules in the two tumor subtypes. Our analysis identified 5 miRNA-mRNA module pairs (MPs) for PPC and MPC, respectively. Each MP includes one positiveconnection (correlation) module and one negative-connection (correlation) module. The number of miRNAs or mRNAs (genes) in each module varies from 2 to 8 or from 6 to 622. The modules discovered for PPC are more informative than those for MPC in terms of the implicated biological insights. In particular, one negative-connection module in PPC fits well with the popularly recognized miRNA-mediated post-transcriptional regulation theory. That is, the 3'UTR sequences of the involved mRNAs (∼620) are enriched with the target site motifs of the 7 modular miRNAs, has-miR-106b, -191, -19b, -92a, - 92b, -93, and -141. About 330 GO terms and KEGG pathways, including TGF-beta signaling pathway that maintains tissue homeostasis and plays a crucial role in the suppression of the proliferation of cancer cells, are over-represented (adj.p<0.05) in the modular gene list. These computationally identified modules provide remarkable biological evidence for the interference of miRNAs in the development of prostate cancers and warrant additional follow-up in independent laboratory studies. [ABSTRACT FROM AUTHOR]