miRNA and lncRNA Expression Networks Modulate Cell Cycle and DNA Repair Inhibition in Senescent Prostate Cells

Cellular senescence is a state of permanent growth arrest that arises once cells reach the limit of their proliferative capacity. It creates an inflammatory microenvironment favouring the initiation and progression of various age-related diseases, including prostate cancer. Non-coding RNAs (ncRNAs) have emerged as important regulators of cellular gene expression. Nonetheless, very little is known about the interplay of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) and how deregulation of ncRNA networks promotes cellular senescence. To investigate this, human prostate epithelial cells were cultured through different passages until senescent, and their RNA was extracted and sequenced using RNA sequencing (RNAseq) and microRNA sequencing (miRNA-seq) miRNAseq. Differential expression (DE) gene analysis was performed to compare senescent and proliferating cells with Limma, miRNA-target interactions with multiMiR, lncRNA-target interactions using TCGA data and network evaluation with miRmapper. We found that miR-335-3p, miR-543 and the lncRNAs H19 and SMIM10L2A all play central roles in the regulation of cell cycle and DNA repair processes. Expression of most genes belonging to these pathways were down-regulated by senescence. Using the concept of network centrality, we determined the top 10 miRNAs and lncRNAs, with miR-335-3p and H19 identified as the biggest hubs for miRNAs and lncRNA respectively. These ncRNAs regulate key genes belonging to pathways involved in cell senescence and prostate cancer demonstrating their central role in these processes and opening the possibility for their use as biomarkers or therapeutic targets to mitigate against prostate ageing and carcinogenesis

miRmapper: A Tool for Interpretation of miRNA–mRNA Interaction Networks

It is estimated that 30% of all genes in the mammalian cells are regulated by microRNA (miRNAs). The most relevant miRNAs in a cellular context are not necessarily those with the greatest change in expression levels between healthy and diseased tissue. Differentially expressed (DE) miRNAs that modulate a large number of messenger RNA (mRNA) transcripts ultimately have a greater influence in determining phenotypic outcomes and are more important in a global biological context than miRNAs that modulate just a few mRNA transcripts. Here, we describe the development of a tool, “miRmapper”, which identifies the most dominant miRNAs in a miRNA–mRNA network and recognizes similarities between miRNAs based on commonly regulated mRNAs. Using a list of miRNA–target gene interactions and a list of DE transcripts, miRmapper provides several outputs: (1) an adjacency matrix that is used to calculate miRNA similarity utilizing the Jaccard distance; (2) a dendrogram and (3) an identity heatmap displaying miRNA clusters based on their effect on mRNA expression; (4) a miRNA impact table and (5) a barplot that provides a visual illustration of this impact. We tested this tool using nonmetastatic and metastatic bladder cancer cell lines and demonstrated that the most relevant miRNAs in a cellular context are not necessarily those with the greatest fold change. Additionally, by exploiting the Jaccard distance, we unraveled novel cooperative interactions between miRNAs from independent families in regulating common target mRNAs; i.e., five of the top 10 miRNAs act in synergy

miRNA and lncRNA Expression Networks Modulate Cell Cycle and DNA Repair Inhibition in Senescent Prostate Cells

Cellular senescence is a state of permanent growth arrest that arises once cells reach the limit of their proliferative capacity. It creates an inflammatory microenvironment favouring the initiation and progression of various age-related diseases, including prostate cancer. Non-coding RNAs (ncRNAs) have emerged as important regulators of cellular gene expression. Nonetheless, very little is known about the interplay of microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) and how deregulation of ncRNA networks promotes cellular senescence. To investigate this, human prostate epithelial cells were cultured through different passages until senescent, and their RNA was extracted and sequenced using RNA sequencing (RNAseq) and microRNA sequencing (miRNA-seq) miRNAseq. Differential expression (DE) gene analysis was performed to compare senescent and proliferating cells with Limma, miRNA-target interactions with multiMiR, lncRNA-target interactions using TCGA data and network evaluation with miRmapper. We found that miR-335-3p, miR-543 and the lncRNAs H19 and SMIM10L2A all play central roles in the regulation of cell cycle and DNA repair processes. Expression of most genes belonging to these pathways were down-regulated by senescence. Using the concept of network centrality, we determined the top 10 miRNAs and lncRNAs, with miR-335-3p and H19 identified as the biggest hubs for miRNAs and lncRNA respectively. These ncRNAs regulate key genes belonging to pathways involved in cell senescence and prostate cancer demonstrating their central role in these processes and opening the possibility for their use as biomarkers or therapeutic targets to mitigate against prostate ageing and carcinogenesis

GAIL: An interactive webserver for inference and dynamic visualization of gene-gene associations based on gene ontology guided mining of biomedical literature

In systems biology, inference of functional associations among genes is compelling because the construction of functional association networks facilitates biomarker discovery. Specifically, such gene associations in human can help identify putative biomarkers that can be used as diagnostic tools in treating patients. Although biomedical literature is considered a valuable data source for this task, currently only a limited number of webservers are available for mining gene-gene associations from the vast amount of biomedical literature using text mining techniques. Moreover, these webservers often have limited coverage of biomedical literature and also lack efficient and user-friendly tools to interpret and visualize mined relationships among genes. To address these limitations, we developed GAIL (Gene-gene Association Inference based on biomedical Literature), an interactive webserver that infers human gene-gene associations from Gene Ontology (GO) guided biomedical literature mining and provides dynamic visualization of the resulting association networks and various gene set enrichment analysis tools. We evaluate the utility and performance of GAIL with applications to gene signatures associated with systemic lupus erythematosus and breast cancer. Results show that GAIL allows effective interrogation and visualization of gene-gene networks and their subnetworks, which facilitates biological understanding of gene-gene associations. GAIL is available at http://chunglab.io/GAIL/

Time- and Oil-Dependent Transcriptomic and Physiological Responses to Deepwater Horizon Oil in Mahi-Mahi (Coryphaena hippurus) Embryos and Larvae

The Deepwater Horizon (DWH) oil spill contaminated the spawning habitats for numerous commercially and ecologically important fishes. Exposure to the water accommodated fraction (WAF) of oil from the spill has been shown to cause cardiac toxicity during early developmental stages across fishes. To better understand the molecular events and explore new pathways responsible for toxicity, RNA sequencing was performed in conjunction with physiological and morphological assessments to analyze the time-course (24, 48, and 96 h post fertilization (hpf)) of transcriptional and developmental responses in embryos/larvae of mahi-mahi exposed to WAF of weathered (slick) and source DWH oils. Slick oil exposure induced more pronounced changes in gene expression over time than source oil exposure. Predominant transcriptomic responses included alteration of EIF2 signaling, steroid biosynthesis, ribosome biogenesis and activation of the cytochrome P450 pathway. At 96 hpf, slick oil exposure resulted in significant perturbations in eye development and peripheral nervous system, suggesting novel targets in addition to the heart may be involved in the developmental toxicity of DHW oil. Comparisons of changes of cardiac genes with phenotypic responses were consistent with reduced heart rate and increased pericardial edema in larvae exposed to slick oil but not source oil