Inferring plant microRNA functional similarity using a weighted protein-protein interaction network

A microRNA-microRNA crosstalk network inferred from genome-wide single nucleotide polymorphism variants in natural populations of Arabidopsis thaliana

To adapt to variable natural conditions, plants have evolved several strategies to respond to different environmental stresses. MicroRNA (miRNA)-mediated gene regulation is one of such strategies. Variants, e.g., single nucleotide polymorphisms (SNPs) within the mature miRNAs or their target sites may cause the alteration of regulatory networks and serious phenotype changes. In this study, we proposed a novel approach to construct a miRNA-miRNA crosstalk network in Arabidopsis thaliana based on the notion that two cooperative miRNAs toward common targets are under a strong pressure to be inherited together across ecotypes. By performing a genome-wide scan of the SNPs within the mature miRNAs and their target sites, we defined a "regulation fate profile" to describe a miRNA-target regulation being static (kept) or dynamic (gained or lost) across 1,135 ecotypes compared with the reference genome of Col-0. The cooperative miRNA pairs were identified by estimating the similarity of their regulation fate profiles toward the common targets. The reliability of the cooperative miRNA pairs was supported by solid expressional correlation, high PPImiRFS scores, and similar stress responses. Different combinations of static and dynamic miRNA-target regulations account for the cooperative miRNA pairs acting on various biological characteristics of miRNA conservation, expression, homology, and stress response. Interestingly, the targets that are co-regulated dynamically by both cooperative miRNAs are more likely to be responsive to stress. Hence, stress-related genes probably bear selective pressures in a certain group of ecotypes, in which miRNA regulations on the stress genes reprogram. Finally, three case studies showed that reprogramming miRNA-miRNA crosstalk toward the targets in specific ecotypes was associated with these ecotypes' climatic variables and geographical locations. Our study highlights the potential of miRNA-miRNA crosstalk as a genetic basis underlying environmental adaptation in natural populations.

3D RNA-seq: a powerful and flexible tool for rapid and accurate differential expression and alternative splicing analysis of RNA-seq data for biologists

RNA-sequencing (RNA-seq) analysis of gene expression and alternative splicing should be routine and robust but is often a bottleneck for biologists because of different and complex analysis programs and reliance on specialized bioinformatics skills. We have developed the '3D RNA-seq' App, an R shiny App and web-based pipeline for the comprehensive analysis of RNA-seq data from any organism. It represents an easy-to-use, flexible and powerful tool for analysis of both gene and transcript-level gene expression to identify differential gene/transcript expression, differential alternative splicing and differential transcript usage (3D) as well as isoform switching from RNA-seq data. 3D RNA-seq integrates state-of-the-art differential expression analysis tools and adopts best practice for RNA-seq analysis. The program is designed to be run by biologists with minimal bioinformatics experience (or by bioinformaticians) allowing lab scientists to analyse their RNA-seq data. It achieves this by operating through a user-friendly graphical interface which automates the data flow through the programs in the pipeline. The comprehensive analysis performed by 3D RNA-seq is extremely rapid and accurate, can handle complex experimental designs, allows user setting of statistical parameters, visualizes the results through graphics and tables, and generates publication quality figures such as heat-maps, expression profiles and GO enrichment plots. The utility of 3D RNA-seq is illustrated by analysis of data from a time-series of cold-treated Arabidopsis plants and from dexamethasone-treated male and female mouse cortex and hypothalamus data identifying dexamethasone-induced sex- and brain region-specific differential gene expression and alternative splicing.

Human microRNA similarity in breast cancer

MicroRNAs (miRNAs) play important roles in a variety of human diseases, including breast cancer. A number of miRNAs are up- and down-regulated in breast cancer. However, little is known about miRNA similarity and similarity network in breast cancer. Here, a collection of 272 breast cancer-associated miRNA precursors (pre-miRNAs) were utilized to calculate similarities of sequences, target genes, pathways and functions and construct a combined similarity network. Well-characterized miRNAs and their similarity network were highlighted. Interestingly, miRNA sequence-dependent similarity networks were not identified in spite of sequence–target gene association. Similarity networks with minimum and maximum number of miRNAs originate from pathway and mature sequence, respectively. The breast cancer-associated miRNAs were divided into seven functional classes (classes I–VII) followed by disease enrichment analysis and novel miRNA-based disease similarities were found. The finding would provide insight into miRNA similarity, similarity network and disease heterogeneity in breast cancer.

Genome-wide analysis of long non-coding RNAs responsive to multiple nutrient stresses in Arabidopsis thaliana

Nutrient stress is the most important environmental stress that limits plant growth and development. Although recent evidence highlights the vital functions of long non-coding RNAs (lncRNA) in response to single nutrient stress in some model plants, a comprehensive investigation of the effect of lncRNAs in response to nutrient stress has not been performed in Arabidopsis thaliana. Here, we presented the identification and characterization of lncRNAs under seven nutrient stress conditions. The expression pattern analysis revealed that aberrant expression of lncRNAs is a stress-specific manner under nutrient stress conditions and that lncRNAs are more sensitive to nutrient stress than protein-coding genes (PCGs). Moreover, competing endogenous RNA (ceRNA) network and lncRNA-mRNA co-expression network (CEN) were constructed to explore the potential function of these lncRNAs under nutrient stress conditions. We further combined different expressed lncRNAs with ceRNA network and CEN to select key lncRNAs in response to nutrient stress. Together, our study provides important information for further insights into the role of lncRNAs in response to stress in plants.

