GateView: A Multi-Omics Platform for Gene Feature Analysis of Virus Receptors within Human Normal Tissues and Tumors.

Viruses are obligate intracellular parasites that rely on cell surface receptor molecules to complete the first step of invading host cells. The experimental method for virus receptor screening is time-consuming, and receptor molecules have been identified for less than half of known viruses. This study collected known human viruses and their receptor molecules. Through bioinformatics analysis, common characteristics of virus receptor molecules (including sequence, expression, mutation, etc.) were obtained to study why these membrane proteins are more likely to become virus receptors. An in-depth analysis of the cataloged virus receptors revealed several noteworthy findings. Compared to other membrane proteins, human virus receptors generally exhibited higher expression levels and lower sequence conservation. These receptors were found in multiple tissues, with certain tissues and cell types displaying significantly higher expression levels. While most receptor molecules showed noticeable age-related variations in expression across different tissues, only a limited number of them exhibited gender-related differences in specific tissues. Interestingly, in contrast to normal tissues, virus receptors showed significant dysregulation in various types of tumors, particularly those associated with dsRNA and retrovirus receptors. Finally, GateView, a multi-omics platform, was established to analyze the gene features of virus receptors in human normal tissues and tumors. Serving as a valuable resource, it enables the exploration of common patterns among virus receptors and the investigation of virus tropism across different tissues, population preferences, virus pathogenicity, and oncolytic virus mechanisms.

A human-specific insertion promotes cell proliferation and migration by enhancing TBC1D8B expression.

Human-specific insertions play important roles in human phenotypes and diseases. Here we reported a 446-bp insertion (Insert-446) in intron 11 of the TBC1D8B gene, located on chromosome X, and traced its origin to a portion of intron 6 of the EBF1 gene on chromosome 5. Interestingly, Insert-446 was present in the human Neanderthal and Denisovans genomes, and was fixed in humans after human-chimpanzee divergence. We have demonstrated that Insert-446 acts as an enhancer through binding transcript factors that promotes a higher expression of human TBC1D8B gene as compared with orthologs in macaques. In addition, over-expression TBC1D8B promoted cell proliferation and migration through "a dual finger" catalytic mechanism (Arg538 and Gln573) in the TBC domain in vitro and knockdown of TBC1D8B attenuated tumorigenesis in vivo. Knockout of Insert-446 prevented cell proliferation and migration in cancer and normal cells. Our results reveal that the human-specific Insert-446 promotes cell proliferation and migration by upregulating the expression of TBC1D8B gene. These findings provide a significant insight into the effects of human-specific insertions on evolution.

The translational landscape of human vascular smooth muscle cells identifies novel short open reading frame-encoded peptide regulators for phenotype alteration.

AIMS: The plasticity of vascular smooth muscle cells (VSMCs) enables them to alter phenotypes under various physiological and pathological stimuli. The alteration of VSMC phenotype is a key step in vascular diseases, including atherosclerosis. Although the transcriptome shift during VSMC phenotype alteration has been intensively investigated, uncovering multiple key regulatory signalling pathways, the translatome dynamics in this cellular process, remain largely unknown. Here, we explored the genome-wide regulation at the translational level of human VSMCs during phenotype alteration. METHODS AND RESULTS: We generated nucleotide-resolution translatome and transcriptome data from human VSMCs undergoing phenotype alteration. Deep sequencing of ribosome-protected fragments (Ribo-seq) revealed alterations in protein synthesis independent of changes in messenger ribonucleicacid levels. Increased translational efficiency of many translational machinery components, including ribosomal proteins, eukaryotic translation elongation factors and initiation factors were observed during the phenotype alteration of VSMCs. In addition, hundreds of candidates for short open reading frame-encoded polypeptides (SEPs), a class of peptides containing 200 amino acids or less, were identified in a combined analysis of translatome and transcriptome data with a high positive rate in validating their coding capability. Three evolutionarily conserved SEPs were further detected endogenously by customized antibodies and suggested to participate in the pathogenesis of atherosclerosis by analysing the transcriptome and single cell RNA-seq data from patient atherosclerotic artery samples. Gain- and loss-of-function studies in human VSMCs and genetically engineered mice showed that these SEPs modulate the alteration of VSMC phenotype through different signalling pathways, including the mitogen-activated protein kinase pathway and p53 pathway. CONCLUSION: Our study indicates that an increase in the capacity of translation, which is attributable to an increased quantity of translational machinery components, mainly controls alterations of VSMC phenotype at the level of translational regulation. In addition, SEPs could function as important regulators in the phenotype alteration of human VSMCs.

