The mitochondrial cytochrome oxidase (CO) genes are involved in complex IV of the electron transport system, and dysfunction of CO genes leads to several diseases.However, no work has been reported on the codon usage pattern of these genes. We used bioinformatic methods to analyze the compositional properties and the codon usage pattern of the COI, COII, and COIII genes in fishes, birds, and mammals to understand the similarities and dissimilarities of codon usage in these genes, which gave an insight into the molecular biology of these genes. The effective number of codons (ENC) value of genes was high in different species of fishes, birds and mammals, which indicates that the codon bias of CO genes was low and the ENC values were significantly different among fishes, birds, and mammals, as revealed from the t test. The overall guanine and cytosine (GC) content in fishes, birds, and mammals was lower than 50% in all genes, indicating that the genes were AT-rich and significantly different among fishes, birds, and mammals. The TCA codon was overrepresented in fishes, birds, and mammals for the COI gene, in birds and mammals for the COII gene, but it was not overrepresented in others. Only three codons, namely CTA, CGA, and AAA, were overrepresented in all three groups for the COI, COII, and COIII genes, repectively. From the neutrality plot in fishes, birds, and mammals, it was observed that the slopes of the regression lines (regression coefficients) in the COI, COII, and COIII genes were <0.5, suggesting that natural selection played a major role, whereas mutation pressure played a minor role.