Abstract
BACKGROUND: Numerous circular RNAs (circRNAs) have been recently identified in porcine tissues and cell types. Nevertheless, their significance in porcine spleen development is yet unelucidated. Herein, we reported an extensive overlook of circRNA expression profile during spleen development in Meishan pigs. RESULTS: Overall, 39,641 circRNAs were identified from 6,914 host genes. Among them, many circRNAs are up- or down-regulated at different time points of pig spleen development. Using WGCNA analysis, we revealed two essential modules for protein-coding genes and circRNAs. Subsequent correlation analysis revealed 67 circRNAs/co-expressed genes that participated in immnue-associated networks. Furthermore, a competing endogenous RNA (ceRNA) network analysis of circRNAs revealed that 12 circRNAs modulated CD226, MBD2, SAMD3, SIT1, SRP14, SYTL3 gene expressions via acting as miRNA sponges. Moreover, the circRNA_21767/miR-202-3p axis regulated SIT1 expression in a ceRNA manner, which is critical for the immune-based regulation of spleen development in Meishan pigs. CONCLUSIONS: Overall, our results demonstrated that the circRNAs were differentially expressed during different stages of porcine spleen development, meanwhile the circRNAs interacted with immune-related genes in a ceRNA-based fashion. Moreover, we presented biomedical researchers with RNAseqTools, a user-friendly and powerful software for the visualization of transcriptome profile data.