Whole-transcriptome RNA sequencing reveals CeRNA regulatory network under long-term space composite stress in Rats

Abstract

To systematically evaluate the effect of simulated long-term spaceflight composite stress (LSCS) in hippocampus and gain more insights into the transcriptomic landscape and molecular mechanism, we performed whole-transcriptome sequencing based on the control group (Ctrl) and the simulated long-term spaceflight composite stress group (LSCS) from six hippocampus of rats. Subsequently, differential expression analysis was performed on the Ctrl and LSCS groups, followed by enrichment analysis and functional interaction prediction analysis to investigate gene-regulatory circuits in LSCS. In addition, competitive endogenous RNA (ceRNA) network was constructed to gain insights into genetic interaction. The result showed that 276 differentially expressed messenger RNAs (DEmRNAs), 139 differentially expressed long non-coding RNAs (DElncRNAs), 103 differentially expressed circular RNAs (DEcircRNAs), and 52 differentially expressed microRNAs (DEmiRNAs) were found in LSCS samples compared with the controls, which were then subjected to enrichment analysis of Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways to find potential functions. PI3K-Akt signaling pathway and MAPK signaling pathway may play fundamental roles in the pathogenesis of LSCS. A ceRNA network was constructed with the predicted 340 DE pairs, which revealed the interaction roles of 220 DEmiRNA-DEmRNA pairs, 76 DEmiRNA-DElncRNA pairs, and 44 DEmiRNA-DEcircRNA pairs. Further, Thrombospondins2 was found to be a key target among those ceRNAs. Overall, we conducted for the first time a full transcriptomic analysis of the response of hippocampus to the LSCS that involved a potential ceRNA network, thus providing a basis to study the underlying mechanism of the LSCS.