tModBase: deciphering the landscape of tRNA modifications and their dynamic changes from epitranscriptome data.

tRNA molecules contain dense, abundant modifications that affect tRNA structure, stability, mRNA decoding and tsRNA formation. tRNA modifications and related enzymes are responsive to environmental cues and are associated with a range of physiological and pathological processes. However, there is a lack of resources that can be used to mine and analyse these dynamically changing tRNA modifications. In this study, we established tModBase (https://www.tmodbase.com/) for deciphering the landscape of tRNA modification profiles from epitranscriptome data. We analysed 103 datasets generated with second- and third-generation sequencing technologies and illustrated the misincorporation and termination signals of tRNA modification sites in ten species. We thus systematically demonstrate the modification profiles across different tissues/cell lines and summarize the characteristics of tRNA-associated human diseases. By integrating transcriptome data from 32 cancers, we developed novel tools for analysing the relationships between tRNA modifications and RNA modification enzymes, the expression of 1442 tRNA-derived small RNAs (tsRNAs), and 654 DNA variations. Our database will provide new insights into the features of tRNA modifications and the biological pathways in which they participate.

SARS-CoV-2 causes a significant stress response mediated by small RNAs in the blood of COVID-19 patients.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a serious impact on the world. In this study, small RNAs from the blood of COVID-19 patients with moderate or severe symptoms were extracted for high-throughput sequencing and analysis. Interestingly, the levels of a special group of tRNA-derived small RNAs (tsRNAs) were found to be dramatically upregulated after SARS-CoV-2 infection, particularly in coronavirus disease 2019 (COVID-19) patients with severe symptoms. In particular, the 3'CCA tsRNAs from tRNA-Gly were highly consistent with the inflammation indicator C-reactive protein (CRP). In addition, we found that the majority of significantly changed microRNAs (miRNAs) were associated with endoplasmic reticulum (ER)/unfolded protein response (UPR) sensors, which may lead to the induction of proinflammatory cytokine and immune responses. This study found that SARS-CoV-2 infection caused significant changes in the levels of stress-associated small RNAs in patient blood and their potential functions. Our research revealed that the cells of COVID-19 patients undergo tremendous stress and respond, which can be reflected or regulated by small non-coding RNA (sncRNAs), thus providing potential thought for therapeutic intervention in COVID-19 by modulating small RNA levels or activities.

tsRFun: a comprehensive platform for decoding human tsRNA expression, functions and prognostic value by high-throughput small RNA-Seq and CLIP-Seq data.

tRNA-derived small RNA (tsRNA), a novel type of regulatory small noncoding RNA, plays an important role in physiological and pathological processes. However, the understanding of the functional mechanism of tsRNAs in cells and their role in the occurrence and development of diseases is limited. Here, we integrated multiomics data such as transcriptome, epitranscriptome, and targetome data, and developed novel computer tools to establish tsRFun, a comprehensive platform to facilitate tsRNA research (http://rna.sysu.edu.cn/tsRFun/ or http://biomed.nscc-gz.cn/DB/tsRFun/). tsRFun evaluated tsRNA expression profiles and the prognostic value of tsRNAs across 32 types of cancers, identified tsRNA target molecules utilizing high-throughput CLASH/CLEAR or CLIP sequencing data, and constructed the interaction networks among tsRNAs, microRNAs, and mRNAs. In addition to its data presentation capabilities, tsRFun offers multiple real-time online tools for tsRNA identification, target prediction, and functional enrichment analysis. In summary, tsRFun provides a valuable data resource and multiple analysis tools for tsRNA investigation.

Genome-wide identification of microRNA targets reveals positive regulation of the Hippo pathway by miR-122 during liver development.

Liver development is a highly complex process that is regulated by the orchestrated interplay of epigenetic regulators, transcription factors, and microRNAs (miRNAs). Owing to the lack of global in vivo targets of all miRNAs during liver development, the mechanisms underlying the dynamic control of hepatocyte differentiation by miRNAs remain elusive. Here, using Argonaute (Ago) high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP) in the mouse liver at different developmental stages, we characterized massive Ago-binding RNAs and obtained a genome-wide map of liver miRNA-mRNA interactions. The dynamic changes of five clusters of miRNAs and their potential targets were identified to be differentially involved at specific stages, a dozen of high abundant miRNAs and their epigenetic regulation by super-enhancer were found during liver development. Remarkably, miR-122, a liver-specific and most abundant miRNA in newborn and adult livers, was found by its targetome and pathway reporter analyses to regulate the Hippo pathway, which is crucial for liver size control and homeostasis. Mechanistically, we further demonstrated that miR-122 negatively regulates the outcomes of the Hippo pathway transcription factor TEAD by directly targeting a number of hippo pathway regulators, including the coactivator TAZ and a key factor of the phosphatase complex PPP1CC, which contributes to the dephosphorylation of YAP, another coactivator downstream of the Hippo pathway. This study identifies for the first time the genome-wide miRNA targetomes during mouse liver development and demonstrates a novel mechanism of terminal differentiation of hepatocytes regulated by the miR-122/Hippo pathway in a coordinated manner. As the Hippo pathway plays important roles in cell proliferation and liver pathological processes like inflammation, fibrosis, and hepatocellular carcinoma (HCC), our study could also provide a new insight into the function of miR-122 in liver pathology.

PERK Signaling Controls Myoblast Differentiation by Regulating MicroRNA Networks.

The unfolded protein response (UPR) plays important roles in various cells that have a high demand for protein folding, which are involved in the process of cell differentiation and development. Here, we separately knocked down the three sensors of the UPR in myoblasts and found that PERK knockdown led to a marked transformation in myoblasts from a fusiform to a rounded morphology, which suggests that PERK is required for early myoblast differentiation. Interestingly, knocking down PERK induced reprogramming of C2C12 myoblasts into stem-like cells by altering the miRNA networks associated with differentiation and stemness maintenance, and the PERK-ATF4 signaling pathway transactivated muscle differentiation-associated miRNAs in the early stage of myoblast differentiation. Furthermore, we identified Ppp1cc as a direct target gene of miR-128 regulated by the PERK signaling pathway and showed that its repression is critical for a feedback loop that regulates the activity of UPR-associated signaling pathways, leading to cell migration, cell fusion, endoplasmic reticulum expansion, and myotube formation during myoblast differentiation. Subsequently, we found that the RNA-binding protein ARPP21, encoded by the host gene of miR-128-2, antagonized miR-128 activity by competing with it to bind to the 3' untranslated region (UTR) of Ppp1cc to maintain the balance of the differentiation state. Together, these results reveal the crucial role of PERK signaling in myoblast maintenance and differentiation and identify the mechanism underlying the role of UPR signaling as a major regulator of miRNA networks during early differentiation of myoblasts.

The cardiac translational landscape reveals that micropeptides are new players involved in cardiomyocyte hypertrophy.

Hypertrophic growth of cardiomyocytes is one of the major compensatory responses in the heart after physiological or pathological stimulation. Protein synthesis enhancement, which is mediated by the translation of messenger RNAs, is one of the main features of cardiomyocyte hypertrophy. Although the transcriptome shift caused by cardiac hypertrophy induced by different stimuli has been extensively investigated, translatome dynamics in this cellular process has been less studied. Here, we generated a nucleotide-resolution translatome as well as transcriptome data from isolated primary cardiomyocytes undergoing hypertrophy. More than 10,000 open reading frames (ORFs) were detected from the deep sequencing of ribosome-protected fragments (Ribo-seq), which orchestrated the shift of the translatome in hypertrophied cardiomyocytes. Our data suggest that rather than increase the translational rate of ribosomes, the increased efficiency of protein synthesis in cardiomyocyte hypertrophy was attributable to an increased quantity of ribosomes. In addition, more than 100 uncharacterized short ORFs (sORFs) were detected in long noncoding RNA genes from Ribo-seq with potential of micropeptide coding. In a random test of 15 candidates, the coding potential of 11 sORFs was experimentally supported. Three micropeptides were identified to regulate cardiomyocyte hypertrophy by modulating the activities of oxidative phosphorylation, the calcium signaling pathway, and the mitogen-activated protein kinase (MAPK) pathway. Our study provides a genome-wide overview of the translational controls behind cardiomyocyte hypertrophy and demonstrates an unrecognized role of micropeptides in cardiomyocyte biology.

ColorCells: a database of expression, classification and functions of lncRNAs in single cells.

Although long noncoding RNAs (lncRNAs) have significant tissue specificity, their expression and variability in single cells remain unclear. Here, we developed ColorCells (http://rna.sysu.edu.cn/colorcells/), a resource for comparative analysis of lncRNAs expression, classification and functions in single-cell RNA-Seq data. ColorCells was applied to 167 913 publicly available scRNA-Seq datasets from six species, and identified a batch of cell-specific lncRNAs. These lncRNAs show surprising levels of expression variability between different cell clusters, and has the comparable cell classification ability as known marker genes. Cell-specific lncRNAs have been identified and further validated by in vitro experiments. We found that lncRNAs are typically co-expressed with the mRNAs in the same cell cluster, which can be used to uncover lncRNAs' functions. Our study emphasizes the need to uncover lncRNAs in all cell types and shows the power of lncRNAs as novel marker genes at single cell resolution.

Novel organization of mitochondrial minicircles and guide RNAs in the zoonotic pathogen Trypanosoma lewisi.

Kinetoplastid flagellates are known for several unusual features, one of which is their complex mitochondrial genome, known as kinetoplast (k) DNA, composed of mutually catenated maxi- and minicircles. Trypanosoma lewisi is a member of the Stercorarian group of trypanosomes which is, based on human infections and experimental data, now considered a zoonotic pathogen. By assembling a total of 58 minicircle classes, which fall into two distinct categories, we describe a novel type of kDNA organization in T. lewisi. RNA-seq approaches allowed us to map the details of uridine insertion and deletion editing events upon the kDNA transcriptome. Moreover, sequencing of small RNA molecules enabled the identification of 169 unique guide (g) RNA genes, with two differently organized minicircle categories both encoding essential gRNAs. The unprecedented organization of minicircles and gRNAs in T. lewisi broadens our knowledge of the structure and expression of the mitochondrial genomes of these human and animal pathogens. Finally, a scenario describing the evolution of minicircles is presented.

Differential impacts of charcoal-stripped fetal bovine serum on c-Myc among distinct subtypes of breast cancer cell lines.

Charcoal-stripped fetal bovine serum (CS-FBS) is frequently used in studies on hormone-responsive cancers to provide hormone-free cell culture conditions. CS-FBS may influence the growth of cancer cells; however, the underlying mechanisms remain unclear. In this study, we aimed to clarify the effects of CS-FBS on distinct subtypes of breast cancer cells. We found that the crucial oncoprotein c-Myc was significantly inhibited in estrogen receptor alpha (ER-α)-positive breast cancer cells when cultured in CS-FBS-supplemented medium, but it was not suppressed in ER-α-negative cells. The addition of 17ß-estradiol (E2) to CS-FBS-supplemented medium rescued the CS-FBS-induced inhibition of c-Myc, while treatment with 5α-dihydrotestosterone (DHT) suppressed c-Myc expression. Our data demonstrated that CS-FBS may impede the growth of ER-α-positive breast cancer cells via c-Myc inhibition, and this was possibly due to the removal of estrogen. These results highlighted that the core drivers of c-Myc expression were subtype-specific depending on the distinct cell context and special caution should be exercised when using CS-FBS in studies of hormone-responsive cancer cells.

MicroRNA-122 supports robust innate immunity in hepatocytes by targeting the RTKs/STAT3 signaling pathway.

MicroRNA-122 (miR-122) is the most abundant microRNA in hepatocytes and a central player in liver biology and disease. Herein, we report a previously unknown role for miR-122 in hepatocyte intrinsic innate immunity. Restoration of miR-122 levels in hepatoma cells markedly enhanced the activation of interferons (IFNs) in response to a variety of viral nucleic acids or simulations, especially in response to hepatitis C virus RNA and poly (I:C). Mechanistically, miR-122 downregulated the phosphorylation (Tyr705) of STAT3, thereby removing the negative regulation of STAT3 on IFN-signaling. STAT3 represses IFN expression by inhibiting interferon regulatory factor 1 (IRF1), whereas miR-122 targets MERTK, FGFR1 and IGF1R, three receptor tyrosine kinases (RTKs) that directly promote STAT3 phosphorylation. This work identifies a miR-122-RTKs/STAT3-IRF1-IFNs regulatory circuitry, which may play a pivotal role in regulating hepatocyte innate immunity. These findings renewed our knowledge of miR-122's function and have important implications for the treatment of hepatitis viruses.

miR-372 and miR-373 enhance the stemness of colorectal cancer cells by repressing differentiation signaling pathways.

miR-372/373, a cluster of stem cell-specific microRNAs transactivated by the Wnt pathway, has been reported to be dysregulated in various cancers, particularly colorectal cancer (CRC); however, the unique role of these microRNAs in cancer remains to be discovered. In the present study, we characterized the upregulation in expression of miR-372/373 in CRC tissues from The Cancer Genome Atlas data, and then showed that overexpression of miR-372/373 enhanced the stemness of CRC cells by enriching the CD26/CD24-positive cell population and promoting self-renewal, chemotherapy resistance and the invasive potential of CRC cells. To clarify the mechanism underlying microRNA-induced stemness, we profiled 45 cell signaling pathways in CRC cells overexpressing miR-372/373 and found that stemness-related pathways, such as Nanog and Hedgehog, were upregulated. Instead, differentiation-related pathways, such as NFκB, MAPK/Erk and VDR, were markedly repressed by miR-372/373. Numerous new targets of miR-372/373 were identified, including SPOP, VDR and SETD7, all of which are factors important for cell differentiation. Furthermore, in contrast to the increase in miR-372/373 expression in CRC tissues, the expression levels of SPOP and VDR mRNA were significantly downregulated in these tissues, indicative of the poor differentiation status of CRC. Taken together, our findings suggest that miR-372/373 enhance CRC cell stemness by repressing the expression of differentiation genes. These results provide new insights for understanding the function and mechanisms of stem cell-specific microRNAs in the development of metastasis and drug resistance in CRC.

A group of tissue-specific microRNAs contribute to the silencing of CUX1 in different cell lineages during development.

Cut-like homeobox 1 (CUX1) is a highly conserved homeoprotein that functions as a transcriptional repressor of genes specifying terminal differentiation. We previously showed that liver-specific microRNA-122 (miR-122) regulates the timing of liver development by silencing CUX1 post-transcriptionally. Since the CUX1 protein is expressed in a subset of embryonic tissues, we hypothesized that it is regulated by specific microRNAs (miRNAs) in each cell type during development. Using a large-scale screening method, we identified ten tissue-specific miRNAs from different cell lineages that directly targeted CUX1. An analysis of the interaction between heart-specific microRNA-208a (miR-208a) and CUX1 in the hearts of developing mouse embryos and in P19CL6 cells undergoing cardiac differentiation indicated that CUX1 is regulated by miR-208a during heart development and cardiomyocyte differentiation. Functional analysis of miR-208a in P19CL6 cells using lentiviral-mediated over-expression showed that it regulates the transition between cellular proliferation and differentiation. These results suggest that these tissue-specific miRNAs might play a common role in timing the progression of terminal differentiation of different cell lineages, possibly by silencing the differentiation repressor CUX1.

Inhibition of the JNK/MAPK signaling pathway by myogenesis-associated miRNAs is required for skeletal muscle development.

Skeletal muscle differentiation is controlled by multiple cell signaling pathways, however, the JNK/MAPK signaling pathway dominating this process has not been fully elucidated. Here, we report that the JNK/MAPK pathway was significantly downregulated in the late stages of myogenesis, and in contrast to P38/MAPK pathway, it negatively regulated skeletal muscle differentiation. Based on the PAR-CLIP-seq analysis, we identified six elevated miRNAs (miR-1a-3p, miR-133a-3p, miR-133b-3p, miR-206-3p, miR-128-3p, miR-351-5p), namely myogenesis-associated miRNAs (mamiRs), negatively controlled the JNK/MAPK pathway by repressing multiple factors for the phosphorylation of the JNK/MAPK pathway, including MEKK1, MEKK2, MKK7, and c-Jun but not JNK protein itself, and as a result, expression of transcriptional factor MyoD and mamiRs were further promoted. Our study revealed a novel double-negative feedback regulatory pattern of cell-specific miRNAs by targeting phosphorylation kinase signaling cascade responsible for skeletal muscle development.

Comprehensive Genomic Characterization of RNA-Binding Proteins across Human Cancers.

RNA-binding proteins (RBPs) regulate the expression of thousands of transcripts, and some are reported to be involved in human tumorigenesis. However, little is known about the dysregulation of RBPs at the genomic level in human cancers. Here, we conducted comprehensive analyses for expression, somatic copy number alteration (SCNA), and mutation profiles of 1,542 RBPs in ∼7,000 clinical specimens across 15 cancer types. We identified markedly dysregulated RBPs and found that downregulation was a predominant pattern in cancer. Combined with recurrent SCNA data, we identified 76 RBPs as potential drivers. We also discovered a set of 139 RBPs that were significantly mutated in cancers. We confirmed the oncogenic property of six RBPs in colorectal and liver cancer cell lines by using in vitro functional experiments. Our study highlights the potential roles of RBPs in carcinogenesis and lays the groundwork to better understand the functions and mechanisms of RBPs in cancer.

RMBase v2.0: deciphering the map of RNA modifications from epitranscriptome sequencing data.

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 (m6A), ∼5400 N1-methyladenosines (m1A), ∼9600 pseudouridine (Ψ) modifications, ∼1000 5-methylcytosine (m5C) 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.

dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease.

Although thousands of pseudogenes have been annotated in the human genome, their transcriptional regulation, expression profiles and functional mechanisms are largely unknown. In this study, we developed dreamBase (http://rna.sysu.edu.cn/dreamBase) to facilitate the investigation of DNA modification, RNA regulation and protein binding of potential expressed pseudogenes from multidimensional high-throughput sequencing data. Based on ∼5500 ChIP-seq and DNase-seq datasets, we identified genome-wide binding profiles of various transcription-associated factors around pseudogene loci. By integrating ∼18 000 RNA-seq data, we analysed the expression profiles of pseudogenes and explored their co-expression patterns with their parent genes in 32 cancers and 31 normal tissues. By combining microRNA binding sites, we demonstrated complex post-transcriptional regulation networks involving 275 microRNAs and 1201 pseudogenes. We generated ceRNA networks to illustrate the crosstalk between pseudogenes and their parent genes through competitive binding of microRNAs. In addition, we studied transcriptome-wide interactions between RNA binding proteins (RBPs) and pseudogenes based on 458 CLIP-seq datasets. In conjunction with epitranscriptome sequencing data, we also mapped 1039 RNA modification sites onto 635 pseudogenes. This database will provide insights into the transcriptional regulation, expression, functions and mechanisms of pseudogenes as well as their roles in biological processes and diseases.

Exo-miRExplorer: A Comprehensive Resource for Exploring and Comparatively Analyzing Exogenous MicroRNAs.

MicroRNAs (miRNAs) are small regulatory RNAs that play important roles in animals, plants, and viruses. Deep-sequencing technology has been widely adopted in miRNA investigations. However, it is still a big mysterious why nearly all sequencing data contain miRNA sequences from exogenous species, called exo-miRNAs. In this study, we developed a novel platform, exo-miRExplorer, for mining and identifying exo-miRNAs from high-throughput small RNA sequencing experiments which originated from tissues and cell lines of multiple organisms. Thousands of exo-miRNAs are characterized with their expression abundance, the RNA families, original organisms and the sequencing platforms presented in exo-miRExplorer. Subsequently, we used exo-miRExplorer to perform further analysis. Comparative analysis of the exo-miRNAs between different sequencing datasets revealed significant correlation of exo-miRNAs between experiments in the same study. The plant-derived exo-miRNAs analysis provided robust evidence for non-diet source of exo-miRNAs. Virus-derived exo-miRNA analysis showed that pathogen RNAs could transfer to host cells and exist in deep-sequencing result at abundance level. In conclusion, exo-miRExplorer provides users with an integrative resource to facilitate detection and analysis of exo-miRNAs. exo-miRExplorer is available at the following URL: http://rna.sysu.edu.cn/exomiRDB/.