Codon use bias is an important characteristic in the process of genetic information transmission of species. With the publication of the genome of C. panzhihuaensis, it is of great significance to investigate the codon use bias for understanding the genetic evolution of this species and for molecular breeding. In this study, Perl language and CodonW 1.4.2 software was used to systematically analyze the coding gene sequences of the C. panzhihuaensis genome, obtain the characteristics of codon use, and identify the sources of variation affecting codon use. The results showed that GC3s content, GCall content, and ENC value were 46.75%, 47.67%, and 52.8 respectively. As the first seed plant, the GC3 content of Cycas is similar to that of most gymnosperms and most dicotyledons, but considerably different from that of monocotyledons.RSCU value was greater than 1, among which 21 codons ended in U or A. Enc-plot analysis, PR2-plot neutral plot analysis, and correlation coefficient analysis of Axis1 with GC3s and CAI revealed that the codon bias was mainly affected by natural selection, but also by mutation pressure and other factors. Thirty-one optimal codons of Cycas genes were screened. This research can provide a reference for the phylogeny and genome codon evolution of Cycas.