BackgroundAlthough microRNAs (miRNAs) are involved in all hallmarks of cancer, miRNA dysregulation in the metastatic process remains poorly understood. We investigated the role of miRNAs in metastatic evolution of colorectal cancer (CRC) by analyzing smallRNA-seq datasets from primary CRC, metastatic locations (liver, lung and peritoneum), and corresponding adjacent tissues. Addressing main challenges of miRNA analysis, a bioinformatics pipeline was developed that contains bona fide miRNA annotations from MirGeneDB, utilizes the quality control software miRTrace, applies physiologically meaningful cutoffs and accounts for contribution of cell-type specific miRNAs and host tissue effects.
ResultsTwo hundred-and-seventy-five miRNA sequencing datasets were analyzed, and after adjusting for the contribution of heterogeneity in cellular composition, strong signatures for primary and metastatic CRC were identified. The signature for primary CRC contained many previously known miRNAs with known functions. Deregulation of specific miRNAs was associated with individual metastatic sites, but the metastatic signatures contained overlapping miRNAs involved in key elements of the metastatic process, such as epithelial-to-mesenchymal transition and hypoxia. Notably, four of these miRNAs (MIR-8, MIR-10, MIR-375, MIR-210) belong to deeply conserved families present in many other organisms, triggering questions about their evolutionary functions and opportunities for experimental validation.
Conclusion:Applying a meticulous pipeline for the analysis of smallRNA-seq data, miRNA signatures for primary and metastatic CRC were identified, contributing novel insights into miRNA involvement in CRC metastatic evolution and site-specific metastatic adaptations. New datasets can easily be included in this publicly available pipeline to continuously improve the knowledge in the field.