Spinal cord injury (SCI) is a devastating neurological disease without effective treatment. To generate a comprehensive view of the mechanisms involved in SCI pathology, we applied RNA-Sequencing (RNA-Seq) technology to characterize the temporal changes in global gene expression after contusive SCI in mice. We sequenced tissue samples from acute and subacute phases (2 days and 7 days after injury) and systematically characterized the transcriptomes with the goal of identifying pathways and genes critical in SCI pathology. The top enriched functional categories include “inflammation response,” “neurological disease,” “cell death and survival” and “nervous system development.” The top enriched pathways include LXR/RXR Activation and Atherosclerosis Signaling, etc. Furthermore, we developed a systems-based analysis framework in order to identify key determinants in the global gene networks of the acute and sub-acute phases. Some candidate genes that we identified have been shown to play important roles in SCI, which demonstrates the validity of our approach. There are also many genes whose functions in SCI have not been well studied and can be further investigated by future experiments. We have also incorporated pharmacogenomic information into our analyses. Among the genes identified, the ones with existing drug information can be readily tested in SCI animal models. Therefore, in this study we have described an example of how global gene profiling can be translated to identifying genes of interest for functional tests in the future and generating new hypotheses. Additionally, the RNA-Seq enables splicing isoform identification and the estimation of expression levels, thus providing useful information for increasing the specificity of drug design and reducing potential side effect. In summary, these results provide a valuable reference data resource for a better understanding of the SCI process in the acute and sub-acute phases.