miRNA-miRNA crosstalk: from genomics to phenomics

The discovery of microRNA (miRNA)-miRNA crosstalk has greatly improved our understanding of complex gene regulatory networks in normal and disease-specific physiological conditions. Numerous approaches have been proposed for modeling miRNA-miRNA networks based on genomic sequences, miRNA-mRNA regulation, functional information and phenomics alone, or by integrating heterogeneous data. In addition, it is expected that miRNA-miRNA crosstalk can be reprogrammed in different tissues or specific diseases. Thus, transcriptome data have also been integrated to construct context-specific miRNA-miRNA networks. In this review, we summarize the state-of-the-art miRNA-miRNA network modeling methods, which range from genomics to phenomics, where we focus on the need to integrate heterogeneous types of omics data. Finally, we suggest future directions for studies of crosstalk of noncoding RNAs. This comprehensive summarization and discussion elucidated in this work provide constructive insights into miRNA-miRNA crosstalk.

Refine gene functional similarity network based on interaction networks

BACKGROUND

In recent years, biological interaction networks have become the basis of some essential study and achieved success in many applications. Some typical networks such as protein-protein interaction networks have already been investigated systematically. However, little work has been available for the construction of gene functional similarity networks so far. In this research, we will try to build a high reliable gene functional similarity network to promote its further application.

RESULTS

Here, we propose a novel method to construct and refine the gene functional similarity network. It mainly contains three steps. First, we establish an integrated gene functional similarity networks based on different functional similarity calculation methods. Then, we construct a referenced gene-gene association network based on the protein-protein interaction networks. At last, we refine the spurious edges in the integrated gene functional similarity network with the help of the referenced gene-gene association network. Experiment results indicate that the refined gene functional similarity network (RGFSN) exhibits a scale-free, small world and modular architecture, with its degrees fit best to power law distribution. In addition, we conduct protein complex prediction experiment for human based on RGFSN and achieve an outstanding result, which implies it has high reliability and wide application significance.

CONCLUSIONS

Our efforts are insightful for constructing and refining gene functional similarity networks, which can be applied to build other high quality biological networks.

Gene Ontology semantic similarity tools: survey on features and challenges for biological knowledge discovery

Gene Ontology (GO) semantic similarity tools enable retrieval of semantic similarity scores, which incorporate biological knowledge embedded in the GO structure for comparing or classifying different proteins or list of proteins based on their GO annotations. This facilitates a better understanding of biological phenomena underlying the corresponding experiment and enables the identification of processes pertinent to different biological conditions. Currently, about 14 tools are available, which may play an important role in improving protein analyses at the functional level using different GO semantic similarity measures. Here we survey these tools to provide a comprehensive view of the challenges and advances made in this area to avoid redundant effort in developing features that already exist, or implementing ideas already proven to be obsolete in the context of GO. This helps researchers, tool developers, as well as end users, understand the underlying semantic similarity measures implemented through knowledge of pertinent features of, and issues related to, a particular tool. This should empower users to make appropriate choices for their biological applications and ensure effective knowledge discovery based on GO annotations.

Significance and Functional Similarity for Identification of Disease Genes

One of the most significant research issues in functional genomics is insilico identification of disease related genes. In this regard, the paper presents a new gene selection algorithm, termed as SiFS, for identification of disease genes. It integrates the information obtained from interaction network of proteins and gene expression profiles. The proposed SiFS algorithm culls out a subset of genes from microarray data as disease genes by maximizing both significance and functional similarity of the selected gene subset. Based on the gene expression profiles, the significance of a gene with respect to another gene is computed using mutual information. On the other hand, a new measure of similarity is introduced to compute the functional similarity between two genes. Information derived from the protein-protein interaction network forms the basis of the proposed SiFS algorithm. The performance of the proposed gene selection algorithm and new similarity measure, is compared with that of other related methods and similarity measures, using several cancer microarray data sets.

A path-based measurement for human miRNA functional similarities using miRNA-disease associations

Compared with the sequence and expression similarity, miRNA functional similarity is so important for biology researches and many applications such as miRNA clustering, miRNA function prediction, miRNA synergism identification and disease miRNA prioritization. However, the existing methods always utilized the predicted miRNA target which has high false positive and false negative to calculate the miRNA functional similarity. Meanwhile, it is difficult to achieve high reliability of miRNA functional similarity with miRNA-disease associations. Therefore, it is increasingly needed to improve the measurement of miRNA functional similarity. In this study, we develop a novel path-based calculation method of miRNA functional similarity based on miRNA-disease associations, called MFSP. Compared with other methods, our method obtains higher average functional similarity of intra-family and intra-cluster selected groups. Meanwhile, the lower average functional similarity of inter-family and inter-cluster miRNA pair is obtained. In addition, the smaller p-value is achieved, while applying Wilcoxon rank-sum test and Kruskal-Wallis test to different miRNA groups. The relationship between miRNA functional similarity and other information sources is exhibited. Furthermore, the constructed miRNA functional network based on MFSP is a scale-free and small-world network. Moreover, the higher AUC for miRNA-disease prediction indicates the ability of MFSP uncovering miRNA functional similarity.