Although microRNAs (miRNAs), other non-coding RNAs (ncRNAs) (e.g. lncRNAs, pseudogenes and circRNAs) and competing endogenous RNAs (ceRNAs) have been implicated in cell-fate determination and in various human diseases, surprisingly little is known about the regulatory interaction networks among the multiple classes of RNAs. In this study, we developed starBase v2.0 (http://starbase.sysu.edu.cn/) to systematically identify the RNA–RNA and protein–RNA interaction networks from 108 CLIP-Seq (PAR-CLIP, HITS-CLIP, iCLIP, CLASH) data sets generated by 37 independent studies. By analyzing millions of RNA-binding protein binding sites, we identified ∼9000 miRNA-circRNA, 16 000 miRNA-pseudogene and 285 000 protein–RNA regulatory relationships. Moreover, starBase v2.0 has been updated to provide the most comprehensive CLIP-Seq experimentally supported miRNA-mRNA and miRNA-lncRNA interaction networks to date. We identified ∼10 000 ceRNA pairs from CLIP-supported miRNA target sites. By combining 13 functional genomic annotations, we developed miRFunction and ceRNAFunction web servers to predict the function of miRNAs and other ncRNAs from the miRNA-mediated regulatory networks. Finally, we developed interactive web implementations to provide visualization, analysis and downloading of the aforementioned large-scale data sets. This study will greatly expand our understanding of ncRNA functions and their coordinated regulatory networks.

N⁶-methyladenosine (m⁶A) is the most prevalent modification in eukaryotic messenger RNAs (mRNAs) and is interpreted by its readers, such as YTH domain-containing proteins, to regulate mRNA fate. Here, we report the insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs; including IGF2BP1/2/3) as a distinct family of m⁶A readers that target thousands of mRNA transcripts through recognizing the consensus GG(m⁶A)C sequence. In contrast to the mRNA-decay-promoting function of YTH domain-containing family protein 2, IGF2BPs promote the stability and storage of their target mRNAs (for example, MYC) in an m⁶A-dependent manner under normal and stress conditions and therefore affect gene expression output. Moreover, the K homology domains of IGF2BPs are required for their recognition of m⁶A and are critical for their oncogenic functions. Thus, our work reveals a different facet of the m⁶A-reading process that promotes mRNA stability and translation, and highlights the functional importance of IGF2BPs as m⁶A readers in post-transcriptional gene regulation and cancer biology.

The abnormal transcriptional regulation of non-coding RNAs (ncRNAs) and protein-coding genes (PCGs) is contributed to various biological processes and linked with human diseases, but the underlying mechanisms remain elusive. In this study, we developed ChIPBase v2.0 (http://rna.sysu.edu.cn/chipbase/) to explore the transcriptional regulatory networks of ncRNAs and PCGs. ChIPBase v2.0 has been expanded with ∼10 200 curated ChIP-seq datasets, which represent about 20 times expansion when comparing to the previous released version. We identified thousands of binding motif matrices and their binding sites from ChIP-seq data of DNA-binding proteins and predicted millions of transcriptional regulatory relationships between transcription factors (TFs) and genes. We constructed ‘Regulator’ module to predict hundreds of TFs and histone modifications that were involved in or affected transcription of ncRNAs and PCGs. Moreover, we built a web-based tool, Co-Expression, to explore the co-expression patterns between DNA-binding proteins and various types of genes by integrating the gene expression profiles of ∼10 000 tumor samples and ∼9100 normal tissues and cell lines. ChIPBase also provides a ChIP-Function tool and a genome browser to predict functions of diverse genes and visualize various ChIP-seq data. This study will greatly expand our understanding of the transcriptional regulations of ncRNAs and PCGs.

More than 100 distinct chemical modifications to RNA have been characterized so far. However, the prevalence, mechanisms and functions of various RNA modifications remain largely unknown. To provide transcriptome-wide landscapes of RNA modifications, we developed the RMBase v2.0 (http://rna.sysu.edu.cn/rmbase/), which is a comprehensive database that integrates epitranscriptome sequencing data for the exploration of post-transcriptional modifications of RNAs and their relationships with miRNA binding events, disease-related single-nucleotide polymorphisms (SNPs) and RNA-binding proteins (RBPs). RMBase v2.0 was expanded with ∼600 datasets and ∼1 397 000 modification sites from 47 studies among 13 species, which represents an approximately 10-fold expansion when compared with the previous release. It contains ∼1 373 000 N6-methyladenosines (m⁶A), ∼5400 N1-methyladenosines (m¹A), ∼9600 pseudouridine (Ψ) modifications, ∼1000 5-methylcytosine (m⁵C) modifications, ∼5100 2′-O-methylations (2′-O-Me), and ∼2800 modifications of other modification types. Moreover, we built a new module called ‘Motif’ that provides the visualized logos and position weight matrices (PWMs) of the modification motifs. We also constructed a novel module termed ‘modRBP’ to study the relationships between RNA modifications and RBPs. Additionally, we developed a novel web-based tool named ‘modMetagene’ to plot the metagenes of RNA modification along a transcript model. This database will help researchers investigate the potential functions and mechanisms of RNA modifications.